AD-A226  022 

RADC-TR-90-42,  Vo!  II  (of  two) 
In-House  Report 
March  1990 


PROCEEDINGS  OF  THE  1989  ANTENNA 
APPLICATIONS  SYMPOSIUM 


Paul  Mayes,  et  al 


DT1C 


ELECTE 
AUG  3  0 1990 


Sponsored  by 

Directorate  of  Electromagnetics 
Rome  Air  Development  Center 
Hanscom  Air  Force  Base 
Air  Force  Systems  Command 

and 


Electromagnetics  Laboratory 
Department  of  Electrical  and  Computer  Engineering 
College  of  Engineering 
University  of  Illinois 
Urbana-Champaign,  Illinois  61801 


i 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


Rome  Air  Develoment  Center 
Air  Force  Systems  Command 
Griffiss  Air  Force  Base,  NY  13441-5700 


This  report  has  been  reviewed  by  the  RADC  Public  Affairs  Division  (PA) 
a->d  is  releasable  to  the  National  Technical  Information  Service  (NTIS).  At 
NTIS  it  will  be  releasable  to  the  general  public,  including  foreign  nations. 

RADC-TR-90-42 ,  Volume  II  (of  two)  has  been  reviewed  and  is  approved  for 
publication. 


ROBERT  J.  MAILLOUX 

Chief,  Antennas  and  Components  Division 
Directorate  of  Electromagnetics 


APPROVED: 


JOHN  K.  SCHINDLER 
Director  of  Electromagnetics 


Directorate  of  Flans  and  Programs 


If  your  address  has  changed  or  if  you  wish  to  be  removed  from  the  RADC 
mailing  list,  or  if  the  addressee  is  no  longer  employed  by  your  organization, 
please  notify  RADC  (EEAS  )  Hanscom  AFB  MA  01731-5000.  This  will  assist  us  in 
maintaining,  a  current  mailing  list. 


Do  hot  return  copies  of  this  report  unless  contractual  obligations  or  notices 
on  a  specific  document  require  that  it  be  returned. 


REPORT  DOCUMENTATION  PAGE 


torm  Approved 
OMb  No  0/04-0188 


Public  fvpd'tmg  <0<  lhi>  collation  ol  .nt*  tt*  jU  n  v  Mui-sw  '  r*  *,i  i  .  ■>,>*!  •  <m  *  i  ,  u  .  i  <1.  i  ,i  ><  i  ol  k  j  j  t(  i  »,»( 

gjlhonng  jnd  momt-ii img  the  d4lo  nttfUed.  jno  wntfju  uu  j  uno  fi %»x  .*  og  im  .JkUiuiivimi  v-tijwinr  hu  u  ^  ngm  j  i>iu  b«iU>  tio’nr,  iu  yi  in,  /iher  *#*>x.-.4  yi  thi> 

cOlleXtlOn  ut  informal  <K  ilWudifT  J  luggoUvii*  It/I  uJk.  t(«  J  U<»*  h  I  j  .t  i.l  >n  ,Uu  Ml  j-.i-jI'  ,  k.v>.  Mum  ilk  !..f  .1.1  .Iduitiun  *>,  •  utol  ,  *uj  Lk.  |<>)  kHum> 

OdviWiighway.  Suite  .04.  Arlington,  VA  12202*4  JQ/,  ai.a  to  th»‘  0*tuo  )t  MjMugeinont  <ma  nuc,  ,t 1  >'ipir>>v‘cl  *ooi  u».i*  K.  tc«,uu/t4  ol«o)  «.niM>.jtjn  wt  ^ot»G) 


1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE 

March  1990 


4.  TITLE  AND  SUBTITLE 

Proceedings  of  the  1989  Antenna  Applications  Symposium 


J.  fttl'Gitr  TYPE  A.xiU  dates  COVLIUO 

In-house 


6.  AUTHOR(S) 

Paul  Mayes  et  al 


S.  FUNOInU  NUMBERS 

PE  62702F 

PR  4600  TA  14  WU  PD 


7.  PERFORMING  ORGANIZATION  NAML(S)  AND  AUGRtSS(ES) 

Rome  Air  Development  Center 
Hanscom  AFB,  MA  01731-5000 


Responsible  Individual:  John  Antonucci  (EEAS) 


9.  SPONSORING/MONITORING  AGfcNCY  NAMt(S)  AND  ADDKEbSiU) 


u.  I'i  IuOiuvIii^L  CrtCANlZ aTiON 
(UPCkT  NUMblK 

RADC-TR-90-42 
Volume  II 


lu.  :*i'UN  CORING.  MONITORING 
aGLnCY  IxlI'OkT  NUMBER 


11.  SUPPLEMENTARY  NOTES 

Volume  I  consists  of  pages  1  through  342;  Volume  II  consists  of  pages  343 
through  611 


12a.  distribution/availability  statement 


12b.  UulMUUTlUN  CODE 


Approved  for  public  release;  distribution  unlimited 


ABSTRACT  (Moxirnum  200  words) 


The  Proceedings  of  the  1989  Antenna  Applications  Symposium  is  a  collection  of 
state-of-the-art  papers  relating  to  phased  array  antennas,  multibeam  antennas, 
satellite  antennas,  microstrip  antennas,  reflector  antennas,  HP,  VHP,  UHF  and 
various  other  antennas. 


4T  SUBJECT  TERMS 


Antennas 

Satellite  antennas)  / 
Broadband  antennas ,  / 


crostripT* 

<r 

^Reflector  , 

CEF'' 

Lvhf  , 


cMultibeam  antennas 
"-Array  antennas,/  h 


IS.  NUMBER  OF  PAGES 


$.  NUMBER  ( 

276 


16.  PtxiCE  CODE 

?IL  - 


17.  SECURITY  CLASSIFICATION 

OF  REPORT 

IB  SCCUUITy  CLASSIFICATION 

Of  THIS  PAGE 

19.  SlCUKll  Y  CLmS j»  iv.aTiUN 

OF  ABSTRACT 

hj.  LIMITATION  OFABSfltAcr 

Unclassified 

Unclassified 

Unclassified 

SAR 

NSN  7b40-01-280-SbOU 


SUf.Ouro  rOfm  2uB  (Rev  2-BV) 

•  i  -«J  u<  Vld  / 1  i  Itt 


Accesion  For  1 

NTIS  CRA&I 
OTIC  TAB 
Unannounced 
Justification 

£ 

□ 

□ 

By 

Distribution/ 

Availability  Codes 

Dist 

Avail  and/or 

Specal 

A'l 

CONTENTS 

*  Keynote:  'Future  deirectica  in  Solid-State  Millimeter- 

Wave  Arrays,"  J.  Kinzel ,  General  Electric,  Electronics 
Laboratory,  Syracuse,  NY 

WEDNESDAY,  20  SEPTEMBER  1989 

ANTENNA  ELEMENTS  AND  CPSS 


1.  “Microstrip  Antennas  in  Composite  Layers,"  D.  R.  Jackson,  1 

N,  G.  Alexopoulos,  and  A.  A.  Oliner 

2.  "A  Compact  Multi-Polarized  Broadband  Antenna,"  Dean  A.  31 

Hofer,  Oren  B.  Kesler,  and  Lowell  L.  Loyet 

3.  "A  Simple  Circularly  Polarized  Primary  Antenna  for  Direct  54 

Broadcast  From  Satellite  Receivers,'  Takashi  Nakamura 
and  Kamel  Maamria 

4.  “Closed  Loop  Performance  of  an  Adaptive  Polarization  69 

Combiner,'  Kyle  A.  Dallabetta,  Bernard  J.  Lamberty, 
and  Kenneth  G.  Voyce 


#  5.  'A  Dual  Polarized  Horn  With  a  Scanning  Beam,"  Z.  Frank 

6.  "A  New  Circular  Polarization  Selective  Surface  (CPSS) ,'  85 

Gilbert  A.  Morin 


REFLECTOR  AND  LENS  ANTENNAS 

*  7.  "Design  and  Analysis  of  a  dual  Reflector  Multibeam  Torus 

Antenna,"  S,  Carrillo 

8.  'The  ARSR-4  Antenna,  a  Unique  Combination  of  Array  and  110 

Antenna  Technology,'  Richard  S.  White,  and  Daniel  Davis 


#  NOT  INCLUDED  IN  THIS  VOLUME 


i  i  i 


Contents 


I  ; 

9.  *A  Segmented  Mirror  Antenna  for  Radiometers,'  S.  W.  Lee,  124 

B.  Houshmand ,  •  M.  Zimmerman,  and  R.  Acosta 

10.  'Reflector  Performance  Degradation  Due  to  an  Arcjet  163 

Plume,'  H.  Ling,  G.  A.  Hallock,  H.  Kim,  B.  W.  Birkner, 
and  A.  Zaman 

11.  'Parametric  Estimation  of  the  J-Band  Radiation  Properties  185 

of  Cylindrical  Feeds  With  or  Without  External  Coaxial 
Cups,’  Constantinos  J.  Papathomas 

12.  'Impedance-Matched  coaxial-Probe  Feed  for  a  Rotman  Lens-  213 

Fed  Array  with  4.5:1  Bandwidth  and  Sector  .Pattern 
Coverage,'  Joseph  A.  Troychak,  Peter  S.  Simon 

13.  "Ku-Band  Antenna  for  VSAT  Applications,"  Nicholas  Moldovan  223 

THURSDAY,  21  SEPTEMBER  1989 
ARRAYS 

14.  ‘Electromagnetic  Coupler,'  Peter  R.  Franchi  and  239 

Nicholas  P.  Kernweis 

15.  .  'Design  of  a  Planar  Array  cf  Parasitic  Microstrip  Patch  251 

Antennas,'  R.  W.  Shaw,  J.  K.  Kovitz,  I.  Paz , 

L.  Johnson  and  G.  D.  Arndt 

16.  "Advances  in  Monolithic  Horn-Antenna  Imaging  Arrays,  271 

Philip  Stimson,  Yong  Guo,  Karen  Lee  and  David  Rutledge 

17.  "Wideband  dual  Polarized  Apertures  Utilizing  Closely  283 

Spaced  Printed  Circuit  Flared  Slot  Antenna  Elements 
for  Active  Transmit  and  Receive  Phased  Array 
Demonstration,"  Mark  J.  Povinelli 

18.  "Design,  Performance  Characterization  and  Hybrid  Finite  311 

Element  Boundary  Element  Analysis  of  a  Linearly 
Polarized  Printed  Circuit  Tapered  Notch  Array,' 

Mark  J.  Povinelli  and  John  D’Angelo 

19.  "EHF  Active  Array,"  Joe  Kobus  343 

20.  'Unique  Joint  Stars  Phased-Array  Antenna,'  H.  Shnitkin  360 

iv 


Contents 


ARRAYS  AND  NUMERICAL  TECHNIQUES 


21.  'An  Octive  Band  High  Gain  Receive  Cylindrical  Array 

Antenna  With  Full  Azimuth  Coverage,”  J.  C.  Herper , 

A.  M.  bucceri  and  J.  J.  Stangel 

22.  'Conformal  Active  Phased  Array  Demonstration , “ 

Jerome  D.  Hanfling 

23.  'High  Precision  Frequency  Locking  technique  for  Active 

Microstrip  Antenna  Arrays,'  Gabriel  Colef , 

George  Eichmann,  Raymand  L.  Camisa  and  Aly  E.  Fathy 

#  24.  ‘ EHF  Polyrod  Array  Antenna.'  P.  W.  Brock  and  R.  W.  Major 

25.  'PHEDRE:  A  Computer  Package  for  the  Study  of 

Electromagnetic  Radiation,"  Christian  Renard ,  Fabrice 
Clerk  and  Louis  Beaulieu 

*  26.  ‘A  Novel  Technique  for  the  Synthesis  of  Separable 

Aperture  Distributions  Over  Arbitrarily-Shaped  Array 
Surfaces,'  W.  W.  Milroy 

27.  ‘Moment  Method  Analysis  of  a  Str ipl ine-Fed  Finite  Slot 

Using  Subdomain  Basis  Functions,'  Peter  S.  Simon, 

Kim  Mclnturff,  Deron  L.  Johnson,  and  Joseph  A. 
Troychak 

28.  (No  Paper) 


381 

411 

441 

463 


FRIDAY,  22  SEPTEMBER  1989 
EXPERIMENT  AND  THEORY 

29.  "A  Technique  to  Measure  Large  Antenna  Arrays  That  Contain 

Random  Amplitude  Errors,'  Peter  R.  Franchi  and 
Harvey  Tobin 

30.  ‘Ah  Incident  Field  Sensor  for  EM?  Measurements,'  Everett  G. 

Farr  and-  Joseph  S.  Hofstra 

31.  'Quantifying  Simultaneous  V  and  H  Coverage  From  an 

Underside  Aircraft  Location,*  George  Monser  and 
Leonard  T.  Surette 

»  NOT  INCLUDED  IN  THIS  VOLUME 


507 

527 

549 


v 


Contents 


Polarization  Determination  in  a  DF  System  at 
Wavelengths  Without  a  Rotary  Joint,"  Harry  B 

The  Effect  of  Feed  on  the  Input  Impedance  of  a 
■.Antenna,"  W.  C.  Chew,  Z.  Nie  and  Y.  T.  Lo 


MM 

.  Sefton,  Jr 
Microstrip 


EHF  ACTIVE  ARRAY 


JOE  KOBUS 


SEPTEMBER  1989 


MOTOROLA  INC. 

Government  Electronics  Group 
Mail  Drop  R1216 
8220  E.  Roosevelt 
Scottsdale,  AZ  85251 


-1- 


343 


UNCLASSIFIED  ABSTRACT 


This  paper  describes  a  circularly  polarized,  minimum  volume, 

solid-state  active  aperture  for  tactical  applications  at  Q-band.  The  design 

consists  of  32  waveguide  modules  arrayed  into  four  groups,  each  group  fed 

fran  an  active  source.  Broadwall  Ciavin  compensated  slots  are  used  to 

♦ 

maximize  gain  and  reduce  mutual  coupling.  Five  Impatt  modules  (one  driver 
and  four  output)  are  incorporated  utilizing  spatial  combining.  Circular 
polarization  is  accomplished  via  a  meander-line  polarizer  suspended  in  front 
of  the  aperture.  The  array  has  50%  minimum  efficiency  over  a  4.5%  band. 
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1.0 


SUMMARY  OF  OBJECTIVES 


Primary  objectives  of  the  transmit  antenna  design  are: 

•  Integrated  Radio/Antenna  Assembly 

•  Compatibility  with  receive  array  in  a  single  assembly 

•  Capable  of  integration  with  four  active  sources 

•  Maximum  gain 

•  Minimum  volume 

•  Minimum  cost 

The  goal  was  a  440  cubic  inch  rectangular  box  containing  all  electronics  and 
antennas.  Although  other  antenna  types  were  considered,  the  most  convenient 
and  volume  efficient  are  planar  arrays  permanently  attached  to  a  face  of 
the  box.  In  addition  to  attitude  positioning,  setting-up  a  terminal  then 
only  requires  removal  of  a  protective  cover.  No  connections  need  be  made  or 
canplex  mechanical  structures  assembled.  An  integrated  array,  furthermore, 
is  obviously  compatible  with  a  rectangular  packaging  form  factor  appropriate 
for  minimum  volume.  To  realize  the  program  objectives,  Motorola  thus 
developed  the  Q-band  active  array  modular  concept  shown  in  Figure  1 ;  a 
photograph  of  the  feasibility  demonstration  hardware  is  presented  in 
Figure  2. 

The  transmit  array  occupies  nearly  half  of  approximately  a  one- 
foot  square  area;  a  K-band  receive  array  occupies  the  other  half  of  the  same 
surface.  To  achieve  maximum  aperture  efficiency,  a  waveguide  slot  array  was 
chosen  for  the  Q-band  antenna.  Minimum  overall  thickness  results  using 
broadwall  slots  grouped  into  32  identical  modules  "window-fed"  from  a 
waveguide  distribution  network.  Each  group  of  8  is  excited  from  one 
quadrant  of  the  antenna  distribution  network,  and  4  active  sources  deliver 
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UNCLASSIFIED 

ANTENNA/TRANSMITTER  MODULAR 

CONCEPT(U) 
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COUPLING  SLOTS  TO 
RADIATING  GUIDES 
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Fi  sure  2.  CU')  Q-Bancl  Array  Feasibility  Demonstration  Model  with  Polarizer  Removed 


paver  to  the  quadrants.  Circular  polarization  of  the  radiated  fields  is 

1  7 

produced  by  ?  rectangular  sandwich  ’  located  in  front  of  the  array. 

Regarding  the  Q-barri  transmit  array,  a  50  percent  overall  aperture 
efficiency  (including  all  feed  network  and  polarizer  losses)  was  achieved  by 
Motorola's  feasibility  model.  This  is  at  least  as  good  as  electrical 
efficiencies  commonly  obtained  with  reflector  antennas. 

Motorola's  plan  for  realizing  the  objectives  in  a  modular  design 
was  devised  to: 

•  Establish  radiation  module  and  feed  network  fabrication  techniques 

early  in  the  program 

•  Optimize  all  components  during  development 

•  Utilize  computer  modeling 

•  Scale  radiation  module  including  polarizer  to  X-band  to  establish 
design  parameters 

•  Integrate  antenna  and  transmitter  modules  in  levels 
These  goals  were  all  realized. 

2.0  BASIC  DESIGN  CONSIDERATIONS 

The  SATCOMA  transmit  array  gain  goal  is  37.5  dBI,  including  all 
losses.  Fifty  per  cent  minimum  efficiency  is  also  a  design  goal.  Thus  the 

3  4  5 

required  uniformly  weighted  aperture  gain  is  40.5  dBI.  A  Cl av in-compensated 
waveguide  slot  was  selected  for  minimum  element-to-element  coupling  and 
maximum  radiator  gain  (7.8  dBI ) .  Under  these  assumptions  the  number  of 
array  elements  required  becomes 


N  =  10 


1  10  ) 


=  1862 


(1) 
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Analysis  reveals,  however,  waveguide  slot-array  technology  limits  the 
maximum  number  of  elements  to  1536,  resulting  in  a  thinning  loss 


LT 


-10  log, - 
a10 


1536  1 

1862 
V  / 


.84  dB 


(2) 


Finally,  the  aperture  area  yielding  40.5  a5I  is 


A  = 


10 


,40.5/10 


62.81 


.  2 


an 


(3) 


at  44.5  GHz. 

Minimum  depth  is  realized  from  an  array  of  slots  formed  in  the 
broadwalls  of  a  number  of  waveguide  sections  or  "sticks".  The  size  of  each 
module,  or  group  of  sticks,  depends,  as  is  well-known,  on  the  bandwidth 
demanded  from  it.  Figures  3  and  4  shew,  respectively,  a  Q-barri  module  and 
its  scaled  X-band  counterpart. 

Refer  to  Figures  1  and  5.  Each  quadrant  of  the  array  consists  of 
8  modules  fed  from  a  waveguide  distribution  network  containing  7  magic-tee 
power  dividers.  An  active  signal  source  input  is  provided  to  each  quadrant. 
These  active  sources  are  driven  fran  a  single  input  to  another  waveguide 
distribution  network. 

In  order  to  demonstrate  effective  radiated  power  (ERP),  it  is 
first  desirable  to  ascertain  array  gain.  Invoking  the  reciprocity 
principle,  this  is  accomplished  through  operation  of  the  Q-band  transmit 
array  as  a  receive  antenna.  Active  sources  in  the  transmit  routing  netwonc 
are  removed  and  replaced  with  waveguide  sections.  Received  power  collected 


349 


-7- 


351 


fin  First  Q-Band  Motorola  Fabricated  Brassboard  Module 
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Figure  4.  (TD  Final  X-Band  Scaled  Radiation  Module 
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by  e3t.'h  of  the  four  quadrant  distribution  networks  is  combined  using  the 
transmitter  routing  network  (in  reverse)  sans  active  sources. 

Calculation  of  ERP  requires  knowledge  of  all  module  gain  and 
insertion  phase  variations,  as  well  as  the  quadrant  network  insertion  loss 
and  phase  variations,  together  with  the  power  and  phase  of  the  four  active 
sources.  This  information  was  determined  and  subsequently  compared  with 
recorded  ERP.  As  component  data  accumulated  during  the  engineering 
development  effort  budgets  were  prepared  and  updated.  The  last  of  these  is 
presented  in  Table  I. 

Our  approach  taken  to  size  the  aperture  is  based  on  40.5  dBLI 
ideal  gain,  together  with  a  50  percent  minimum  aperture  efficiency. 

Thinning  loss  of  .84  dB  was  found  to  result  for  a  practical  waveguide  slot 
array  design  consisting  of  32  modules  of  48  slots  each.  Thus  at  44.5  GHz, 
at  least,  individual  module  gains  should  be  nominally 

40.5  -  .84  -  10  log  (32)  =  24.6  dBI,  (4) 

excluding  impedance  mismatch  loss.  Recorded  module  gain  values  listed  in 
Table  I  include  mismatch  loss.  Consequently  we  may  conclude  radiation 
module  gain  is  consistent  with  our  overall  array  design  goal  of  37.5  dBI. 
3.0  FINAL  TEST  DATA 

Input- to-output  amplitude  and  phase  measurements  were  performed 
for  each  of  the  quadrant  distribution  networks  of  Figure  5,  as  well  as  for 
the  transmitter  routing  network  with  active  sources  removed.  These  data 
indicate  channel  line  losses  in  excess  of  those  expected  for  idea:! 
waveguides;  in  addition  to  magic  tees  nearly  1  dB  of  quadrant  network  loss 
appears.  Transmitter  routing  network  loss  is  not  relevant  to  array 
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Table  I.  (U)Final  Transmit  Array  Gain  Budget 


GAIN  OR  LOSS  PARAMETER 
(dB) 

FREQUENCY  F(GHz) 

43.5 

44.5 

:  45.5 

Average  Module  Gain  G  (dBLI) 

Including  Impedance 

Mismatch  Loss 

24.08 

24.35 

24.38 

Array  Factor  Gain 

10  log1Q  (32) 

15.05 

15.05 

15.05 

Module  Uniformity  Loss 
for  ±.25  dB  Amplitude 
and  ±10  Degree  Peak 

Phase  Deviations 
-10  log1Q  (nM) 

.06 

.06 

.06 

Combined  Loss  of  Four 

Quadrant  Networks 
-1°  !og10  („QN) 

1.30 

1.36 

i 

1.50 

Polarizer  Loss 

.35 

.35 

.35 

Axial  Ratio  Loss 

Relative  to  Perfect 

Circular  Polarization 

.037 

.024 

.021 

Expected  Gain  of 

Active  Aperture  (dBCPI) 

37.38 

37.61 

37.50 

Recorded  Passive  Array  Gain  (dBCPI) 

37.6 

37.6 

37.6 

354 


-12- 


performance  because  only  the  final  outputs  of  the  four  active  sources  affect 
ERP.  It  is  nevertheless  necessary  to  know  the  transmitter  network  less 
since  it  is  usta  as  a  passive  combiner  to  establish  array  gain.  Quadrant 
network  ohmic  losses,  on  the  other  hand,  do  impact  our  budget. 

Recorded  amplitude  and  phase  data  ODtained  for  the  quadrant 
networks  yield  the  composite  efficiency  of  the  antenna  distribution  networx 
(line  4  of  Table  I).  Loss  contributors  include  both  magic  tees  and  guide 
channels.  Analysis  reveals  phase  errors  are  so  small  as  to  produce  an 
almost  negligible  contribution  to  network  efficiency.  This  was  confirmed 
from  recorded  array  patterns  (e.g.  Figures  6  and  7)  which  generally  exhibit 
deep  nulls  between  sidelobes.  Table  II  summarizes  pattern  performance. 

Thus  it  maty  be  concluded: 

•  Magic  tees  produce  well-balanced  power  divisions  and  adequate 

isolation  of  radiation  module  mismatches 

•  All  waveguide  channels  are  precision  phased 
Machining  tolerances  must  therefore  be  on  target. 

In  order  to  perform  passive  array  tests,  active  sources  were 
removed  from  the  transmitter  routing  network,  and  replaced  with  four 
waveguide  feed-throughs  directly  to  the  quadrant  network  inputs.  These  are 
precision-phased  using  shims.  Amplitude  and  phase  data  recorded  for  the 
transmit  network  are  not  given  here,  but  were  incorporated  in  calculation  of 
passive  array  gain. 

Waveguide  covers  in  the  demonstration  model  are  attached  with 
screw- type  fasteners.  Lower  waveguide  channel  resistive  losses  can 
undoubtedly  be  obtained  by  laser  welding,  the  technique  proven  for  efficient 
radiation  module  fabrication.  Lower  waveguide  loss  can  conceivably  raise 

overall  transmit  array  efficiency  by  .8  dB  to  60  percent. 
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FIGURE  6  PASSIVE  MODE  H-PL ANE  PATTERN  FOR  Q-B AND  ARRAY 
INCLUDING  POLARIZER,  RECORDED  AT  44.5  GHz 
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UNIQUE 

JOINT  STARS  PHASED-ARRAY  ANTENNA 

H.  Shnitkin 
Norden  Systems,  Inc. 

Norwalk,  CT  06856 


1.  ABSTRACT 


The  Joint  STARS  radar  antenna  is  a  24-foot  wide  by  2-foot  high 
planar  X-band  slotted  array,  scanned  electronically  via  ferrite  phase 
shifters  to  +/-60  degrees  in  azimuth  and  scanned  mechanically  to  +/-100 
degrees  in  elevation.  (See  Figures  1  and  2.)  It  is  mounted  inside  a 
radome  underneath  the  fuselage  of  an  E-8A  aircraft,  (See  Figure  3),  and 
operates  in  a  mapping  as  well  as  in  a  three-port  interferometer  azimuth 
target  tracking  mode. 


Figure  1.  Joint  STARS  Antenna  on  Pattern  Range 
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Figure  3.  Antenna  Installed  in  EA-8  Aircraft 
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After  a  description  of  the  antenna's  architecture,  performance,  and 
construction,  its  five  unique  features  are  presented,  namely: 

1 .  Colinear  slotted  array  aperture 

2.  Switchable  dual-mode  aperture 

3.  Multifunction  electronic  phasing 

4.  Automated  phase  alignment 

5.  Triple  temperature  compensation 

2.  ANTENNA  DESCRIPTION 

The  Joint  STARS  antenna  is  a  horizontally  polarized,  rectilinear 
lattice  slot  array,  exhibiting  extremely  low  azimuth  and  elevation 
sidelobes,  and  is  capable  of  high  precision  azimuth  tracking.  Its  aperture 
illumination  is  switchable  via  ferrite  phase  shifters  so  that  it  can  generate 
a  single,  low  sidelobe  beam  or  three  triple  width  low  sidelobe  beams 
from  each  one-third  section  of  the  aperture.  The  latter  mode,  called 
interferometer  tracking,  essentially  divides  the  aperture  into  three  parts 
and  permits  individual  receive  phase  measurements  to  ascertain 
direction  of  arrival  from  differential  time  delay  or  phase  measurements. 
In  addition,  a  multichannel  microwave  receiver  with  Transmitter/Receiver 
limiters  and  RF  down  converters  is  part  of  the  antenna. 

2.1  Performance  -  Because  of  security  restrictions,  no  detailed 
performance  data  can  be  devulged,  except  that  the  antenna  operates 
over  3%  in  X-band,  exhibits  very  low  sidelobes  in  all  modes,  supports  a 
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high  power  radar  transmitter  and  possesses  fractional  milliradian 
boresight  and  beam  pointing  accuracies. 

2.2  Antenna  Architecture  -  The  aperture  consists  of  456  vertical 
linear  28-slot  array  radiating  elements,  which  are  connected  to  twelve 
waveguide  dual-mode  feeds  (DMF-i-38)  via  ferrite  phase  shifters, 
compensating  delay  lines,  and  twist  and  match  sections,  as  shown  in 
Figure  4.  Each  vertical  slotted  array  consists  of  a  four-way  divider  (EU-4) 
and  four  seven-slot  resonant  subarrays.  The  ferrite  phase  shifters  are  of 
the  latching,  nonreciprocal  type,  require  8-BIT  digital  phase  commands, 
and  can  switch  modes  or  scan  angles  within  7  ps  at  the  PRF.  The  twist 
and  match  sections  (quantity  of  1824)  are  rectangular,  horizontally 
polarized  waveguide  to  double-ridge,  vertically  polarized,  waveguide 
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transitions  with  a  high-Q  resonant  iris.  These  sections  interconnect  from 
the  aperture's  ridge  waveguide  to  the  four-way  power  divider’s 
rectangular  waveguide  and  provide  optimum  impedance  match  to  the 
array  at  42°  of  azimuth  scan. 

The  compensating  delay  lines  are  semirigid,  silicon  dioxide 
coaxial  cables  whose  group  delay  corrects  for  the  differential  phase 
delays  of  the  waveguide  dual  mode  feed  to  maintain  accurate  aperture 
phase  over  the  instantaneous  frequency  band.  Each  dual-mode  feed 
couples  38  radiators  to  either  the  sum  corporate  feed  in  the  mapping 
mode  or  to  one  of  three  interferometer  corporate  feeds  (AZ-*-4)  in  the 
tracking  mode.  Both  corporate  feeds  are  time  delay  compensated 
waveguide  directionalcoupler  feeds  and  therefore  achieve  broadband 
phase  and  amplitude  uniformity. 

In  addition,  a  loosely  coupled  traveling  wave  feed  is  added  to  the 
rear  of  the  aperture  with  a  -45  dB  coupling  path  to  each  of  the  456 
elevation  slot  arrays.  Its  purpose  is  to  simulate  a  received  wave  from 
about  47°  in  azimuth  by  injecting  an  RF  test  signal  into  the  BIT  INJECT 
port.  This  BIT  feed  is  used  for  aperture  phase  alignment  and  in-flight 
pattern  performance  checks. 

2.3  Construction  -  The  mechanical  construction  is  shown  in 
Figure  5.  A  four-sided  structural  frame  is  formed  by  a  three-sided 
aluminum  trestle  and  the  radiating  aperture.  The  large  triangular 
windows  in  the  trestle  surfaces  allow  easy  access  to  install  and  service 
all  electronic  components  while  assuring  structural  rigidity  and  aperture 
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flatness  within  .030  inch  TIR.  Semirigid  cables  were  chosen  to  connect 
the  phase  shifters  to  the  feed  so  as  to  allow  for  mechanical  alignment 
tolerances  and  to  facilitate  phase  shifter  replacement. 


3.  COLLINEAR  SLOTTED  ARRAY  APERTURE 

The  aperture  consists  of  456  vertical  28-slot  arrays,  each  a  single¬ 
ridge  waveguide  with  broadwall  shunt  slots.  Traditional  broadwall 
slotted  arrays  require  alternating  slot  displacement  from  the  waveguide 
centerline,  as  shown  in  Figure  6.  This  configuration  can  be  construed  as 


t  _ _ _  -.n  »»  A  »  r\|i  10  O  O  O  I  of  pH  "R11  c  Irvtc 
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Both  the  "A"  and  the  "B"  array  possess  grating  lobes  due  to  their  large 
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Figure  6.  Conventional  Waveguide  Broadwall  Slotted  Array 


vertical  slot  separation.  But  together  at  0°  azimuth  the  grating  lobes  are 
changed  into  ordinary  sidelobes  by  the  combination  array  factor.  This  is 
no  longer  true  during  azimuth  scan  causing  a  gradual  transition  from  side 
lobes  to  intercardinal  grating  lobes  as  scan  angle  is  increased.  [Ref.  1] 
Since  these  secondary  beams,  or  grating  lobes,  exceed  the  required 
sidelobe  level,  a  truly  collinear  slotted  array  radiator  was  devised.  [Ref.  2 
and  3] 

All  slots  are  placed  on  a  common  centerline,  while  the  trough 
depths  of  the  ridge  waveguide  are  altered  periodically  to  achieve  a 
waveguide  field  asymmetry,  as  illustrated  in  Figure  7.  Thus,  the  electrical 
centerline,  where  the  longitudinal  H-field  vanishes,  car.  be  displaced  by 
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Figure  7.  Collinear  Slotted  Array  Configuration 

varying  amounts  from  the  physical  centerline.  Slot  conductance,  guide 
wavelength,  and  07180,  phase  can  therefore  be  controlled  by 
appropriate  choice  of  the  ridge  waveguide  trough  asymmetry.  This 
results  in  an  easily  manufacturable  array  design  with  all  slots  on  a 
common  centerline. 

The  success  of  the  l  slot  collinear  array  was  demonstrated  by 
measured  array  elevation  patterns,  which  achieved  sidelobe  levels  to 
meet  the  specified  maximum  sidelobe  envelope. 
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4. 


SWITCHABLE  DUAL-MODE  APERTURE 

By  appropriate  commands  to  the  ferrite  phase  shifters,  the  twelve 


dual-mode  feeds  convert  one  total  low  sidelobe  aperture  with  tapered 
illumination  for  the  mapping  mode  into  three,  one-third  aperture,  tapered 
illuminations  for  the  tracking  mode  (see  Figure  8).  Each  dual  mode  feed 
consists  of  two  cascaded  traveling  wave  directional  coupler  distribution 
feeds  (Figure  9),  similar  to  a  two-beam  Blass-type  multibeam  feed.  [Ref. 
4  and  5].  However,  instead  of  the  similar  illumination  and  different  linear 
phase  tapers  of  the  Blass  feed,  two  different  illuminations  and  two  phase 
functions,  specifically  designed  for  maximum  decoupling  between  the 
two  traveling  wave  feeds,  are  produced.  Mode  switching  via  the  phase 
shifters  is  accomplished  by  selectively  adding  one  of  two  phase 
collimating  functions  to  the  scan  angle  phase  taper.  Each  phase 
collimating  function  produces  a  linear  aperture  phase  taper  for  only  one 
mode  and  a  decollimated,  and  thus  decoupled,  aperture  phase  for  the 
other. 


FORE 


( 'ENTER 


AFT 


V 


i  i 

INTERFEROMETER  •  FEED  AMPLITUDE  TAPER 

Figure  8.  Aperture  Illumination  Function 
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Figure  9.  Dual-Mode  Feed  (DMF+38) 

The  RF  delays  between  the  two  rows  of  couplers  have  been 
adjusted  at  the  factory  to  precisely  yield  the  required  output  port  power 
levels  for  the  rear  (interferometer)  feed  by  controlling  the  phase 
difference  between  the  voltage  coupled  from  the  rear  feed  and  that  from 
the  front  feed.  The  latter  is  produced  by  unavoidable  coupling  of  the  front 
feed  directional  couplers. 

While  the  dual-mode  feed  produces  the  two  prescribed  aperture 
magnitudes,  one  for  total  aperture  and  one  for  one-third  aperture,  the 
transmission  phases  to  each  radiator  are  allowed  to  be  arbitrary.  This 
flexibility  was  purposely  provided  to  allow  the  component  manufacturers 
to  concentrate  on  very  accurate  amplitude  control.  The  arbitrary 
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transmission  phases  are  subsequently  compensated  by  the  ferrite  phase 
shifters,  using  a  phase  alignment  table  as  described  in  Section  6. 

5.  MULTIFUNCTION  ELECTRONIC  PHASING 

5.1  Description  of  the  Ferrite  Phase  Shifters  -  The  Joint  STARS 
antenna  utilizes  228  latching,  nonreciprocal,  dual  ferrite  phase  shifter 
assemblies  manufactured  by  Electro-Magnetic  Sciences,  Inc.  (See 
Figure  10.)  Nineteen  of  these  are  serially  loaded  at  a  5-MHz  clock  rate 
with  8-Bit  phase  shift  commands.  A  redundant  set  of  data  shift  registers 
are  provided  to  avoid  catastrophic  antenna  failure  should  a  single  shift 
register  in  any  chain  of  nineteen  fail.  Single  saturated  current  reset  is 
used  to  erase  previous  magnetic  history.  A  1 .5-degree  RMS  differential 


Figure  10.  Ferrite  Phase  Shifter  Assembly 
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phase  accuracy  is  mainly  achieved  via  a  magnetic  flux  sampling  coil,  a 
volt  xtime  =  flux  integrator  circuit,  and  an  analog  flux  command 
comparator.  The  latter  also  performs  a  BIT  function  by  confirming  all 
successful  ferrite  flux  drive  sequences.  Another  important  contributor  to 
phase  accuracy  is  an  E-prom  containing  measured  microwave  phase 
data  at  specific  frequencies  and  temperatures.  This  data  is  used  to 
construct  phase  linearization  tables  needed  to  modify  input  phase 
commands  to  obtain  the  desired  RF  phase,  as  explained  in  Section  7.1. 

5.2  Phase  Shifter  Tasks  -  Because  of  the  high  cost  per  ferrite 
phase  shifter,  six  tasks  have  been  assigned  to  it  in  addition  to  the 
obvious  beam  scanning  function,  as  explained  below. 

1.  Mode  Switching  -  Either  mapping  beam  or  the  three 
interferometer  tracking  beams  are  selected  via  the  phase 
shifters.  Specifically,  the  dual-mode  feed  (DMF-r38)  phase 
delays  are  phase  compensated  for  aperture  collimation  for 
either  mapping  or  tracking  mode,  while  randomly  phased 
for  the  unwanted  mode.  Since  mode  switching  can  be 
accomplished  within  7  ps,  the  radar  can  transmit  a  full 
aperture  SUM  beam  while  receiving  on  three,  one-third 
aperture,  INTERFEROMETER  beams. 

2.  Duplexing  -  Adding  180°  to  the  commands  for  phase 
shifters  No.  1  through  228,  switches  the  aperture  collimation 
trom  the  Sum-Receive  port  of  the  AZ-sum  feed  to  the 
difference  or  Transmit  port.  (See  Figure  4.)  For  Sum- 
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Receive  the  E-plane  arm  of  the  input  magic  tee  is  used  and 
for  Transmit,  the  H-plane  arm.  Thus  a  circulator/duplexer 
handling  the  full  transmit  power  is  avoided. 

3.  Aperture  Phase  Alignment  -  The  combined  feed  assembly 
is  permitted  a  200°  random  phase  tolerance  for  ease  of 
manufacturing  and  for  minimizing  phase  shifter  quantization 
lobes.  Every  phase  shifter  command  contains  a  phase 
compensation  term  to  collimate  the  above  random  phase. 
This  compensation  phase  is  obtained  by  the  automated 
phase  alignment  procedure  described  in  Section  6  and  is 
stored  in  eight  alignment  tables,  for  four  frequencies  and 
two  modes. 

4.  Thermal  ( ompensation  -  Aperture  phase  errors  due  to 
thermal  distortion  (bow  and  expansion)  of  either  the 
aperture  or  the  corporate  feed  are  compensated  for 
temperature  via  corrective  phase  shifter  commands  to 
maintain  low  sidelobes  in  flight.  (For  details  see  Section  7). 

5.  Beam  Broadening  -  A  quadratic  or  a  cosine  phase  error  is 
added  to  the  phase  shifter  commands  to  broaden  the  main 
lobe  for  the  rapid  search  radar  mode  without  raising  the 
random  sidelobe  level.  The  amount  of  quadratic  or  cosine 
phase  error  can  be  adjusted  for  as  high  as  a  hundred-fold 
increase  in  beamwidth. 
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6.  Self-Test  -  Sequentially  increasing  the  aperture  phase 
taper  commands  to  the  ferrite  phase  shifters,  while  applying 
a  test  signal  to  the  BIT-Inject  port  (Figure  4),  produces  an 
electrically  scanned  antenna  pattern  at  the  radar  ■'eceiver. 
This  pattern  can  be  analyzed  for  conformance  to  the 
specification  and  for  antenna  trouble  shooting. 

Summary  -  The  total  resulting  phase  command  to  all  ferrite 
phase  shifters  consists  of  the  sum  total  of  the  following  four  456-number 
arrays: 

a.  Ideal  scan  angle  phase  taper 

b.  Insertion  phase  alignment  correction 

c.  Thermal  compensation  phase 

d.  Beam  broadening  phase 

6.  AUTOMATED  PHASE  ALIGNMENT 

6.1  Description  -  The  feed  system  exhibits  up  to  200°  of  random 
phase  error  due  to  manufacturing  tolerances  of  the  waveguide  corporate 
feeds,  dual-mode  feeds,  delay  cables,  ferrite  phase  shifters,  four-way 
power  dividers,  and  aperture.  To  maximize  antenna  gain,  minimize 
sidelobes,  and  maintain  boresight  accuracy,  phase  compensation  must 
be  supplied  through  the  ferrite  phase  shifters.  This  requires  a 
measurement  of  aperture  phase  error  from  which  an  insertion  phase 
alignment  table  is  generated.  Essentially,  this  phase  alignment  table  is 
defined  as  the  456  negated  measured  aperture  phase  values  of  the 
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assembled  antenna  system.  To  achieve  perfect  aperture  phase,  i.e. 
perfectly  linear  phase  taper,  always  add  the  above  negated  measured 
phase  values  to  their  respective  scan  phase  commands. 

The  radar  provides  a  computerized  method  for  phase  aligning  the 
aperture  in  less  than  half  a  second,  using  either  a  received  radiated  RF 
signal  or  a  test  signal  applied  to  the  BIT-Inject  port,  shown  in  Figure  4. 
(for  details  see  Ref.  6) 

Three  sequential  phase  alignments  are  generally  performed; 
these  are  defined  as  follows: 

1.  Range  Alignment  (RA):  This  is  performed  at  the  factory, 
using  a  far-field  radiated  signal  from  the  zero  azimuth 
direction,  at  temperature  profile  T1. 

2.  First  BIT  Calibration  ( BC1 ):  This  is  performed  at  the  factory 
using  an  RF  test  signal  applied  to  the  BIT-Inject  port  at 
temperature  profile  T1. 

3.  Second  BIT  Calibration  (BC2):  This  is  performed  in-flight, 
using  an  RF  test  signal  applied  to  the  BIT-Inject  port,  but  at 
ambient  temperature  profile  T2. 

Applying  the  negated  RA  values  as  phase  shifter  commands 
generates  constant  phase  across  the  aperture  at  temperature  T1  in  the 
Transmit  mode.  Subtracting  the  two  alignment  tables,  RA  minus  BC1, 
produces  a  BIT-Feed  Correction  table.  This  table  provides  for 
temperature  insensitive  phase  conversion,  allowing  a  BIT-Feed  aligned 
antenna  to  be  converted  into  an  aperture-aligned  antenna.  Applying  the 
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negated  RA  -  BC1  +  BC2  values  as  phase  shifter  commands  produces 
constant  aperture  phase  in-flight  at  temperature  profile  T2,  by  correcting 
RA,  the  alignment  at  temperature  II,  by  the  change  in  alignment  due  to 
in-flight  temperature  change  from  T1  to  T2. 

The  following  procedure  will  align  the  antenna  aperture: 

1.  Apply  RF  excitation  (far-field  or  BIT-injection) 

2.  Set  phase  shifters  to  456  sequential  linear  phase 
commands  [456  linear  phase  differentials,  from  0°  to 
(455/456)  x  360°].  This  is  equivalent  to  456  progressive 
scan  angles. 

3.  Measure  received  voltage  vectors  at  the  antenna's  sum  and 
three  interferometer  ports  for  each  of  the  above  456  phase 
commands. 

4.  Perform  a  Fast  Fourier  Transform  (FFT)  computation  on  the 
456  voltages  measured  in  (3)  to  obtain  complex  aperture 
excitation  voltages. 

5.  Negated  phase  of  excitation  voltages  of  (4)  constitute  an 
alignment  table  (RA,  BC1,  or  BC2)  for  addition  to  the  phase 
shifter  scan  commands. 

6.  To  verify  the  above  procedure  repeat  steps  (1)  through  (4) 
after  adding  the  alignment  table  of  (5)  to  all  phase 
commands.  Examine  the  measured  pattern  in  (3)  and  the 
phase  and  amplitude  deviations  in  (4)  to  ascertain  the 
alignment  quality.  A  3°  rms  phase  error  and  a  0.34  dB  rms 
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amplitude  error  will  assure  adequately  low  random 
sidelobe  level. 

6.2  Experimental  Results  -  Simulated  computer  print-outs  from 
the  above  alignment  procedures  are  shown  in  Figure  1 1  to  illustrate  the 
precision  of  the  aperture  amplitude  and  phase  alignment  as  well  as 
theelectronically  scanned  antenna  sidelobes.  The  upper  plot  starts  at  the 
left  with  the  peak  to  the  main  beam,  followed  by  all  the  right-side  and 
then  the  left-side  sidelobes,  and  ends  at  the  right  just  short  of  the  peak  of 
the  main  beam.  The  other  two  plots  are  generated  after  performing  an 
FFT  and  show  amplitude  and  phase  precision.  Precision  is  highest  near 
the  array  center  and  poorest  at  the  array  edges,  where  the  array  voltage 
is  only  about  9%  that  at  the  center. 

The  measured  patterns  of  the  five  completed  antenna  systems 
illustrate  the  quality  of  the  alignment.  In  all  cases  the  specified  maximum 
and  92%  sidelobe  limits  were  achieved  for  the  24-foot  array. 
Furthermore,  the  measured  patterns  are  consistent  with  sidelobe 
predictions  based  on  aperture  phase  and  amplitude  tolerances  as  well 
as  systematic  sidelobes. 

6.3  BIT-Feed  Diagnostics  -  As  a  by-product  of  BIT-Calibration  or 
Range  Alignment,  the  quality  of  the  antenna  pattern,  generated  by 
measuring  received  signal  strength  while  commanding  the  phase  shifters 
to  sequentially  increasing  aperture  phase  tapers,  and  the  causes  of 
pattern  deficiencies  can  be  diagnosed,  as  follows. 
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Examine  the  electronic  scan  pattern  which  was  measured  after 
applying  alignment  phase  correction.  If  the  sidelobes  are  excessive,  look 
for: 

1.  Exceptionally  low  amplitude  of  an  individual  radiator, 
indicating  a  faulty  RF  connector,  waveguide,  or  phase 
shifter. 

2.  Excessive  random  phase  variation,  indicating  a  faulty 
phase  shifter  or  a  noisy  receiver. 

3.  Excessive  random  amplitude  variation,  indicating  a  faulty 
corporate  feed,  a  faulty  dual-mode  feed,  or  a  noisy  receiver. 

7.  TRIPLE  TEMPERATURE  COMPENSATION 

7.1  Phase  Shifter  Phase  Linearity  -  Both,  a  temperature  sensor 
and  an  E-prom  containing  a  measured  RF  phase  data  table,  are  resident 
on  every  phase  shifter  for  the  purpose  of  instructing  the  beam  steering 
computer  how  to  calculate  the  appropriate  phase  linearization  command 
table.  This  table  is  used  to  look  up  the  compensated  or  linearized  phase 
command  to  the  phase  shifter  needed  to  achieve  the  desired  RF  phase 
shift  at  the  current  temperature  and  frequency. 

7.2  Aperture  Distortion  -  Data  from  temperature  sensors  at  the 
aperture  and  at  the  rear  of  the  antenna  structure  (trestle)  allow 
computation  of  phase  shifter  command  compensation  to  correct  for  a 
thermal  aperture  bow.  This  compensation  is  handled  just  like  the  beam 
broadening  described  in  Section  5.  In  addition,  thermal  aperture  growth 
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is  computed  and  differential  phase  between  adjacent  radiating  slots 
recalculated  at  the  new  temperature  before  scanning  the  beam 
electronically. 

7.3  Feed  Distortion  -  Ambient,  phase  shifter  dissipation,  and 
nonuniform  transmitter  power  heating  of  the  waveguide  corporate  feed 
cause  severe  aperture  phase  errors  in  flight.  These  are  corrected  by  the 
three  sequential  phase  alignment  procedures  described  in  Section  6, 
namely:  Range  Alignment,  BIT  Calibration  at  same  temperature  at  the 
factory,  and  a  second  in-flight  BIT  Calibration.  Every  several  minutes  an 
updating  in-flight  BIT  Calibration  (BC2)  is  performed.  Using  the  factory 
supplied  tables  of  RA  and  BC1,  an  updated  aperture  alignment  table  of  - 
(RA  -  BC1  +  BC2)  is  computed  and  added  to  each  phase  shifter 
command  every  nine  minutes.  (See  Section  6.) 

8.  CONCLUSION 

The  design,  performance,  and  unique  features  of  the  Joint  STARS 
phased  array  antenna  were  presented.  These  included: 

1 .  A  novel  collinear  slot  array  configuration 

2.  Switching  between  low  sidelobe  sum  and  low  sidelobe 
interferometer  modes  via  ferrite  phase  shifters 

3.  A  rapid  in-flight  phase  alignment  technique 

4.  A  method  of  maintaining  phase  alignment  and  low 
sidelobes  over  a  wide  temperature  range. 
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AN  OCTAVE  BAND  HIGH  GAIN  RECEIVE  CYLINDRICAL 
ARRAY  ANTENNA  WITH  FULL  AZIMUTH  COVER  E* 
by  J.C.  Herper,  A.  M.  Bucceri  and  J.  J.  Stangel 
UNISYS  CORPORATION 
UNISYS  S&GSG 
GREAT  NECK,  NEW  YORK 

Abstract  -  An  octave  bandwith,  compact  high  gain  receive  antenna 
providing  multiple  beams  with  360-degree  azimuth  coverage  has 
been  developed.  A  test  antenna  has  been  fabricated  consisting  of 
a  Luneburg  lens  collimator  feeding  a  cylindrical  array  of  128 
column  networks. 

Multiple  simultaneous  beams  are  realized  in  azimuth  by 
focusing  the  received  signals  through  the  Luneburg  lens  to  a 
focal  aperture.  A  circulator  at  each  lens  port  provides  the 
isolation  between  the  column  networks  and  the  receiver.  Adjacent 
lens  ports  are  combined  after  the  circulators  to  provide  the 
desired  focal  aperture  to  which  a  receiver  is  connected.  The 
lens  focal  aperture  design  controls  antenna  beamwidth  to  be 
approximately  constant  with  frequency  for  an  octave  bandwidth. 

The  column  networks  generate  a  CSC^  fan  beam  in  elevation. 
Each  printed  stripline  circuit  consists  of  a  broadband  16  to  1 
power  combiner  with  baluns  and  dipoles  incorporated  into  a  com- 

*  This  work  was  supported  in  part  by  the  USAF,  Air  Force  Systems 


Command,  Hanscom  AFB  under  contract  number  F19628-85-0018 . 


pact  and  lightweight  design*  Broadband  Vilkenson  tee  type  power 
combiners  and  transmission  line  lengths  are  configured  to  realize 
the  appropriate  amplitude  and  phase  distribution.  Schiffman  type 
phase  shifters  constitute  the  balun  design  providing  the  trans¬ 
formation  from  unbalanced  to  balanced  two  wire  line  feed  for  the 
wideband  dipole. 

Full  360  degree  azimuth  coverage  is  demonstrated  for  simul¬ 
taneous  beams  showing  agreement  with  computer  model  predictions 
for  both  beamwidth  and  sidelobe  level  control.  These  patterns 
also  confirm  the  relationship  between  the  Luneburg  lens  focal 
aperture  and  antenna  beamwidth  invariance  with  frequency. 


INTRODUCTION 


This  paper  describes  the  design,  development  and  experimen¬ 
tal  evaluation  of  an  octave  bandwidth  cylindrical  array  for  an 
advanced  radar  application.  The  application  requires  that  the 
array  efficiently  receive  signals  simultaneously  in  a  cluster  of 
contiguous  high  gain  beams  which  may  be  electronically  scanned  to 
any  azimuth  sector.  The  combined  requirement  for  high 
efficiency,  high  cross-over,  low  sidelobes  and  octave  bandwidth 
represents  a  challenging  antenna  design  problem. 

The  design  approach  is  illustrated  by  the  functional  diagram 
of  Figure  1.  N  column  networks  with  integrated  radiating 
elements  are  arranged  to  form  the  cylindrical  array.  Signals 
incident  on  the  array  exit  from  the  column  networks  through 
single  pole,  double  throw  (SPDT)  switches  to  a  two-dimensional 
Luneburg  lens  network.  The  Luneburg  lens  has  a  circular  geometry 
with  ports  along  its  perimeter.  Signals  incident  from  a  given 
azimuth  are  focused  toward  the  opposite  ports  of  the  lens.  The 
focused  signals  exit  through  the  SPDT  switches  to  the  electronic 
switching  matrix.  The  switching  matrix  accesses  the  M  Luneburg 
lens  ports  associated  with  the  azimuth  sector  designated  for  beam 
formation.  The  switching  matrix  contains  single  pole,  multi¬ 
throw  switches  where  the  number  of  throws  equals  N/M  and 
determines  the  angular  width  of  the  sector  for  simultaneous 
coverage . 
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The  signals  emanating  from  the  switch  ports  are  amplified 
and  combined  in  a  beam  weighting  network  to  produce  K  beams.  To 
achieve  low  sidelobes  and  high  cross-over  between  adjacent  beams, 
a  lossy  beam  weighting  network  is  generally  needed,  and  hence  the 
low  noise  amplifiers  (LNA’s)  preserve  signal-to-noise  ratio. 
With  proper  choice  of  beam  weights,  the  antenna  produces  beams 
whose  gain,  beamwidth,  sidelobes  and  pointing  direction  are 
substantially  independent  of  frequency. 

These  design  principles  and  performance  characteristics  are 
demonstrated  by  the  experimental  model  antenna  pictured  in  Figure 
2.  The  model,  which  operates  in  S-through  C-band,  is  a 
cylindrical  array  formed  from  128  columns  of  vertically  polarized 
dipoles.  The  unit  is  capable  of  simultaneously  producing  eight 
contiguous  beams  over  any  electronically  selectable  45  degree 
angular  sector.  As  illustrated  in  the  block  diagram  of  Figure  3, 
the  configuration  and  operation  of  the  demonstration  antenna  is 
similar  to  that  of  the  baseline  design  approach  except  that  (1) 
circulators  are  used  in  place  of  SPDT  switches,  (2)  3  dB  power 

dividers  are  used  after  the  circulators  to  simplify  the  switching 
matrix  and  eliminate  the  beam  weighting  network  and  (3)  the  LNA’s 
are  deleted.  The  design,  construction  and  measured  performance 
of  the  demonstration  model  are  detailed  in  the  following  sub¬ 
sections  of  this  paper. 

DEMONSTRATION  MODEL  DESIGN 

The  cylindrical  array/Luneburg  lens  antenna  represents  a 
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SEAM  PORTS 
TO  RECEIVERS 

Figure  1  -  CYLINDRICAL  ARBAY  DESIGN  APPROACH  -  A  two-dimensional 
Luneburg  lens  beamformer  allows  wideband  simultaneous  reception 
on  a  cluster  of  directive  beams  covering  an  electronically 
selectable  azimuth  sector. 


Figure  2  -  DEMONSTRATION  MODEL  -  The  key  principles  of  the 
approach  are  verified  in  an  experimental  model  which  operates 
over  the  octave  band  from  3.0  to  6.0  GHz. 
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near  optimal  design  for  an  octave  band,  high  gain  receive  antenna 
with  360  degree  coverage  in  azimuth.  This  configuration  requires 
25  percent  fewer  elements,  is  more  compact  and  utilizes  a  single 
lens  beamformer  compared  to  a  multifaced  planar  array  approach. 
The  Luneburg  lens  coupled  with  a  multiple  beamforming  switching 
network  is  a  simpler  and  more  reliable  design  than  the 
complicated  commutator  transfer  switch  commonly  used  in  cylin¬ 
drical  arrays. 

The  cylindrical  array  lattice  design  is  based  on  the  scan 
requirements  determined  by  elevation  coverage  for  the  cosecant- 
squared  fan  beam  and  the  projected  illuminated  aperture  on  the 
cylinder  in  azimuth.  For  octave  bandwidth  operation,  the  rect¬ 
angular  lattice  inter-element  spacing  was  calculated  at  the 
highest  operating  frequency  to  conform  to  a  main  beam  elevation 
coverage  of  25  degrees  and  an  azimuthal  scan  of  £60  degrees.  The 
demonstration  model  antenna  is  comprised  of  a  cylindrical  array 
of  2,048  vertically  polarized  dipoles  arranged  in  128  columns. 
The  array  is  40  inches  in  diameter  and  22-1/2  inches  high. 
Dipoles  are  integrated  with  a  stripline  elevation  distribution 
network  to  comprise  a  single  column  network  assembly.  A 
corporate  power  divider  with  wideband  Wilkinson  tee  couplers 
generates  a  vertical  illumination  consistent  with  realizing 
cosecant  squared  elevation  coverage  to  25  degrees  over  the 
specified  octave  band  (3.0  to  6.0  GHz}. 

Azimuth  collimation  of  signals  entering  the  antenna  through 
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the  column  networks  is  implemented  with  a  two-dimensional 
Luneburg  lens1,  whose  interconnection  with  the  column  networks  is 
illustrated  in  Figure  3.  The  Luneburg  lens  shows  superior 
performance  characteristics  over  an  octave  bandwidth  because  of 
the  continuity  of  its  index  of  refraction.  Simultaneous  colli- 
mation  in  all  azimuth  directions  is  achieved  by  radially  varying 
the  propagation  constant  in  the  lens  medium.  This  is  accom¬ 
plished  by  varying  the  effective  dielectric  constant  to  provide 
real  time  delay  for  wideband  operation. 

Signals  entering  the  Luneburg  lens  along  a  given  arc  come  to 
a  focus  at  a  point  opposite  to  this  arc.  Thus,  the  same  port  of 
the  Luneburg  lens  is  an  entry  point  of  signals  from  the  column 
network  and  an  exit  point  for  signals  to  the  switch  matrix  where 
these  signals  impinge  on  the  cylindrical  array  from  substantially 
opposite  azimuth  directions.  Adjacent  focal  ports  may  be 
weighted  and  combined  to  provide  the  appropriate  lens  excitation 
to  meet  antenna  beamwidth  and  sidelobe  requirements.  In  general, 
for  beamwidth  invariance  control,  focal  aperture  lengths  greater 
than  one  wavelength  will  contract  azimuth  illumination  with 
increasing  frequency  thereby  reducing  the  variation  in  beamwidth. 
The  insertion  of  wideband  circulators  between  the  column  networks 
and  Luneburg  lens  ports  isolates  exiting  and  entering  signals  by 
virtue  of  the  circulators  non-reciprocal  properties.  Use  of  the 
circulators  enables  the  use  of  a  single  beamformer  lens  and 
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results  in  a  more  compact  beam  feed  network  for  octave  band 
operation. 

LUNEBURG  LENS 

The  design  of  the  Luneburg  lens  is  based  on  a  radially 
symmetric  lens  whose  dielectric  constant  is  determined  as  a 
function  of  the  radius  to  any  point  within  the  lens 
circumference.  The  relationship  governing  this  principal  was 
derived  by  R.  K.  Luneburg^  using  geometric  optics: 

n  - Je 7  »  Jz  -  r2/R2  0  <  r  <  R 

where  n  is  the  refractive  index,  R  is  the  radius  to  the  perimeter 
of  the  lens  and  r  is  the  radius  to  any  point  within  the  lens 
region.  The  amount  of  dielectric  fill  in  a  region  is  directly 
analogous  to  the  refractive  index  of  that  region.  By  varying  the 
height  of  the  dielectric  as  a  function  of  the  radius ,  radiation 
of  a  point  source  at  the  edge  of  the  lens  produces  a  linear  phase 
front  exiting  on  the  diametrically  opposite  side  of  the  lens  as 
illustrated  by  the  ray  trace  diagram  of  Figure  4.  A  computer 
program  was  implemented  to  provide  lens  design  and  milling 
machine  data  for  the  fabrication  of  a  piece-wise  linear, 
contoured  lens. 

The  general  function  for  this  contour  is  derived  from 
Snell’s  law  and  ray  tracing  methods (3)  showing  that  at  a  point  on 
a  ray  path: 

JTV  (r/R)sin^»  sin  0O 

where  cf  is  the  angle  between  the  ray  path  and  the  radial  line  and 
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Figure  3  -  DEMONSTRATION  MODEL  BLOCK  DIAGRAM  -  The  model  is  an 
electronic  scanning  cylindrical  array  incorporating  all  the  key 
component  hardware  components  inherent  to  the  design  approach. 


Figure  4  -  LUNEBURG  LENS  BEAMFORMER  COLLIMATING  PRINCIPLE  -  The 
radially  symmetric  varying  refractive  medium  in  the  two-dimen¬ 
sional  lens  focuses  plane  waves  incident  from  any  azimuth 
angle. 
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60  is  the  entering  ray  angle. 

The  demonstration  Luneburg  lens  has  a  lid  and  a  housing  as 
shown  in  Figure  5.  Together  they  form  a  cylindrical  cavity, 
41.16  inches  in  diameter  and  0.50  inch  high.  This  height  is  such 
that  all  but  the  dominant  parallel  plate  waveguide  mode  are  cut¬ 
off.  Inside  the  cavity  a  Noryl  lens  is  permanently  bonded  to  the 
lid.  This  convex  lens  is  40.00  inches  in  diameter  and  0.40  inch 
high  at  the  center.  The  cavity  formed  by  the  lid  and  housing 
protects  the  dielectric  lens.  SMA  connectors  are  mounted  on  the 
outside  of  the  housing.  The  center  conductor  with  Teflon 
insulation  penetrates  through  the  plate  into  the  cavity  where 
probe  elements  are  affixed  by  means  of  a  small  setscrew.  There 
are  128  probe  elements  mounted  circumferentially  along  the  outer 
perimeter  of  the  cavity.  The  probes  are  equally  spaced  at  one- 
half  wavelength  for  the  high  frequency. 

The  design  of  the  Luneburg  lens  probes  employs  a  grating 
lobe  series  analysis  for  the  determination  of  the  input  imped¬ 
ance.  This  assumes  that  the  curvature  of  the  lens  is  large 
enough  for  the  probes  to  appear  as  though  they  are  in  a  linear 
array.  From  this  an  initial  probe  configuration  is  obtained. 
Using  a  planar  array  simulator  test  fixture,  the  probe  design  is 
adjusted  for  a  good  match.  Figure  6  shows  the  final  probe  match 
measured  in  the  actual  Luneburg  lens  environment  and  indicates 
better  than  a  1.6  VSWR  over  the  octave  band. 

The  excitation  of  the  probes  can  be  designed  to  produce  far 
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Figure  5  -  DEMONSTRATION  MODEL  LUNEBURG  LENS  -  The  Luneburg  lens 
used  in  the  demonstration  model  is  a  cylindrical  cavity  whose 
internal  refractive  index  variation  is  determined  by  a  contoured 
dielectric  sheet. 


FREQUENCY  (GHZ) 


Figure  6  -  LUNEBURG  LENS  PROBE  MATCH  -  The  input  impedance  of  the 
Luneburg  lens  is  matched  over  an  octave  bandwidth  with  a  maximum 
VSWR  less  than  1.5. 
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field  patterns  whose  azimuth  beamwidths  remain  approximately 
constant  over  the  octave  bandwidth.  The  normal  variation  of 
bearawidth  for  a  fixed  aperture  over  an  octave  bandwidth  is  2:1. 
Consequently,  the  effective  aperture  of  the  antenna  must  vary 
with  frequency  to  maintain  beamwidth  invariance.  As  illustrated 
in  Figure  7,  the  effective  aperture  produced  by  the  Luneburg  lens 
changes  with  frequency  to  minimize  beamwidth  variation.  If  the 
probes  are  excited  to  form  a  line  source  of  length  L,  the  width 
of  the  effective  aperture,  W,  approximately  equals  2R  / L  where 
is  the  wavelength  and  R  is  the  lens  radius.  This  results  in  a 
far  field  half-power  beamswidth  of  approximately  L/2R,  which  is 
independent  of  frequency.  If  the  probe  excitation  is  such  that  L 
is  greater  than  a  wavelength  at  the  lowest  operating  frequency, 
the  far  field  beamwidth  is  substantially  constant  over  the  band. 
CYLINDRICAL  ARRAY  ANALYSIS 

The  performance  of  the  cylindrical  array /Luneburg  lens 
antenna  is  simulated  using  a  computer  program  developed  for  this 
specific  antenna  configuration.  The  simulation  computes  the  far 
field  patterns  via  a  diffraction  calculation  using  the  cylindri¬ 
cal  array  amplitude  and  phase  distribution.  The  far  field  pat¬ 
tern  can  be  calculated  using  Luneburg  lens  voltages  obtained  from 
a  theoretical  model  or  from  experimental  data. 

The  Luneburg  lens  is  designed  to  give  the  desired  far  field 
pattern  performance  by  optimizing  the  focal  illumination.  The 
appropriate  output  ports  of  the  Luneburg  lens  are  assigned  volt- 
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age  phases  and  amplitudes  by  superposing  probe  excitation  on  the 
opposite  side  of  the  lens.  These  probes  are  in  the  focal  region 
and  are  spaced  and  electrically  weighted  to  control  the  focal 
illumination.  The  focal  illumination  includes  the  effects  of 
mutual  coupling  by  employing  element  patterns  for  the  probes 
computed  using  modal  analysis.  For  this  case  a  planar  dipole 
array  modal  analysis  program  approximates  the  local  region  of  the 
probes  with  the  probe  spacing  used  as  the  column  spacing.  The 
losses  through  the  inhomogeneous  lens  cavity  are  included  in  the 
amplitude  determinations. 

Computed  patterns  for  the  demonstration  antenna  are  pres¬ 
ented  in  Figure  8*  The  antenna  simulation  is  for  64  pairs  of 
equally  weighted  focal  probes  forming  64  beam  ports.  From  these 
results,  the  demonstration  model  Luneburg  lens  was  designed, 
fabricated  and  tested.  Far  field  patterns  generated  from  mea¬ 
sured  lens  port  voltages  are  presented  in  Figure  9.  There  is 
good  agreement  between  the  theoretical  and  experimental  patterns 
with  regard  to  3dB  beamwidth  and  near-in  sidelobes  as  shown  over 
the  entire  octave  frequency  band.  High  backlobes  present  in  the 
generated  patterns  for  the  measured  lens  result  from  inter-probe 
coupling5  which  causes  energy  to  be  reflected  back  through  the 
circulators  and  to  the  array  face.  These  backlobes  can  be  elimi¬ 
nated  by  replacing  the  circulators  with  SPDT  switches. 

For  a  202  beamwidth  variation  over  an  octave  band,  an  opti¬ 
mum  design  utilizing  three  lens  probes  in  the  focal  aperture  with 
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Figure  7  -  DESIGN  FOR  CONSTANT  BEAMWTDTH  -  The  Luneburg  lens 
facilitates  an  antenna  design  whose  azimuth  beamwidth  remains 
substantially  invariant  with  frequency. 


LOW  FREQUENCY 
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Figure  8  -  COMPUTED  DEMONSTRATION  MODEL  PATTERNS  -  The  azimuth 
patterns  of  the  cylindrical  array  are  predicted  using  a  computer 
model  which  incorporates  analytically  derived  coupling  coeffi¬ 
cients  for  the  Luneburg- lens . 
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a  power  weighting  of  (0.5,  1.0,  0.5)  is  implemented  as  shown  in 
Figure  10,  This  design  requires  low  noise  amplifiers  in  conjunc¬ 
tion  with  these  lossy  feed  networks  since  beam  orthogonality 
requirements  for  all  passive  simultaneous  beam  antennas  limit  the 
efficiency  for  multiple  probe  feed  networks.  For  this  design 
receive  signals  pass  through  circulators  to  the  switch  matrix 
which  selects  a  sector  consisting  of  16  contiguous  lens  ports  and 
are  amplified  before  being  combined  through  the  beamforming 
network  resulting  in  7  contiguous  beamports. 

A  summary  of  the  predicted  performance  generated  from  both 
computed  and  measured  lens  port  voltages  is  given  in  Figure  11 
which  illustrates  the  3  dB  beamwidth  control  and  the  peak 
sidelobe  level  for  various  lens  feed  configurations  over  the 
octave  band.  These  generated  far  field  characteristics  confirm 
the  lens  feed  aperture  relationship  derived  in  Figure  7. 

ELEVATION  COLUMN  NETWORK 

A  column  network  illumination  is  synthesized  to  realize  a 
cosecant-squared  fan  beam  to  25  degrees  in  elevation  at  the  high 
frequency.  The  complex  voltages  at  the  elements  are  determined 
by  Woodward  synthesis  using  complex  beam  weights  to  match  the 
desired  elevation  beam  pattern.  The  voltage  distribution  trans¬ 
forms  into  a  fan  beam  in  the  far  field  with  the  peak  2.5  degrees 
above  the  horizon  and  -3  dB  or  more  energy  on  the  horizon  over 
the  entire  frequency  band.  The  highest  sidelobe  is  at  5  degrees 
below  the  horizon  and  15  dB  below  the  beam  peak  at  the  high 
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Figure  9  -  COMPUTED  PATTERNS  WITH  MEASURED  COUPLING  DATA  -  Mutual 
coupling  in  the  Luneburg  lens  causes  the  backlobes  of  the  demon¬ 
stration  model  to  increase  by  5  to  12dB.  This  effect  is  virtu¬ 
ally  eliminated  when  SPDT  switches,  rather  than 
circulators,  interconnect  the  Luneburg  lens  and  the  column  net¬ 
work. 


1  2  •  •  •  •  7 


Figure  10  -  BASELINE  FEED  NETWORK  SCHEMATIC  -  A  feed  network 
which  uses  three  Luneburg  lens  probes  for  each  beam  with  a  power 
weighting  of  (0.5, 1.0, 0.5)  reduces  sidelobes  and  improves  fre¬ 
quency  invariance. 

-  14  - 


396 


THEORY 


Figure  11  -  LUNEBURG  LENS  EXCITATION  TRADE-OFFS  -  Comparing  the 
azimuth  beamwidths  and  peak  sidelobes  produced  when  one,  two  or 
three  Luneburg  lens  probes  are  excited  for  beam  formation  shows 
that  three-probe  excitation  limits  beamwidth  changes  to  202  with 
lower  sidelobes  over  the  octave  band* 


frequency.  Figure  12  shows  the  cosecant-squared  fan  beam  cover¬ 
age  to  25  degrees  at  the  high  frequency  and  50  degrees  at  the  low 
frequency.  Ripple  compared  to  the  ideal  cosecant  squared  power 
drop  is  better  than  0.5  dB.  Element  amplitudes  are  modified 
slightly  to  simplify  coupler  design.  Taper  efficiency  ranges 
from  -1.8  dB  at  the  high  frequency  to  -1.1  dB  at  the  low 
frequency. 

The  desired  illumination  taper  is  generated  by  a  power 
divider  network  which  feeds  sixteen  vertically  polarized 
dipoles.  The  power  divider,  dipoles  and  wideband  baluns  are 
integrated  into  a  single  stripline  assembly  pictured  in 
Figure  13. 


ELEVATION  ANGLE  (DEGREES) 


Figure  12  -  CALCULATED  ELEVATION  PATTERNS  -  The  antenna  provides 
cosecant-squared  elevation  coverage  to  greater  than  25  degrees 
over  the  octave  band. 
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•  DISTRIBUTION 

DIPOLE  BALUN  NETWORK 


Figure  13  -  COLUMN  NETWORK  PRINTED  CIRCUIT  -  Sixteen  dipoles  and 
baluns,  and  a  wide  band  distribution  network  are  printed  on  a 
single  layer  simplifying  the  column  network. 
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The  power  divider  network  consists  of  15  compensated 
Wilkinson  tees  ranging  in  coupling  values  from  3.0  dB  to  4.5  dB. 
The  various  coupling  values  are  determined  by  a  top-down  power 
coupling  technique  from  the  element  ports.  The  VSWR  over  most  of 
the  frequency  band  was  measured  to  be  better  than  1.2  with  1.35 
maximum  at  the  band  edges.  These  values  are  very  close  to  the 
optimum  theoretical  values.  The  measured  loss  through  the  column 
network  is  then  1.8  dB  over  the  frequency  band. 

The  wideband  balun  is  comprised  of  three  cascaded 
60°  Schiffman  phase  shifters  and  a  transformer  to  140  ohm  bal¬ 
anced  two  wire  line.  This  design  allows  for  a  planar  topography 
whereby  networks  are  printed  on  a  single  sheet  resulting  in  a 
mechanically  simple  structure.  The  phase  shift  in  the  Schiffman 
phase  shifter  is  obtained  by  a  quarter-wave  resonant  coupled 
line.  The  measured  phase  error  and  VSWR  for  a  single  Schiffman 
stage  are  plotted  in  Figure  14  indicating  a  maximum  balun  phase 
imbalance  of  10  degrees  over  an  octave  band. 

The  radiating  element  is  a  simple  dipole  with  a  reactive 
wedge.  Starting  from  calculations  of  active  dipole  impedances 
obtained  from  a  cylindrical  array  modal  analysis^  computer  pro¬ 
gram,  a  design  for  the  launching  (two  wire  line)  section  was 
generated  matching  these  load  impedances  into  the  balun  section. 
A  single  element  waveguide  simulator  was  fabricated  and  the 
dipole  design  was  validated  and  refined  empirically  to  optimize 
the  match  over  the  octave  band.  An  average  VSWR  of  less  than 
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1,4:1  was  measured  over  the  operating  band  of  the  simulator. 

With  the  dipole  design  optimized,  element  patterns  for  the  single 
element  were  taken  in  both  principal  planes  at  sample  frequencies 
over  the  octave  band.  These  measurements  show  clean  smooth 
patterns  with  no  detected  resonances. 

The  integrated  stripline  network  is  printed  on  a  (KAPTON) 
sheet  as  previously  shown  in  Figure  13,  which  is  sandwiched 
between  two  foam  spacer  sheets.  Aluminum  spacers  0.125"  in 
height  are  used  to  maintain  a  continuous  spacing  between  the 
outer  aluminum  ground  planes.  The  printed  circuit  has  50  to 
70  ohm  stepped  transformers  at  the  connectors  to  match  into  the 
70  ohm  system  of  the  16:1  power  divider.  This  entire  network 
provides  a  light  weight  compact  design  which  measures  23.5  inches 
high  by  9  inches  wide  by  0.25  inch  thick. 

ANTENNA  TEST 

Antenna  pattern  data  on  the  composite  cylindrical  array 
assembly  demonstrate  key  performance  characteristics  over  an 
octave  band.  These  include:  360°  azimuthal  coverage,  relative 
azimuthal  beamwidth  invariance  and  sidelobe  level  control,  and  an 
elevation  fan  beam  with  CSC^  power  drop-off. 

The  cylindrical  array /Luneburg  lens  demonstration  antenna  is 
configured  with  64  azimuthal  beamoorts  which  can  be  accessed  one 
at  a  time  using  the  remote  beam  selector.  In  Figure  15,  each 
main  beam  is  recorded  sequentially  and  referenced  through 
360  degrees  of  azimuth.  The  data  shows  less  than  +0.5dB  varia- 
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Figure  14  -  60-DEGREE  SCHIFFMAN  PHASE  SHIFTER  PERFORMANCE  -  Each 
balun  contains  three  60-degree  Schiffman  phase  shifters  to  pro¬ 
vide  a  balanced  two-wire  feed  for  the  dipoles.  Measured  data 
shows  that  the  VSWR  is  less  than  1.2  and  the  phase  error  is  less 
tnan  3  degrees  for  each  Schiffman  stage  over  the  octave  band. 
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Figure  15  -  DEMONSTRATION  MODEL  AZIMUTH  COVERAGE  -  The  measured 
beam  patterns  of  the  demonstration  model  verify  uniform  360 
degrees  azimuth  coverage  within  +0.5dB  with  64  beams. 
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tion  in  gain  over  the  full  azimuth  coverage. 

To  demonstrate  beamwidth  control,  beamport  No.  13  is 
configured  for  a  triple  probe  excitation  with  (0.5,  1.0,  0.5) 
power  weighting.  As  shown  in  Figure  16,  62  beamports  are  con¬ 
figured  for  a  double  probe  excitation  and  beamport  No.  14  pro¬ 
vides  single  probe  excitation.  The  3  dB  beamwidth  variation  with 
frequency  for  the  measured  and  computed  far  field  patterns  for 
the  demonstration  antenna  are  plotted  in  Figure  17  for  the  vari¬ 
ous  lens  feed  configurations.  This  plot  shows  that  the  3  dB 
beamwidth  for  the  measured  and  computed  patterns  correlate  wel' 
over  the  octave  band. 

Since  the  3  dB  beamwidth  is  inversely  proportional  to  fre¬ 
quency  for  a  fixed  antenna  aperture,  the  beamwidth  control  can 
only  be  achieved  by  changes  in  antenna  aperture  distribution. 
For  the  cylindrical  array /Luneburg  lens  this  can  only  be  effected 
by  controlling  the  lens  aperture  distribution  via  a  combination 
of  the  focal  aperture  length  and  element  illumination.  The  probe 
element  spacing  in  the  Luneburg  lens  is  one  half  wavelength  at 
6.0  GHz.  When  two  lens  probes  are  fed  simultaneously,  the  focal 
aperture  is  increased  and  therefore  some  3  dB  beamwidth  control 
is  achieved  near  the  high  end  of  the  band.  When  three  lens 
probes  are  fed  simultaneously,  with  a  (0.5,  1.0,  0.5)  power 
weighting,  less  than  20  percent  variation  in  the  3  dB  beamwidth 
is  achieved  over  more  than  the  octave  band,  demonstrating  the 
desired  control. 
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Figure  16  -  LENS  PORT  FEED  MODIFICATION  -  To  experimentally 
verify  trade-offs  involving  azimuth  sidelobes  and  beamwidth 
invariance,  the  feeds  for  beams  13  and  14  on  the  demonstration 
model  are  modified  to  allow  single-probe  and  weighted  three-probe 
excitation  of  the  Luneburg  lens. 


THEORETICAL 
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Figure  17  -  DEMONSTRATION  MODEL  BEAMWIDTHS  -  The  measured 
beamwidth  for  three  different  Luneburg  lens  probe  excitations 
substantiates  that  beamwidth  variations  over  the  octave  band  are 
less  for  the  weighted  three-probe  excitation  as  expected. 


Recorded  azimuthal  and  cross-polarized  patterns  are  illus¬ 
trated  in  Figure  18  for  the  beamport  configured  with  the  triple 
probe  feed.  Sidelobe  levels  are  plotted  over  the  octave  band  in 
«  Figure  19  for  two  lens  feed  configurations. 

The  predicted  backlobes  from  inter-probe  coupling  in  the 
Luneburg  lens  is  also  observed  and  can  be  easily  improved  by 
20  dB  with  the  replacement  of  the  circulators  with  single  pole 
double  throw  (SPDT)  switches.  Experiments  were  conducted  to 
simulate  this  approach  over  a  wide  frequency  band,  and  Figure  20 
demonstrates  the  effectiveness  of  this  design. 

The  superimposed  cross-polarized  patterns  are  normalized  to 
the  peak  of  the  corresponding  co-polarization  patterns  recorded 
on  the  same  chart.  The  patterns  shown  are  in  the  azimuth  prin¬ 
cipal  plane.  In  all  cases,  the  cross-polarized  levels  are  about 
20  dB  below  the  peak  of  the  co-polarized  pattern  main  beam.  This 
result  indicates  that  the  demonstration  model  antenna  wideband 
element  achieves  good  element  radiation  characteristics  while 
minimizing  element  cross-polarization  level  for  more  than  an 
octave  operating  band.  Measured  elevation  patterns  are  in  good 
agreement  with  the  computed  patterns  over  the  octave  band.  Both 
the  computed  and  measured  patterns  are  illustrated  in  Figure  21 
for  comparison  at  the  design  frequency  of  6.0  GHz.  The  maximum 
deviation  from  a  cosecant  squared  power  drop  off  is  approximately 
1.0  dB  for  the  measured  elevation  pattern  at  6.0  GHz  compared  to 
0.5  dB  for  the  ideal  case.  These  small  deviations  from  the 
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Figure  18  -  DEMONSTRATION  MODEL  AZIMUTH  PATTERNS  -  Measured  co¬ 
polarized  and  cross-polarized  patterns  show  good  agreement  with 
theory.  The  cross-polarized  patterns  are  typica’.ly  more  than  20 
dB  below  the  co-poiarized . 
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Figure  19  -  DEMONSTRATION  MODEL  SIDELOBES  -  Measured  peak  and 
average  azimuth  sidelobes  for  two-probe  and  weighted  three-probe 
Luneburg  lens  excitations  indicates  that  the  sidelobes  for  the 
three-probe  excitation  are  about  4dB  lower. 
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Figure  20  -  AZIMUTH  PATTERN  WITH  SPDT  SWITCHES  -  The  effect  of 
using  SPDT  switches,  rather  than  circulators,  is  evaluated  on  the 
demonstration  model  by  terminating  beam  ports  in  the  vicinity  of 
the  ports  excited  for  the  desired  beam.  Comparison  of  the  resul¬ 
tant  pattern  at  4.5  GHz  with  that  of  Figure  18  shows  that  SPDT 
switches  reduce  backlobes  by  25  to  30dB. 
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ideal  patterns  for  this  broadband  network  can  be  attributed  to 
the  fabrication  tolerance  errors  and  grating  lobe  incipience. 
CONCLUSIONS  AND  SUMMARY 

An  octave  bandwidth,  electronic  scanning  cylindrical  array 
has  been  designed,  and  a  demonstration  unit  has  been  built  and 
tested.  The  array  receives  signals  simultaneously  on  a  cluster 
of  contiguous  directive  beams.  Using  a  novel  two-dimensional 


FIELD  ANGLE  (DEGREES) 

Figure  21  -  DEMONSTRATION  MODEL  ELEVATION  PATTERN  -  The  measured 
elevation  pattern  at  6.0  GHz  matches  the  desired  cosecant-squared 
coverage  to  25  degrees  within  +  1  dB. 
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Luneburg  lens  beamformer  with  a  simple  switching  network  con¬ 
taining  single  pole,  multi- throw  switches,  the  beam  cluster  is 
positioned  electronically  to  cover  any  azimuth  sector.  By  prop¬ 
erly  weighting  the  Luneburg  lens  excitation,  low  azimuth 
sidelobes  are  obtained,  and  beam  shapes  and  pointing  angles 
remain  substantially  invariant  over  the  octave  frequency  band. 
The  radiating  elements  are  axial  dipoles  which  are  integrated 
with  wideband  baluns  and  power  distribution  networks  into  single¬ 
layer  stripline  column  networks.  These  column  networks  produce 
cosecant  -squared  shaped  elevation  coverage  to  25  degrees  or  more 
in  elevation.  A  versatile  and  comprehensive  computer  simulation 
addressing  all  critical  aspects  of  the  design  has  been  generated 
and  validated  by  experimental  data  on  the  demonstration  model* 


-  27  - 


409 


REFERENCES : 

(1)  PEELER,  G.D.M.  and  ARCHER,  D.H.  (1953)  A  Two-Dimensional 
Microwave  Luneburg  Lens.  I.R.E.  Transactions  Antennas  and 
Propagation,  July  1953,  pp.  12-23. 

(2)  LUNEBURG,  R.K.  (1944)  Mathematical  Theory  of  Optics,  Brown 
University  Graduate  School,  Providence,  R.l. 

(3)  MATHIS,  H.F.  (1960)  Checking  Design  of  Stepped  Luneburg 
Lens .  IRE  Trans.  Antennas  and  Prop.,  May  1960, 
pp.  342-343. 

(4)  JASIK,  H.  (1954)  The  Electromagnetic  Theory  of  the  Luneburg 
Lens .  Report  Number  AF  CRC-TR-54-121 . 

(5)  HERPER,  J.C.,  HESSEL,  A.  and  TOMASIC,  B.  (1985)  Element 
Patterns  of  an  Axial  Dipole  in  a  Cylindrical  Phased  Array, 
Part  I:  Theory.  IEEE  Transactions  Antennas  and  Propagation, 
AP-33 ,  No.  3,  pp.  259-272. 

(6)  HESSEL,  A.  and  SHAPIRA,  J.  (1977)  Mutual  Coupling  Coeffi¬ 
cients  in  Collector  Arrays  on  Circular  Cylindrical  Concave 
Surfaces .  IEEE  Transactions  Antennas  and  Propagation,  AP-25, 
No.  6,  pp.  741-747. 


-  28  - 


410 


CONFORMAL  ACTIVE  PHASED  ARRAY  DEMONSTRATION 


Jerome  D.  Hanfling 

Raytheon  Company,  Equipment  Division 
Wayland,  MA 

Abstract 

The  development  and  demonstration  of  a  conformal  active  phased 
array  has  been  sponsored  by  the  Electromagnetics  Directorate  of 
Rome  A.ir  Development  Center  at  Hanscom,  AFB.  The  array  has 
transmit/receive  modules  which  utilize  state-of-the-art  gallium 
arsenide  monolithic  microwave  integrated  circuit  technology  for 
amplifiers  and  phase  shifters. 

The  array  aperture  is  linear  with  64-microstrip  patch  radiators 
flush  mounted  conformal  to  a  cylindrical  surface.  The  patch 
radiators  and  T/R  modules  combine  to  provide  a  low  profile, 
broadband  phased-array  element  for  low-noise  reception  and  high- 
power  transmit  operation  at  S-band  frequencies. 

The  control  of  the  active  devices  is  enabled  by  a  unique  beam 
steering  controller  VLSI  gate  array  resident  in  each  T/R  module. 
This  digital  controller  performs  phase  computations,  controls  the 
RF  paths  through  the  amplifiers  and  provides  built-in  test 
features  for  the  array. 

After  integration  and  checkout  of  the  array  system  indoors  in  the 
laboratory,  the  array  receive  and  transmit  performance  was 
demonstrated  by  outdoor  far-field  antenna  pattern  measurements. 
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1.0  INTRODUCTION 


A  conformal  antenna  experiment  called  "Array/Module  Integration 
and  Test"  (AMIT)  has  been  developed  for  the  Air  Force  Electronic 
Systems  Division  at  Hanscom  AFB,  under  Contract  No.  F19628-84-C- 
0137.  The  antenna  is  a  linear  array  of  64 -microstrip  patch 
radiators  flush  mounted  conformal  to  a  cylindrical  surface.  This 
surface  forms  the  front  of  a  rectangular  enclosure,  as  shown  in 
the  photograph  of  Figure  1 .  The  array  contains  Transmit/Receive 
(T/R)  modules  furnished  by  the  government  out  of  the  Conformal 
Array  Module  (CAM)  program  sponsored  by  The  Signal  Generation  and 
Control  Branch  of  the  Surveillance  Directorate  of  RADC  at 
Griffiss,  AFB  [1] .  The  patch  radiators  and  T/R  modules  combine 
to  provide  a  low-profile,  broadband  phased-array  element.  For 
the  AMIT  program,  the  patch  radiators  are  mounted  external  to  the 
skin  of  the  conformal  array  enclosure  and  the  T/R  modules  are 
fastened  on  the  inside,  as  shown  in  Figure  2.  Much  of  the  array 
system,  as  represented  in  the  block  diagram  of  Figure  3,  is 
contained  within  this  enclosure.  The  hardware  in  the  enclosure 
is  shown  in  Figure  4;  here  the  array  is  undergoing  bench-top 
tests  prior  to  outdoor  installation. 

This  paper  first  describes  the  array  system,  the  hardware  details 
and  the  command  and  control  system.  Then  the  low-profile 
microstrip-patch  radiator  is  discussed.  The  remainder  of  the 
paper  concentrates  on  some  unique  aspects  of  active  array  testing 
and  corresponding  test  results . 


2.0  ARRAY  HARDWARE  DESCRIPTION 


2 . 1  Array  System 

A  block  diagram  of  the  active  array  system  is  shown  in  Figure  3. 
The  heart  of  the  array  system  is  the  CAM  T/R  module,  as  shown  in 
Figure  5,  which  consists  of  an  RF  assembly  (containing  a  phase 
shifter,  and  a  power  amplifier  and  a  low-noise  amplifier,  see 
Figure  6) ,  and  a  Distributed  Module  Controller  (DMC)  gate  array 
separated  by  a  logic  interface  circuit.  Several  subsystems  are 
necessary  to  operate  the  active  array.  Inside  the  conformal 
array  enclosure  is  an  RF  feed  network,  an  ethylene  glycol /water 
cooling  manifold,  power  sxapplies,  energy  storage  capacitors  and  a 
junction  box  which  contains  the  array  interconnections.  Outside 
the  enclosure  are  the  beam  steering  controller,  HP-1000  computer, 
interlock  chassis,  and  monitors.  The  data  and  DC  connections  to 
the  modules  are  made  by  64  separate  cables  and  corresponding  15 
pin  D-sub  connectors.  The  array  system  is  powered  by  eight  10 
volt  power  supplies  (8  T/R  modules  per  supply) ,  one  5  volt  supply 
and  one  -10  volt  supply. 

The  RF  feed  network,  as  presented  in  Figure  7,  is  an  equal  power 
split  1:64  divider  network.  It  is  configured  from  commercial 
couplers  and  dividers.  The  sum  and  difference  amplitude  tapers, 
transmit  uniform  illumination,  and  compensation  for  different 
insertion  losses  >  the  network  through  the  T/R  modules  are 
obtained  by  a  variety  of  fixed  attenuators .  The  pads  are 
inserted  between  the  dividers  and  the  T/R  modules .  The  physical 
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layout  of  the  junction  box,  cables,  power  supplies  and 
interconnections  are  pictured  in  Figure  4 . 

2 . 2  Patch  Radiator 

The  patch  radiator  being  used  in  the  AMIT  program  is  balanced  fed 
providing  good  quality  horizontal  linear  polarization  (along  roll 
axis  of  cylinder)  over  a  wide  bandwidth  at  S-band  [2] .  Th^ 
element  spacing  has  been  selected  to  meet  a  +60  deg.  azimuth  and 
+40  deg.  elevation  scan  coverage.  The  goal  is  to  achieve  a  less 
than  2:1  VSWR  over  the  band  in  a  low  profile,  lightweight,  and 
low  loss  implementation. 

In  accordance  with  the  above  objectives  a  high  performance, 
compact  microstrip  patch  radiator  has  been  developed  [3]  [4] .  A 

mechanical  radiator  design  was  formulated  that  placed  the  double- 
tuned  matching  network/balun  structure  underneath  the  patch  and 
resulted  in  a  low  profile  structure  with  a  single  coaxial  input, 
as  shown  in  Figure  8 . 

The  basis  of  the  development  effort  for  this  Duroid  substrate 
patch  radiator  is  a  recently  developed  software  design  tool  that 
provides  a  numerical  code  that  accurately  models  patch  radiators 
in  an  infinite  array.  The  code  models  patch  radiators  mounted  on 
a  substrate  of  dielectric  constant,  e ,  each  patch  being  driven  by 
up  to  four  probe  feeds.  The  effect  of  the  matching  network  is 
not  included  in  this  code.  However,  the  code  allows  a 
determination  of  the  scattering  matrix  of  the  patch,  each  probe 
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being  considered  an  input  port.  Using  the  methodology  of 
references  5  and  6,  this  scattering  matrix  can  be  coupled  to  that 
of  the  matching  network  so  as  to  arrive  at  the  radiator's 
behavior  at  each  scan  angle  selected  [7] . 

The  design  was  carried  out  using  infinite  array  waveguide 
simulators .  These  measurements  and  the  above  code  were  combined 
to  predict  and  verify  performance.  The  results  are  summarized  in 
the  Smith  chart  plot  of  active  impedance  in  Figure  9 .  Both  the 
E,  H,  and  diagonal  (D)  plane  performance  are  close  to  ideal,  with 
the  exception  of  the  E-plane  performance  for  0  >50  deg.  and  for 
the  upper  end  of  the  frequency  band.  This  blind  spot  is  due  to 
the  generation  of  a  TM  surface  wave.  It  can  be  eliminated  by 
reducing  the  element  spacing  in  the  array  lattice.  A  better  than 
2:1  VSWR  out  to  45  deg.  scan  angle  has  been  achieved.  This  is 
just  beyond  the  required  40  deg.  scan  coverage. 

Sixty-four  patch  radiators  were  built,  tested,  and  installed  on  a 
10  foot  radius  of  curvature  section  representing  the  fuselage  of 
the  aircraft,  as  shown  in  Figure  1.  The  embedded  element 
patterns  were  measured  using  this  conformal  array  fixture  for 
various  element  locations  on  the  cylinder.  The  test  results,  as 
plotted  in  Figure  10,  compare  favorably  to  the  cos^/^q  and  cos  0 
ideal  gain  fall-off  curves. 
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3.0  ARRAY  COMMAND/ CONTROL  SYSTEM 


Each  monolithic  T/R  module  in  the  array  requires  a  digital  Beam 
Steering  Unit  (BSU)  to  perform  such  functions  as  phase 
computations  and  to  set  the  proper  mode  of  operation  of  the 
module.  A  distributed  approach  is  used  for  AMIT,  where  serial 
data  busses  send  to  every  DMC  the  same  steering  angle  and 
frequency.  These  are  then  clocked  into  each  T/R  module 
corresponding  to  row  and  column  intersections,  as  shown  in  Figure 
11,  where  the  phase  shift  setting  are  calculated  in  the  VLSI  DMC 
chip.  The  module  controller  accomplishes  the  steering  function 
by  computing  the  phase  shift  for  each  element  and  then  commanding 
the  phase  shifter  at  the  corresponding  location  to  the  computed 
value . 

In  addition  to  beam  steering,  the  module  controller  must  also 
control  both  the  Power  Amplifier  (PA)  and  the  Low-Noise  Amplifier 
(LNA) .  For  example,  during  the  transmit  mode,  as  shown  by  the 
switch  positions  in  Figure  6,  the  power  amplifier  is  energized 
and  the  T/R  switch  connects  the  PA  input  to  the  phase  shifter  and 
the  output  to  the  antenna.  At  this  time,  the  LNA  is  isolated 
from  the  transmit  amplifier. 

3 „ 1  Array  Operation 

Control  of  the  array  is  achieved  by  sending  serial  message 
signals  to  the  DMC  chips  in  response  to  the  commands  from  the 
host  computer.  The  RP68000  single -board  '  -  s'  in  the  beam 
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steering  controller  unit  executes  the  software  that  generates  the 
waveforms  to  program  and  control  the  DMC  chips  in  response  to  the 
commands  from  the  host  computer,  as  described  in  the  block 
diagram  of  Figure  12 .  The  DMC  interface  card  is  a  special 
purpose  peripheral  to  the  RP68000  board.  It  contains  shift 
registers  for  transmitting  MODE  and  DATA  words  and  for  receiving 
BITE  messages.  It  performs  clock  generation  both  for  the  shift 
registers  and  the  DMC  CLOCK  lines  and  it  performs  the  array 
element  addressing  function,  allowing  control  of  a  single  element 
or  simultaneous  control  of  all  the  T/R  modules  in  the  array. 

3 . 2  Steering  Calculation 

The  module  controller  accomplishes  the  steering  function  by 
computing  the  phase  shift  for  each  element  and  then  commanding 
the  phase  shifter  located  in  each  antennule  to  the  computer 
value.  The  phase  shift  that  must  be  applied  to  each  element  in 
order  to  properly  form  the  radiated  beam  is  (in  radians)  given  by 
equations  (1)  and  (2)  as  follows: 

2n  2  jt  2jt 

<j>  =  -  Y-l  sinp  +  -  Zj_  cosp  +  A4»i  -  -  Z^  (1) 

A<{> j_  =  Determined  from  calibration 
Zj_  =  Path  length  factor 

"  (  (2n  -  1)  As  ^  N 

Zj,  =  R  1  -  cos -  (2) 

l  2R  J  J 
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where 

R  =  radius  of  cylinder 
As  =  element  arc  length  spacing 

n  =  element  position  from  center  of  array 

A  digital  adaption  of  these  equations  is  used  in  the  BSU 
software . 

3 . 3  Self  Test 

The  controller  has  a  self  test  capability  that  can  be  used  to 
check  if  it  is  correctly  performing  the  beam  steering  computation 
as  well  as  its  other  control  functions .  Steering  data  (with  a 
known  result)  can  be  applied  to  the  controller.  After  the 
controller' s  phase  shift  is  computed,  its  value  is  returned  to 
the  command  and  control  interface  where  a  comparison  is  made  of 
this  result  with  the  known  result  and  a  pass/fail  determination 
is  made.  The  pass/fail  result  is  then  sent  to  the  main  display 
computer  in  the  form  of  a  status  word.  Control  data,  such  as  T/R 
switching  can  also  be  sent  to  the  controller  and  returned  for 
comparison.  The  controller  can  be  used  to  perform  RF  tests  on 
any  individual  T/R  module;  and  active  bit-wiggle  test  of  a 
specified  T/R  module  can  be  performed  by  sending  a  test  phase 
shift  word  from  the  command  and  control  interface  to  the  module 
control] er.  A  horn  external  to  the  array  is  then  used  to  verify 
the  proper  RF  operation  of  the  tested  module. 


The  68000  microprocessor  assembly  language  program  for 
controlling  an  array  of  elements  containing  DMCs  contains 
subroutines  for  the  various  array  modes  and  functions  as  follows: 
initialization,  diagnostics,  resets,  clears,  loads,  phase 
compute,  BITE,  receive  and  transmit  operation,  pulse  train 
protection  and  module  test  routines . 

4.0  ARRAY  CALIBRATION 

4 . 1  Calibration  Techniques 

After  the  checkout  of  the  array  in  the  laboratory  the  radiators 
were  installed  on  the  surface  of  the  conformal  housing.  The 
array  was  then  hoisted  onto  an  outdoor  pattern  test  pedestal . 

The  setup  for  operating  the  array  is  shown  in  Figure  13  and 
includes  various  equipments  for  mount  control,  array  electronic 
control,  DC  power,  air  conditioning,  heat  exchange,  RF  excitation 
and  reception,  signal  analysis,  temperature  and  current 
monitoring,  power  interlocks,  beam  steering  and  host  computer 
control . 

The  apriori  knowledge,  from  bench-top  tests,  of  the  insertion 
amplitude  and  phase  from  the  input  of  the  feed  to  the  output  at 
each  module  was  not  adequate  to  form  a  beam.  The  addition  of  "he 
radiators  and  the  outdoor  environment  necessitated  the  use  of  a 
free-space  calibration  technique. 
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A  calibration  procedure  was  developed  which  involved  receiving  a 
wave  front  from  a  far  field  source  and  sampling  the  amplitude  and 
phase  at  each  element.  The  key  to  the  technique  is  to  be  able  to 
turn  on  one  module  at  a  time  while  all  other  modules  are  turned 
off  and  all  radiators  are  switched  into  loads  internal  to  the  T/R  * 

module.  This  enables  the  insertion  amplitude  and  phase  to  be 
determined  in  the  mutual  coupling  environment  of  neighboring  ' 

radiators,  i.e.  the  actual  embedded  element  pattern  is  present. 

Once  this  signal  path  is  established  the  mount  can  be  rotated 
producing  an  embedded  element  pattern.  This  will  enable  the 
transmission  amplitude  and  phase  to  be  observed  for  all  angles  of 
interest.  Any  irregularities  in  the  element  pattern  may  indicate 
a  fault  in  the  microstrip  radiator  or  at  the  array  surface.  This 
would  have  to  be  corrected  before  the  calibration  can  proceed. 

Typical  good  element  patterns  for  radiators  at  locations  are 
shown  in  Figure  10.  This  gain  fall  off  approaches  the  cos  n  0 
function  where  3/2>n>l  which  indicates  a  well  matched  radiator. 

For  the  receive  case,  the  technique  using  an  external  remote 
source  in  the  far  field  for  measuring  amplitude  and  phase  is 
shown  in  Figure  14.  Similar  considerations  exist  for  transmit, 
but  with  the  turned-on  element  transmitting;  the  received  signal 
is  detected  at  the  source  location  filtered  and  brought  back  to  „ 

the  receiver,  as  shown  in  Figure  15.  To  carry  out  receiving  and 
transmitting  calibrations,  the  64  element  positions  were  * 

traversed  and  the  digital  phase  shifters  stepped  through  each  bit 
setting  with  phase  and  amplitude  (Aj_)  being  recorded. 

420 


-10- 


4.2  Phase  Calibration 


For  phase  calibration  the  objective  is  to  have  equal  insertion 
phase  from  the  input  of  the  source  to  the  output  of  the  feed  or 
vice  versa  when  each  element  position  is  interrogated  one  at  a 
time.  This  focuses  the  array  to  the  far  field.  The  geometry  is 
then  translated  to  the  cylindrical  surface  of  the  array  so  that 
the  beam  steexing  alogrithm,  equation  (Eq.  1)  stored  in  the 
discrete  module  controller  (DMC)  can  be  utilized.  The  geometric 
modification  factor  is  the  path  lengths  (Z^)  from  the  cylindrical 
surface  to  a  plane,  equation  (Eq.  2) . 

In  the  process  of  phase  calibration,  both  the  monitoring  of  the 
component  performance  and  the  isolation  and  correction  of  faults 
can  be  exercised.  For  example,  if  a  39  degree  phase  shift  is 
required  the  6,  11  and  22  degree  bits  would  be  set/  a  failed  11 
degree  bit  can  be  partially  compensated  for  by  choosing  a  45 
degree  command  reducing  the  error  from  11  to  6  degrees.  Here  a 
fault  has  been  located  and  partially  corrected. 

4 . 3  Amplitude  Calibration 

For  amplitude,  the  objective  would  be  to  achieve  the  specified 
illumination  functions  such  as  Taylor  for  sum  and  Bayliss  for 
difference.  On  AMIT,  where  only  a  one-step  attenuator  was 
implemented  in  the  module,  electronic  control  of  amplitude  to 
fine  increments  was  not  available;  instead,  fixed  attenuators 
were  inserted  in  series  with  the  module  to  obtain  the  desired 
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illumination.  The  amplitude  was  measured  and  the  results  used  to 
find  faults  in  the  element  radiator-to-output  port  transmission 
paths.  Corrections  were  made  by  attenuator  substitution  or 
insertion  and  T/R  module  replacement.  Where  electronic  BITE  or 
current  monitoring  indicated  that  the  modules  were  probably  good, 
the  fault  was  isolated  to  the  radiator.  The  end  result  was  phase 
and  amplitude  calibration  files  to  which  future  calibrations  were 
referenced  and  appropriately  adjusted.  Using  these  measured 
phase  and  amplitude  data  files,  antenna  patterns  were  calculated. 

4 . 4  Measured  Results 

An  example  of  such  a  calibration  for  the  receive  sum  port  for  a 
scanned  beam  was  taken  by  rotating  the  pedestal  by  40  degrees 
from  the  down  range  boresight.  This  resulted  in  a  plane  wave 
source  at  40  degrees  from  the  array  broadside,  as  shown  in  Figure 
16.  The  resultant  amplitude  illumination  as  compared  to  the 
ideal  illumination  for  an  n  =  10,  -50  dB  Taylor  is  shown  in 
Figure  17.  Correction  for  errors  was  limited  by  the  use  of  the 
fixed  attenuators.  Larger  errors  such  as  at  position  no.  42  may 
mean  a  low  gain  amplifier  or  other  component  deficiency. 

The  corresponding  far  field  pattern  as  computed  from  the  measured 
calibration  file  is  compared  to  the  actual  measured  far-field 
pattern,  as  shown  in  Figure  18. 

During  the  test  program,  receive  and  transmit  beams  were  formed 

and  steered  over  a  wide  scan  sector.  The  free  space  calibration 
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technique  was  used  to  derive  the  array  illuminations  and  produce 
collimated  sum,  difference,  and  transmit  beams.  The  performance 
characteristics  associated  with  the  beams  were  determined  as 
follows:  pattern  sidelobes,  aperture  gain,  amplifier  gain, 

effective  radiated  power  and  power  output.  The  measured  values 
of  these  parameters  were  in  agreement  with  performance 
expectations . 


5.0  CONCLUSION 

A  conformal  phased  array  containing  state-of-the-art  gallium 
arsenide  MMIC  T/R  modules  has  been  successfully  operated  in 
receive  and  transmit  both  in  the  laboratory  and  on  the  pattern 
test  range.  The  array  has  been  calibrated  for  CW  receive  and 
pulse  transmit  operation  using  free-space  far-field  range 
techniques.  It  is  planned  to  add  12-bi.t  digital  attenuators  to 
the  array  to  permit  more  accurate  amplitude  illumination 
calibration  and  to  provide  the  precision  needed  to  compensate  for 
amplitude  illumination  distortion  as  a  function  of  beamsteering 
[8]  . 
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Figure  4 .  Hardware  Inside  Conformal  Array  Enclosure 
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Figure  8 .  Linear  Polarized  Patch  Radiator  Configuration 
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Figure  12.  AMIT  Beam  Steering  Controller  Block  Diagram 
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Figure  14.  AMIT  Receive  Calibration  Setup  Using 
a  Remote  Source 
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Figure  15.  AMIT  Transmit  Calibration  Setup  Using 
a  Filtered  Signal  Link 
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Figure  17.  Measured  Illumination  for  Free  Space  Calibration 
40°  Scan  Angle  (Compared  to  ideal  n  =  10,  -50  dB 
Taylor) 
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ABSTRACT 

In  tlrs  paper,  a  new  frequency  locking  approach  for  active 
microstrip  antenna  arrays  is  introduced.  For  better  inter- injection 
locking,  a  two-dimensional  mutual  coupling  arrangement  of  the 
individual  elements  is  used.  This  spatial  arrangement  may  be  com¬ 
plemented  with  a  phase -lock  loop  to  provide  a  fine  tuning  of  the 
frequency- lock.  Experimental  results  obtained  with  Gunn  oscillators 
mounted  on  the  patch  antennas  are  presented.  The  measured  frequency 
spectrum  of  the  received  locked  signals  are  also  given. 

1.  INTRODUCTION 

Recent  developments  in  microwave  and  milrlimeter-wave  integrated 
circuits  and  solid-state  devices  make  it  possible  to  incorporate 
oscillators  on  a  microstrip  patch  antenna  [1-3].  Active  antenna 


arrays,  which  can  be  obtained  by  properly  combining  individual  active 
elements,  can  be  used  in  applications  where  large  power  solid-state 
microwave  and  millimeter-wave  sources  are  needed.  The  advantage  of: 
using  microstrip  technology  is  its  low  cost,  and  conformal  config¬ 
uration.  The  investigation  of  the  possibilities  of  combining  the 
power  from  such  large  arrays  is  an  important  task.  Since  each  active 
array  element  contains  an  independent  oscillator,  the  control  of  the 
coherence  of  the  large  power  source  becomes  a  major  issue.  Recently, 
several  ways  of  frequency  locking  of  individual  oscillators  were 
reported.  Both  spatial  and  mutual  injection  locking  have  been 
investigated  [1-3].  However,  the  injection  by  mutual  coupling  has 
the  advantage  over  injecting  a  signal  from  space  in  that  it  does  not 
require  an  additional  locking  signal. 

Hummer  and  Chang  investigated  both  active  patch  antennas  and 
two-element  arrays  [1].  Results,  obtained  at  10  GHz,  using  10  mW 
Gunn  diodes  as  active  elements,  have  shown  combining  efficiencies 
of  80  percent.  The  two- element  array  had  a  much  narrower"  pattern 
than  a  single  patch  pattern.  For  frequency  locking  of  the  active 
array  elements  results  were  reported  using  both  spatial  injection 
and  mutual  coupling  as  the  machanisms .  The  spatial  injection  of  a 
signal  from  space  proved  to  produce  a  narrower  bandwidth  than  the 
one  obtained  by  mutual  coupling  injection  locking.  Later  results 
[2],  reported  a  combining  efficiency  of  90  percent  for  both 
aperture -coupled  patch  antenna  arrays  as  well  as  directly  connecting 
the  active  elements  to  the  patches.  For  a  certain  bias  on  the  Gunn 
diodes,  the  pattern  exhibited  a  break  of  the  boresight  beam  into  two 
robes.  This  break  was  explained  as  a  loss  of  phase  lock. 
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Camillcri  and  Bayraktoroglu  [3]  presented  a  monolithic  active 
antenna  using  a  millimeter-wave  IMPATT  oscillator  operating  at  A3. 3 
GHz  with  a  CW  output  power  of  27  mW  and  a  7.2  percent  conversion 
efficiency.  Linear  array  radiating  elements  were  combined  and 
radiation  patterns  were  measured  as  functions  of  the  number  of 
elements  and  the  spacing  between  elements.  An  efficient  power 
combining  from  independent  sources  eliminated  the  need  for  lossy 
power  distribution  networks.  Inexpensive  millimeter-wave  power 
sources  can  be  built  by  incorporating  into  the  array  a  waveguide 
feed.  Because  the  radiation  pattern  of  a  single  element  transmitter 
is  similar  to  that  of  an  open-ended  waveguide,  there  was  an  efficient 
coupling  of  the  source  to  the  waveguide.  The  reported  efficiencies 
(of  the  order  of  10  percent,  which  is  the  efficiency  of  the  same 
IMPATT  structure  without  the  radiating  element)  indicate  that  the 
structure  worked  well  both  as  an  oscillator  and  a  radiator.  The 
arrays  of  oscillators  can  be  injection  locked  by  two  mechanisms; 
either  injection  locking  by  radiation  coupling,  or  injection  locking 
by  means  of  coupling  via  transmission  lines  between  elements .  Here, 
the  results  reported  were  obtained  using  the  second  mechanism.  An 
advantage  of  the  transmission  line  coupling  method  is  that  the 
coupling  strength  could  be  adjusted  by  adjusting  the  characteristic 
impedance  of  the  lines.  The  results  showed  that  for  two-  and 
three -element  arrays  the  injection  locking  mechanism  worked  very 
well  and  the  radiation  patterns  sharpened. 

Mutual  coupling  is  one  of  the  mechanisms  of  injection  locking 
of  adjacent  elements  in  an  active  array.  For  the  modeling  of  these 
effects,  Benalla  and  Gupta  [4]  investigated  a  multiport  network 
approach.  The  model  was  based  on  an  equivalent  mutual-coupling 
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admittance  matrix  network.  To  model  the  edge  fields,  the  elements 
of  these  admittance  matrices  were  found  by  using  equivalent  magnetic 
current  line  sources.  Using  this  model,  mutual  coupling  effects 
between  radiating  edges  of  the  same  patch  antenna  as  well  as 
neighboring  element  effects  may  be  incorporated.  The  model  results 
showed  good  agreement  with  the  measurement  data.  For  spacing  of  a 
quarter  of  free-space  wavelength  between  patches,  the  coupling 
coefficient  reached  values  of  20  dB.  One  of  the  important  conclusions 
was  that  the  coupling  between  nonradiating  edges  is  much  smaller 
than  that  from  radiating  edges.  This  effect  is  due  to  the  cosine 
voltage  distribution  along  the  nonradiating  edges,  thus  making  it 
possible  to  partially  cancel  the  positive  portion  by  the  negative 
portion  of  this  voltage  distribution. 

This  paper  introduces  a  novel  approach  of  frequency  locking  for 
active  microstrip  antenna  array.  It  consists  of  two  mechanisms.  The 
dominant  locking  method  is  the  interinjection  locking  by  mutual 
coupling. 


m  □ 
□  □  □ 
□  □ 


Figure  1.  Proposed  two- dimensional  arrangement 
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The  two-dimensional  arrangement,  shown  in  Fig.  1,  provides  for 
the  spatial  interaction  between  elements  in  a  two-dimensional  fashion. 
This  mechanism  is  complemented  with  a  phase -lock  loop  to  provide 
additi  nal  frequency  lock  fine  tuning. 

2.  TECHNICAL  BACKGOUND 

2 . 1  RECTANGULAR  MICROSTRIP  ANTENNAS 

The  rectangular  microstrip  antennas  have  been  extensively  studied 
and  analyzed  [5-6].  The  various  investigated  models  have  proven 
useful  for  different  applications  and  frequency  ranges. 

The  rectangular  patch  antenna  can  be  viewed  as  two  radiating 
slots,  whose  admittances  are  then  used  to  determine  the  radiation 
impedance  and  other  parameters.  The  distance  between  the  radiating 
slots  is  taken  to  be  half  of  the  guide  wavelength.  The  physical 
distance  L  is  adjusted  to  account  for  the  fringing  fields  that  occur 
at  the  radiating  edges.  These  fields  tend  to  increase  the  effective 
distance.  It  is  Lf//  that  is  to  be  half  a  wavelength.  Expressions  for 
the  dimensions  of  the  rectangular  patch  are  given  below: 


(D 


-2A  l 


(2-) 
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where  W  is  the  length  of  the  slot  and  L  is  the  spacing  between 
slots .The  correction  factor  A/  accounts  for  the  fringing  effects  at 
the  two  radiating  edges  of  the  antenna,  and  is  given  as: 

A  /  _ (c-.+  0.3)(l///i*0.264) 

It  m  °‘‘  U(e,-0.258)(l///i  +  0.8) 

2.2  INPUT  IMPEDANCE 

To  calculate  the  patch  input  impedance,  a  patch  model  must  be 
established.  One  of  the  models  which  proved  to  be  quite  accurate  and 
it  is  also  easy  to  implement  mathematically  is  the  cavity  modal - 
expansion  model,  where  the  patch  is  viewed  as  a  thin  cavity  with 
magnetic  walls  [6-8].  A  series  of  cavity  resonant  modes  is  then  used 
to  express  the  field  between  the  patch  and  the  ground  plane.  For  the 
rectangular  patch  shown  in  Fig.  2, 


Figure  2.  Rectangular  patch  antenna 


the  z- directed  electric  field  is  expanded  as  a  series  of  TM^  modes: 


*'«(*•  y)-  Z  Z/,»««Bm»(-v*y)  ('0 

m  n 

These  mode  functions  eM(t(\\ y)  satisfy  a  two-dimensional  wave-equation 
while  the  corresponding  eigenvalues  satisfy  the  equation: 


kn*  kl 


For  a  nonradiating  cavity 


(^) 


,  /nl  .  .  win 

knml “•  «/id  *M-— .  (6) 

For  the  radiating  cavity  A:n  and  A:m  are  slightly  different  than  the 

above  values  and  hence  kmn  becomes  complex.  The  complex  resonant 
frequency  becomes  u>mn  =  u>r  +  ;ca>,. 

The  input  impedance  becomes: 


y  tn  ,  ,,  .  .  V  V*  (•vo»  y0)^ 
=  =  kh  ^  L—rr—j—G 

•  0  m-On^O  ft- 


where  /i  is  the  height  of  the  dielectric,  k”ix>/c.  Also: 
>iii(/u(ctx/2li/)  sin(mjtdy/21) 


Gmn  =  ' 


/ikcZ,v/2I/ 


mndy/2L 


(7) 


and 


'flM 
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cos(  fcn.v/)  cos(  fcm  y ) 


where  x»tn  is  a  constant. 
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In  Fig.  3,  an  equivalent  electrical  network  is  shown. 


R 


lO 


Figure  3.  Equivalent  circuit  of  patch  antenna 

If  the  operation  of  the  patch  antenna  is  confined  to  a  small 
region  around  the  TM10,  the  network  can  be  simplified  and  the  input 
impedance  becomes: 


/(co/Cio) 

(A)2  -  (U)r  +  /(A),)2 


(8) 


where  C,0  «  ().SC4frcoo(iiy0//.)  and  Cttc  »  The  term  y  A'  <  represents 

a  net  series  inductance  that  accounts  for  the  contributions  from  the 
other  modes. 

The  complex  eigenvalue  A: l0  can  be  found  through  an  iterative 
process  which  solves  the  equation: 


in  n(A •,(,/.) 


?.k  loot  jo 
k  to  "  «  io 


(9) 
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where  : 


(10) 


o  h 

*  a,°” /T7”i7>  “ 

and  Yw  represents  the  radiating  wall  admittance  and  is 

Y  v *  (h/377)(I//X.0)+  ;O.OI668(A///i)(I‘//^oK  (H) 

To  account  for  different  W/L  ratios,  the  wall  admittance  may  be 
modified  by  a  factor: 

/•  y ( 1/  .  I. )  -  0 .7747  +  0 .5977 (  W t  L  -  1  )  -  0 . 1 638(  1/ / /.  -  1  )?  (12) 

The  solution  for  ki0  is  given  by: 

//.)~(A,//.)  (13) 

where 


2(a,0/.)(Jt- A.,)  A?, 

A  ,  - - -  ■  ■  — - - (\  44 

(«l0/.)^  2A;,u-A?,-uz  3  '  J 

For  A0,  a  starting  value  of  zero  is  recommended. 

Another  microstrip  patch  model  is  the  transmission  line  model 
which  assumes  that  the  two  radiating  edges  are  separated  by  a  low 
impedance  transmission  line  of  length  L.  The  expression  of  the  input 
impedance  of  a  microstrip  patch  with  arbitrary  feed  point  is: 
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?  V  r^fC0S(P/.,)-«-y^„8ill(P/.,)i  ^fCOS(lW.z)*^w»8ill(P/.2)"1 
^wcos((i/.|)+yJft8in(p/M)  i'wcos((W.2)+/^f8i»(|i/.2)_ 

+  /.V,  (IK) 

The  two  input  impedance  expressions  ,  Eqs .  (7)  and  (15),  are  in 
good  agreement. 


2.3  MUTUAL  CONDUCTANCE 


The  mutual  conductance  between  two  patch  antennas  is  used  to 
determine  the  array  directivity  and  its  input  power.  The  expression 
of  the  mutual  conductance  is  derived  from  the  total  radiated  power. 
Since  microstrip  planar  arrays  radiate  in  half  space,  the  mutual 
conductance  between  slots  of  length  W  placed  in  the  X-Y  plane  and 
aligned  in  the  Z  direction  is  given  by  [9]: 


(16) 


where  J0  is  the  Bessel  function  of  first  kind  and  zero  order.  Eq. 

(16)  can  be  extended  to  patches  with  arbitrary  position  in  the  X-Y 
plane  by  considering  contributions  from  both  edges  of  each  patch. 
A  single  expression  for  an  odd  mode  excitation  of  the  patch  antennas 
can  be  shown  to  be : 


W 


* 
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The  terms  in  Che  sharp  parenthesis  represent  mutual  conductances 
from  two  slots  spaced  at  distances  of  X,  X-L,  and  X-lL  apart.  Plots 
for  the  mutual  conductance,  obtained  by  numerical  integration  of  Eq. 
(17),  for  various  positioning  of  the  two  patches  are  shown  in  Fig 
4  (a) -(c). 


Figure  4(a) .  Mutual  conductance  of  two  rectangular  patches  when 
either  X  or  Y  are  set  to  zero. 


Figure  4(b) .  Mutual  conductance  of  two  ractangular  patches  when 


r  =  0.598\0. 
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Figure  4(c),  Mutual  conductance  of  two  rectangular  patches  when 

,Y«-0.511S\0. 

The  results  in  Fig.  4  (a)  are  similar, to  the  ones  shown  in  Ref.  [9], 
As  the  spacing  between  the  elements  increases  (see  Fig. 4  (b)-(c)), 
the  mutual  conductance  exhibits  a  slower  decrease  in  magnitude  with 
the  relative  position. 


2.4  TWO -PORT  PARAMETERS 


The  use  of  the  moment  method  enables  an  elegant  calculation  of 
the  two-port  parameters  of  two  microstrip  patch  antennas  [10).  The 
expressions  of  the  two-port  impedance  matrix  are: 
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(19) 


=  z,2=-L'X2> 

n*  I 

where  V is  the  voltage  due  to  source  i,  i  «  1,2,  and  /„  represents 
the  expansion  mode  currents. 

3.  EXPERIMENTAL  RESULTS 

Microstrip  patches  and  arrays  have  been  fabricated  for  the  X-band 
frequencies.  The  Gunn  diodes  used  operated  around  8  GHz.  The  placement 
for  the  position  of  the  Gunn  diodes  was  such  that  oscillation  condition 
occurred  [2]. 

In  Fig.  1,  the  geometrical  positioning  of  the  array  elements  is 
shown.  The  spacing  between  elements  was  a  half  of  guide -wavelength 
in  both  directions.  The  performance  of  wider  patche  arrays,  yielding 
better  coupling  between  elements  ,  was  also  investigated.  In  particular 
a  W  of  0.75  times  the  guide -wavelength  was  chosen. 

In  Fig.  5  (a) -(b) ,  the  measured  patterns  for  a  four  element  array 
excited  from  a  single  common  source  is  shown.  The  half  power  beamwidths 
were  less  than  50  degrees  fL r  the  11-plane  pattern,  and  less  than  35 
degrees  for  the  E-plane  pattern.  The  spatial  arrangement*  shows  no 
main  lobe  break  in  the  patterns.  These  results  show  that  this 
arrangement  is  suitable  for  active  arrays. 


Figure  5(a).  Relative  power  H-plane  pattern  for  a  four-element 

array . 
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Figure  5(b).  Relative  power  E-plane  pattern  for  a  four-element 

array. 


In  Fig.  6  (a)  and  (b) ,  the  measured  coupling  coefficient  between 
patches  of  different  positions  is  shown.  The  results  show  no  sig¬ 
nificant  difference  in  mutual  coupling  between  elements  regardless 
of  their  position  with  respect  to  the  radiating  edges.  Uniform 
two-dimensional  mutual  coupling  is  made  possible  this  way. 
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Figure  6(a).  Coupling  coefficient  for  two  patches  arranged  as 

shown. 
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Figure  6(b).  Coupling  coefficient  for  two  patches  with  adjacent 

non- radiating  edges. 

In  Fig.  7,  measured  input  impedance  an  array  element  where  the 
feed-point  is  approximately  .32L  from  one  of  the  radiating  edges  is 
shown.  The  position  of  the  feed  point  makes  it  possible  that  the 
resonant  condition  for  the  Gunn  diode  is  achieved. 
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Figure  7 .  Measured  input  impedance  of  an  array  element  in  a  seven 

element  array. 
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The  results  obtained  for  three  active  microstrip  patches 


positioned,  as  shown  in  Fig.  8,  are  given  in  Fig.  9  -  12. 


Figure  8,  Two-dimensional  arrangement  of  a  three -element  active 


array.  KKn  ?.ddh  gh* 
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Figure  9.  Frequency  spectrum  of  element  #1  in  the  array. 


In  Fig.  9,  the  frequency  spectrum,  obtained  when  only  the  Gunn 
diode  of  element  ttl  was  biased,  is  shown.  The  radiation  frequency 
is  7.996  GHz. 
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Figure  10.  Frequency  spectrum  of  element  # 2  in  the  array. 

In  Fig.  10,  the  frequency  spectrum,  obtained  when  the  Gunn  diode 
of  element  #2  was  biased,  is  shown.  The  radiation  frequency  is  8.006 
Gllz.  MKH  7. five  OMz 
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Figure  11.  Frequency  spectrum  of  elements  #3  in  the  array. 


In  Fig.  II,  the  frequency  spectrum  for  the  biased  Gunn  diode  of 
element  3  is  shown.  The  radiation  frequency  is  7.972  GHz.  The 
differences  in  radiation  frequencies  are  less  than  34  MHz. 

When  all  the  Gunn  diodes  were  biased,  the  three  oscillators 
injection  locked  to  the  frequency  of  7.948  GHz,* as  shown  in  Fig. 

12.  MKR  7.Q-IU  Gllz 


Figure  12.  Frequency  spectrum  of  injection- locked  three-element 

array. 


The  calculated  combining  efficiency  is  around  87.5  percent.  These 
results  show  that  the  new  arrangement  provides  for  efficient  injection 
locking  by  mutual  coupling.  Additional  measurements  involving  larger 
arrays  are  in  progress. 


453 


4.  SUMMARY  AND  CONCLUSIONS 


In  this  paper,  a  new  spatial  arrangement  for  a  two-dimensional 
microstrip  active  antenna  array  is  introduced.  The  measurements  of 
the  array  elements  input  impedance  and  coupling  coefficient  show 
that  this  arrangement  is  suitable  for  two-dimensional  mutual-coupling 
injection-locked  active  array.  This  approach  is  superior  to  others 
for  its  two-dimensional  array  capabilities.  A  three-element  active 
array  was  fabricated  and  tested.  The  active  elements  were  Gunn  diodes 
with  oscillation  frequencies  around  8  GHz.  The  three  active  elements 
locked  efficiently  via  mutual-coupling.  The  introduction  in  the  array 
of  an  element  provided  with  phase-lock  loop  will  give  a  high  precision 
frequency  locked  array.  This  additional  element  has  to  be  coupled 
to  only  one  of  the  elements  of  the  array  in  order  to  injection  lock 
all  array  elements.  Therefore,  the  mutual-coupling  injection-lock 
mechanism  together  with  a  phase-lock  loop  mechanism  provide  for  a 
high  precision  frequency  locking  method.  Results  to  support  this 
statement  are  in  progress. 
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ABSTRACT 


1 


The  constant  demand  for  more  sophisticated  antennas  requires 
the  use  of  powerful  simulation  facilities  . 

By  analogy  with  the  circuit  field,  we  have  developped  PHEDRE,  a 
CAD-CAM  software  tool  to  design,  characterize  and  optimize  the 
behavior  of  single  antennas,  arrays  and  phased  arrays  with  their 
radome  and  mounted  on  structures  (eg.  airborne  antennas). 

This  paper  describes  the  structure  and  possibilities  of  PHEDRE. 

The  close  interconnection  of  4  modules  allows  the  user  to  choose  the 
appropriate  path  to  conceive  and  describe  the  antenna  . 

These  modules  have  been  devised  with  the  idea  of  being  adaptable 
and  easy  to  use. 

The  computer  software  has  been  enlarged  to  facilitate  exchanges  : 

-  with  other  antenna  computer  programs  (eg.  wire  antennas  through 
the  use  of  EFIE  to  define  new  radiators,  antennas  with  plates,..) 

-  with  automatic  measurement  and  automatic  manufacturing  tools. 

The  oral  presentation  will  be  illustrated  with  the  designs  of: 

-  a  very  low  sidelobe  slot  array  antenna  and  its  radome 

-  a  phased-array  antenna  including  the  study  of  active  elements  . 
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INTRODUCTION 


2 


The  study  of  the  association  of  radiating  elements  with  a 
radome  requires  the  use  of  powerful  simulation  facilities. 

By  analogy  with  microwave  circuits  design  softwares,  we  have 
developped  a  Computer-Aided  Design  software  tool  for  antennas, 
taking  into  account  the  following  features  : 

-  modularity  :  to  allow  the  adaptation  of  the  softwares  to  the 

various  antenna  problems, 

-  possibility  of  updating  in  relation  to  future  problems, 

-  possibility  of  linkage  to  other  facilities  through  the  exchange 
of  data  files  . 

The  basic  concept  of  this  software  lies  in  the  computation 
of  the  response  of  a  set  of  radiating  elements  to  an  arbitrarily 
polarized  incident  wave  . 

On  the  side  of  this  main  axis,  specialized  modules  have  been 
designed  to  master  the  influence  of  external  and  internal  couplings, 
of  the  environment  and  of  the  radomes. 

Figure  1  shows  the  general  structure  of  ti.e  software  and  the 
associated  "peripherals"  described  in  the  following  sections. 


2.  CONNECTION  TO  EXTERNAL  FACILITIES. 


3 


Connections  to  external  facilities  are  : 

-  The  exchange  of  datas  with  test  facilities,  in  order  : 

.  to  facilitate  the  comparison  of  results, 

.  to  use  the  power  of  the  central  computer  when  the  local  ones 
are  not  powerful  enough, 

.  to  easily  obtain  program  inputs  from  measured  datas  . 

The  test  facilities  can  be  : 

.  circuit  test  facilities,  to  characterize  the  feeding  system, 

.  antenna  test  facilities,  to  get  far-field  or  near-field 
(  planar  or  cylindrical  )  measured  data  available. 

-  The  use  of  computer  results  from  circuit  simulation  programs 
(  eg.  SUPER  COMPACT)  to  anali^e  the  feeding  of  antennas 

-  The  link  with  numerically  controlled  manufacturing  machines,  in 
order  to  facilitate  the  exchange  of  the  mechanical  dimensions 

(  eg.  tne  link  to  an  automatic  milling  machine  for  a  slot  array 
antenna  )  and  to  collect  the  datas  from  the  mechanical  checking 
machines  after  the  manufacturing  of  the  antenna,  as  a  feed-back 
to  the  programs  . 
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3.  MAIN  FRAME  OF  THE  PROGRAM. 


4 


The  basic  idea  lies  in  the  expression  of  the  far  field  of  an 
array  as  a  function  of  the  sum  of  fields  originated  from  the 
various  elementary  radiators,  identical  or  different  in  nature, 
placed  in  whatever  manner  in  space,  to  which  different  complex 
weights  are  assigned,  when  the  whole  geometry  is  known. 

In  order  to  do  this,  two  libraries  have  been  established: 

-  a  library  of  radiators  defined  by  : 

.  analitical  formulations, 

.  files  created  by  the  use  of  the  Method  of  Moments  to  describe 
wire  antennas  or  antennas  made  of  plates 
.  files  issued  from  measurements  . 

-  a  library  of  illuminating  functions  with  : 

."standard"  laws  suitable  for  creating  sum,  difference  patterns 
for  rectangular  or  elliptical  antennas,  cosecant-squared,..) 
."advanced"  laws,  by  the  use  of  specific  algorithms 
(  eg.  adaptive  antennas  ). 

From  the  previous  datas,  a  certain  amount  of  parameters  can  be 
be  obtained,  for  example  : 

-  the  response  of  each  radiator  to  an  arbitrarily  polarized 
incident  wave , 

-  the  global  field  radiated  by  the  array  (  figure  2  ). 


A  specific  module  completes  the  previous  results  in  order  to  5 
give,  besides  the  radiation  pattern,  informations  about  : 

-  the  gain, 

-  the  effects  of  errors  on  the  amplitudes  and  the  phases,  by 
successive  random  samplings  or  by  the  use  of  a  statistical  model 

-  the  ecartometry  parameters  for  monopulse  antennas. 


4.  COMPLEMENTS  FOR  THE  PARASITIC  PHENOMENA 


The  good  working  condition  of  an  antenna  requires  the  control 
of  the  following  elements  : 


4.1.  COUPLINGS. 


-  internal  couplings  associated  with  the  feeding  elements  of 
the  antenna.  They  are  easily  linked  to  the  external 
couplings  in  a  matrix  form,  in  order  to  proceed  to  the 
analysis  of  the  antenna; 

-  external  couplings  between  the  radiators,  resulting  in  the 
modification  of  their  weights  or  their  active  patterns. 

In  some  cases,  as  for  the  slot  array  antennas,  a  process  of 
synthesis  can  be  carried  out  in  order  to  determine  the  geometry 
of  the  antenna,  to  manufacture  a  mock-up  and  to  execute  an 
operation  of  Analysis  with  the  real  mechanical  dimensions, 
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in  order  to  validate  the  analitycal  design.  6 

Aii  example  of  comparison  between  simulation  and  measurement 
for  an  antenna  composed  of  five  waveguides,  each  with  30  slots 
fed  in  progressive  mode  (  figure  3  )  ,  with  first  sidelobes  at 
15  dB  below  the  isotropic  level,  is  shown  on  figure  4. 

4.2.  ENVIRONMENT. 


The  influence  of  a  structure  on  the  radiation  pattern  of  an 
antenna  is  analyzed  by  the  use  of  models  based  on  the  Uniform 
Theory  of  Diffraction  (UTD)  in  the  following  steps  (figure  5); 

-  definition  of  the  exact  geometry  of  the  problem  (  through  the 
use  of  appropriate  software  tools,  eg.  CATIA  developped  by 
DASSAULT  SYSTEMES)  or  a  geometry  approximated  by  canonical 
surfaces  ; 

-  definition  of  the  different  rays  to  take  into  account  between 

a  source  and  an  observation  point,  according  to  the  phenomenon 
(  reflection,  diffraction,  multiple  interaction  ) 

-  computation  of  the  transfer  functions  along  the  different  paths 
with  the  use  of  a  library  of  electromagnetic  models  ; 

-  re-composition  of  the  total  field  from  the  various  components 

4.3.  RADOME. 

The  presence  of  a  radome  in  front  of  an  antenna  can  be  studied 


by  specific  programs  using  : 


-  Geometrical  Optics  formulation  for  large  antennas  7 

-  a  near-field  technique  for  medium-sized  antennas 

-  the  use  of  the  reciprocity  theorem  for  small  or  medium-sized 
antennas,  to  take  into  account  the  reaction  antenna- radome . 

The  figure  6  (right)  shows  an  example  of  the  influence  of  a 
radome  optimized  for  the  antenna  shown  above. 

5.  CONCLUDING  REMARKS 

This  paper  has  shown  the  possibility  to  extend  the  principles 
of  Computer  Aided  Design  in  order  \o  obtain  a  general  treatment 
method  for  current  and  future  array  antennas. 
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Abstract 

The  input  impedance  of  a  finite  slot  in  the  groundplane  of  a  stripline 
is  numerically  studied  using  the  method  of  moments.  The  problem  is 
cast  so  that  the  unknowns  to  be  determined  are  the  electric  currents  on 
the  stripline  center  conductor  and  the  equivalent  magnetic  currents  in 
the  slot.  The  cases  where  the  slot  is  present  in  only  one  (the  “unilateral” 
case)  or  both  (the  “bilateral”  case)  of  the  groundplanes  are  each  included 
in  the  analysis. 

The  basis  functions  used  in  this  study  are  the  triangle  subdomain 
functions  introduced  by  Rao,  Wilton  and  Glisson.  The  advantage  of 
this  choice  is  that  the  slot  and  strip  geometries  which  can  be  analyzed 
include  any  shape  whose  boundary  can  be  adequately  approximated  by  a 
polygon.  This  class  includes  practically  all  shapes  of  engineering  interest. 

The  accuracy  of  the  method  is  verified  by  comparison  to  measured 
and  theoretical  results.  The  application  of  this  analysis  technique  to  the 
design  of  broad  band  stripline/slotline  transitions  is  also  discussed. 


1  Introduction 

Tapered  notch  antennas,  otherwise  known  as  Vivaldi  or  LTSA  (Linearly  Ta¬ 
pered  Slotline  Antennas),  have  been  described  by  a  number  of  workers[l]-[5]. 
The  recent  popularity  of  these  elements  for  use  in  phased  arrays  is  due  to  their 
wide  impedance  bandwidth,  ease  of  manufacture,  low  cost,  light  weight,  and 
small  size.  Current  designs  are  mostly  empirical  due  to  the  lack  of  an  adequate 
theoretical  model  for  notch  elements.  Due  to  the  complexity  of  the  elements, 
most  analyses  to  date[6,7,8]  have  had  to  resort  to  rather  crude  approxima¬ 
tions,  which  are  not  sufficiently  accurate  for  many  design  purposes.  Raytheon 
is  engaged  in  an  ambitious  effort  to  rigorously  analyze  an  array  of  single-  or 
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Notch 


Figure  1:  The  geometry  treated  in  this  note  consists  of  a  stripline  feeding  a 
notch  in  one  or  both  of  its  ground  planes. 

dual-polarized  notch  elements  excited  by  stripline  feeds.  The  analysis  will  be 
able  to  handle  arbitrarily  shaped  stripline  and  notch  regions. 

This  paper  reports  on  some  preliminary  work  which  is  the  first  step  in  the 
direction  of  the  full  notch  array  analysis.  The  geometry  treated  here  consists 
of  a  stripline  of  arbitrary  shape  exciting  an  arbitrarily  shaped,  finite  aperture 
in  one  or  both  of  its  groundplanes.  The  motivation  for  this  work  was  to 
develop  and  demonstrate  the  capability  to 

1.  model  electric  and  magnetic  currents  using  triangle  subdomain  basis 
functions  on  arbitrarily  shaped  domains. 

2.  treat  the  coupling  between  strip  and  aperture  currents. 

3.  accurately  compute  the  input  impedance  in  the  stripline  feed. 

These  goals  have  been  attained.  In  addition,  we  describe  in  Section  6  how 
a  slight  modification  to  the  code  results  in  a  useful  design  tool  for  the  strip¬ 
line/slotline  transition  portion  of  the  notch  element. 

The  structure  to  be  analyzed  is  shown  in  Figure  1.  It  consists  of  a  balanced 
stripiine  with  a  finite  aperture  (or  slot)  cut  in  one  (the  “unilateral”  case)  or 
both  (the  “bilateral”  case)  of  its  ground  planes.  Both  the  stripline  center 
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conductor  and  the  slot  may  be  of  irregular  polygonal  shape,  although  the 
striplinc  assumes  a  uniform  width  at  some  distance  from  the  slot.  If  slots 
are  present  in  both  ground  planes  then  they  must  be  identical  in  shape  and 
location.  Notice  also  that  the  stripline  extends  only  a  finite  distance  beyond 
the  slot  in  the  +2  direction. 


A 


2  Source  Modeling  and  Input  Impedance 


An  e>wl  time  convention  is  assumed  and  suppressed.  RMS  phasors  are  used  for 
all  fields  and  sources.  The  quantity  of  interest  is  the  input  reflection  coefficient 
at  s  =  0  in  the  stripline,  which  is  assumed  to  be  far  enough  into  the  uniform 
width  region  and  far  enough  from  the  slot  so  that  the  fields  here  in  the  stripline 
can  be  written  as1 


E(®,  i/,  z)  =  {ae~]kz  +  e(®,  y) 

H(®,  y,  z)  =  {ae~:k:  -  be}kz^j  h (x,y) 


(1) 

(2) 


where  e  is  the  transverse  electric  field  of  the  TEM  stripline  mode,  h  =  zXe/77, 
and  ?/  is  the  intrinsic  impedance  of  the  dielectric,  a  and  b  are  the  usual  travel¬ 
ling  wave  variables  whose  ratio  gives  the  input  reflection  coefficient,  Tin  =  b/a. 
To  determine  T;,,  we  first  place  a  Perfect  Magnetic  Conducting  (PMC)  wall 
at  z  =  0  and  impress  on  it  at  z  =  0+  an  electric  surface  current 


J'  =  zxh. 


(3) 


This  current  impressed  on  the  PMC  wall  is  the  rigorous  field  theory  realization 
of  an  open  circuit  current  source  which  will  ensure  that  the  magnetic  field  at 
2  =  0+  is  exactly  h.  One  could  then,  in  principle,  find  the  input  impedance 
by  computing  the  electric  field  in  the  source  plane.  However,  knowledge  of  the 
stripline  current  I(z)  is  also  sufficient  to  determine  the  reflection  coefficient 
since,  in  the  vicinity  of  2  =  0 


(4) 


Once  the  current  on  the  strip  is  computed  using  the  moment  method  procedure 
described  below,  equation  (4)  is  sampled  at  N  >  3  locations  in  2.  This  results 
in  an  overdetermined  set  of  linear  equations  in  the  unknowns  ij  and  T-ml£. 
A  least  squares  solution  then  gives  the  input  reflection  coefficient. 


'In  equations  (1)  and  (2)  we  have  ignored  the  contribution  of  a  propagating  TEM  parallel 
plate  mode  which  would  be  excited  in  the  unilateral  case.  Due  to  opposing  symmetries  this 
mode  cannot  couple  with  the  striplinc  mode  and  thus  its  inclusion  would  have  no  effect  on 
either  the  induced  stripline  currents  or  the  input  impedance,  which  is  defined  with  respect 
to  the  stripline  mode. 


3 


479 


3  Equivalent  Currents 

First  the  center  conductor  (assumed  to  be  infinitely  thin)  is  removed  and 
replaced  with  an  equivalent  electric  surface  current,  Js. 

If  only  one  ground  plane  contains  a  slot,  we  assume  that  it  is  the  one 
located  at  y  =  h .  If  both  groundplanes  contain  slots  (the  bilateral  case)  then 
by  symmetry  we  restrict  consideration  to  the  region  y  >  0  by  enforcing  a  PMC 
boundary  condition  at  y  =  0.  The  induced  electric  surface  current  flows  just 
above  this  PMC  wall  at  y  =  0+.  Therefore  in  both  the  unilateral  and  bilateral 
cases  we  deal  with  only  a  single  slot. 

As  was  done  with  the  center  conductor,  the  slot  is  removed  (filled  in  with 
'"'EC  or  Perfect  Electric  Conductor)  and  the  necessary  magnetic  surface  cur¬ 
rents  are  defined  so  that  the  total  fields  are  unchanged.  For  fields  in  the  region 
|t/|  <  h  (interior)  the  required  source  is  a  magnetic,  surface  current  flowing  at 
y  =  Ir  given  by 

My‘W)  =  yxE(*,M).  (5) 

For  fields  in  the  region  y  >  h  the  equivalent  source  is  a  magnetic  surface 
current  flowing  in  the  y  =  h+  plane  given  by 

Mfl(s,*)  =  -Mf(x,2).  (6) 

The  analysis  now  reduces  to  determining  the  induced  electric  and  mag¬ 
netic  surface  currents,  Js  and  M"u ,  which  is  accomplished  using  the  method 
of  moments.  The  method  of  moments  formulation  involves  solving  integral 
equations  whose  kernels  are  Green’s  functions;  the  next  section  describes  the 
required  Green’s  functions  and  their  efficient  evaluation. 

4  Green’s  Functions 

Green’s  functions  for  three  separate  geometries  are  required.  The  first  geom¬ 
etry  is  the  PEC/PEC  parallel-plate  waveguide  structure  with  a  PMC  wall  at 
z  =  0.  Only  the  case  where  source  and  observation  point  are  in  the  region 
z  >  0  is  of  interest.  This  geometry  will  be  needed  to  analyze  the  unilateral 
slot  case.  The  second  geometry  is  the  PMC/PEC  parallel-plate  waveguide 
with  a  PMC  wall  at  z  =  0,  needed  for  the  case  of  a  bilateral  slot.  The  final 
geometry  is  the  half-space  y  >  li  with  a.  PEC  boundary  condition  at  y  =  h. 

4.1  PEC/PEC  Parallel-Plate  Green’s  Functions 

We  first  treat  the  case  of  an  infinite,  dielectric  loaded,  PEC/PEC  parallel- 
plate  waveguide  with  plate  separation  b  =  2h  which  lacks  the  PMC  wall  at 
z  -  0.  The  effect  of  the  PMC  wall  is  easily  accommodated  later  using  images. 
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The  material  between  the  plates  is  dielectric  of  permeability  no  and  per¬ 
mittivity  e  =  erco,  with  corresponding  wavenumber  k  =  oc  at  the  radian 
frequency  u.  The  electric  and  magnetic  fields  in  the  waveguide  satisfy 

E  =  -ju>A  -  V$  +  e-1VxF  (7) 

H  =  ju  F-V$+/io1VxA  (8) 

where  A  is  the  magnetic  vectoi  potential,  $  the  electric  scalar  potential,  F  the 
electric  vector  potential,  and  $  the  magnetic  scalar  potential.  The  potentials 
are  taken  to  satisfy  the  Lorentz  gauge,  so 


(V2  +  fc2)A  = 

-MoJ 

(9) 

(V2  +  fc2)$  = 

-a/e 

(10) 

(V2  +  &2)  F  = 

eM 

(11) 

(V2  +  &2)$  = 

-<rm/n  0 

(12) 

where  J  and  a  are  the  electric  current  and  charge  densities,  respectively,  and 
M  and  <rm  are  the  corresponding  magnetic  quantities.  The  currents  and 
charges  are  related  via  the  equations  of  continuity 

V  •  J  +  jua  =  0,  (13) 

V-M  +  jWm  =  0.  (14) 

The  boundary  conditions  on  the  potentials  are  stated  in  Table  1.  A  derivation 
can  be  found  in  [9]. 


PEC  Wall 

.PMC  Wall  | 

$  =  0 
af  _  n 

FK  ~  0 

nxA  =  0 

V„  •  A  =  0 
nx(VxF)  =  0 
n  >F  =  0 

$  =  0 
a<j»  _  n 

nxF  =  0 

V„  •  F  =  0 
nx(VxA)  =  0 
n  •  A  =  0 

Table  1:  Boundary  Conditions  on  Potentials 

The  potential  Green’s  functions  which  are  needed  are  those  for  Ax>  Az, 
Fx,  Fz,  and  'J'.  Note  from  Table  1  that  all  these  functions  satisfy  either 
Dirichiet  ( U  =  0)  or  Neumann  (f^  =  0)  boundary  conditions  at  the  paral¬ 
lel  plates.  Therefore,  we  define  Dirichiet  (Gj)  and  Neumann  (Gn)  Green’s 


functions  as  the  solutions  to 

(V2  +  &2)G(r;r')  =  -6(r  -  r')  (15) 

satisfying  the  radiation  condition  and  the  appropriate  boundary  condition 
on  the  parallel  plates.  In  (15)  r  =  ( x,y,z )  is  the  observation  point  and 
r'  =  (x',y',z')  is  the  source  point. 


4.1.1  Dirichlet  Boundary  Condition 


A  modal  series  expansion  for  G(i  is 


r')  =  ^  f;  sin  (g(«  +  /.))  sin  (f  (v'  +  *))  Ao(<W(*  -  W  +  (-'  -  *')2) 

(1G) 


ns  1 


where  Kq  is  the  modified  Bessel  function  of  order  0  and 


CVn  = 


(17) 


Equation  (16)  is  useful  when  the  source  and  observation  points  are  well  sep¬ 
arated.  An  alternative  series  can  be  obtained  by  imaging  the  source  through 
the  parallel-plate  walls  and  summing  the  contributions  of  the  images,  which 
may  now  be  considered  to  be  radiating  in  a  medium  without  boundaries.  This 
method  yields  a  second  series  expansion  of  Gd  to  be 


(18) 


where  P+  and  Rn  are  distances  between  observation  point  and  nth  order 


image  points  given  by 

Rn  =  \/(x  -  a')2  +  (z  -  2')2  +  (y-y'  -  4n/i)2  (19) 

R-  =  y/(x  -  s')2  +  {z  -  s')2  +  (y  +  y'  -  (2 n  +  1)2 K)\  (20) 

A  technique  due  to  Richards[10]  for  combining  (16)  and  (18)  to  arrive  at  a 
single,  numerically  efficient  formula  valid  for  all  r  and  r'  is  described  in  [11] 
and  used  here.  For  c  >  0,  define  pc,  P+,  and  P~  as 

Pc  =  x  -  a-')2  +  {z-  z'f  +  c2"  (21) 

pn  =  \jp\  +  {y-  y'  -  4n/02  (22) 

P~  =  \fpfHyTy'  -(2w  +  l)2/i)2  (23) 
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which  amounts  to  replacing  (x-x1)2  in  (19)  and  (20)  by  (a:— x')2  +  c2*  Clearly, 
making  this  substitution  in  the  scries  found  in  (16)  and  (18)  preserves  the 
equality  of  the  two  series.  An  additional  benefit  is  that  for  any  c  >  0  the 
series 


1  00 

G(Jn\rr,  O  =  xzr  £  sin  ($(y  +  h))  sin  (%(y'  +  h))  K0(anPc)  (24) 

1  1  71  =  1 

is  convergent  even  when  r  and  r'  coincide.  Application  of  a  Kummer  trans¬ 
formation  to  (18)  then  results  in  the  final  formula 

Gd(  r;r')  =  (/5m)(r;r')  + 

1  fe-JkR*  e-jkP,t  e-jkR~  e-jkP~  \ 

4?  n£o  ~~ +  flr  j  ’  (2o) 


for  the  unbounded,  parallel-plate,  Dirichlet  Green’s  function.  The  constant  c 
is  adjusted  so  that  the  series  in  (2d)  and  (25)  converge  rapidly. 

The  derivative  of  6',;  with  respect  to  y  is  also  required,  and  is  given  by 


06',, ,  agoo 

17(r;r)=-aT(r;r)  + 


= [•"  (#  *  b>) 

( y  +  y'  -  (2 n  +  l)2/i) 


+ 


-jkRt  j 


jk 


J - —  g-jkPn 

-2  T  3  C 


IfiT 


jfe  +  _i 

.Pn  2  P7 


(26) 


where 


1  00 

^2  S  ncos  (if  (y  +  *))  si»  (if  (s'  +  &))  p'o(ttupc). 

n=l 


(27) 


The  Dirichlet  Green’s  function  and  its  y  derivative  are  then  evaluated  for 
source  and  observation  points  in  the  planes  y  =  0  and  ij  =  h.  Note  that  of 
the  special  cases  of  interest  only  Gd(x,0,3;x',0,s')  and  ^^-{x,h,z\x' ,Q,z') 
are  nonzero. 
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4.1.2  Neumann  Boundary  Condition 

A  similar  procedure  can  be  followed  for  the  Neumann  Green’s  function.  The 
modal  expansion  is 

6»m)(r;  r')  =  —itfHkp o) 

+  :r-r  £  cos  (|f(J/  +  /0)  cos  +  /())  A'o(an/>o)  (28) 


where  //q2*  is  the  zero  order  Hankel  function  of  the  second  kind.  The  multiple 
image  expansion  is 


...  1  e-it-a- 


Combining  these  expressions  as  before  we  arrive  at  the  accelerated  formula 
for  the  Neumann  Green’s  function 


Gn(  r;r')  =  C7<m)(r;r0  + 

iff2 


+  fiT" 


where 


f(*%) 

+  ^  £cM§(j+,0)c°MSW  +  M)tfo(a»/>.)-  (31) 

‘  7(=l 

The  derivative  of  G'71  with  respect  to  y  is  given  by 

[,-« ( £  .  £)  - .-»«  (f  .  j, 

[•““  (l>*  5>) "  iw'w).. 
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where 


f)  fi(m)  1  o° 

-fir (r>r')  =  -^2  £nsin  W»  +  fe))cos(^(»'  +  h))K0(anpe). 

y  n=l 

(33) 

Again  we  evaluate  Gn  and  its  y  derivative  for  y  and  y'  only  at  0  or  h. 
Among  these  special  cases,  those  yielding  nonzero  values  are  Gn(®,  0,  z\  x', 0,  z')> 
®  » h, ^  ),  Gn(®, /i,  2;  ®  ,  0, 2  ),  G'n(®) h, 2; ®  , /i, 2  )  and  -g^4,(x, 0, 2;- 

x‘,h,z'). 

4.1.3  Potential  Green’s  Functions  with  PMC  wall  at  2  =  0 

We  define  the  location  of  the  image  of  the  source  point  in  the  PMC  wall  at 
2  =  0  as 

r"  =  (x',y',-z').  (34) 

Then  the  various  potential  Green’s  functions  for  a  PEC/PEC  parallel  plate 
waveguide  with  a  PMC  wall  at  2  =  0  are  given  in  terms  of  Gn  and  Gd  by  the 
following  expressions. 


G»(r; 

r') 

=  r1  {Gd(r 

r1)  +  Gd(r-,r")) 

(35) 

GA,(r\ 

*0 

=  Ho(Gd(r\ 

r')  +  Gd(r;r")) 

(36) 

G^„(r; 

r') 

=  Ho{Gn{r\ 

r')  +  Gn(r;r")) 

(37) 

*0 

=  Mo  (Gd(r; 

r')-Gd(ry)) 

(38) 

<3!f(r; 

*0 

=  Mo1  (^«(* 

•;r >')-Gn(r;r")) 

(39) 

G'g(r; 

r') 

=  -e(G„(r; 

r')-G„(r;r")) 

(40) 

V) 

=  -e  (Gd{  r; 

r')-Gd(r;r")) 

(41) 

V) 

=  -e  (Gn(i; 

r')  +  G„(r;r")). 

(42) 

The  superscript  “int”  used  above  distinguishes  these  interior  magnet  source 
Green’s  functions  from  the  exterior  region  Green’s  functions  described  in  Sec¬ 
tion  4.3.  The  subscript  describes  which  potential  quantity  the  Green’s  function 
represents. 

4.2  PEC/PMC  Parallel-Plate  Green’s  Functions 

We  now  treat  a  parallel-plate  waveguide  consisting  of  the  dielectric  as  before 
between  a  PMC  wall  at  y  =  0  and  a  PEC  wall  at  y  =  h.  One  sees  from 
examination  of  Table  1  that  all  the  potentials  satisfy  either 
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(nd)  Neumann  boundary  conditions  at  y  =  0  and  Dirichlet  boundary  condi¬ 
tions  at  y  =  ft 


or 

(dn)  Dirichlet  boundary  conditions  at  y  =  0  and  Neumann  boundary  condi¬ 
tions  at  y  =  ft. 


4.2.1  Neumann-Dirichlet  Boundary  Conditions 


The  Green’s  function  satisfying  Equation  (15)  with  Neumann  boundary  con¬ 
ditions  at  y  =  0  and  Dirichlet  boundary  conditions  at  y  =  ft  is  denoted  by 
Gmi(r;  r').  The  following  formulas  give  Gnd  and  its  y  derivative  (when  they 
are  not  identically  zero)  for  the  special  cases  of  interest. 


Gnd(x,0,z\x',0,z')  =  2Gd(x,0,z‘,x\0,z') 
-^(*,ft,z;*',0,*')  =  2-^(x,  ft,  z;  x'y  0,  z1). 


(43) 

(44) 


4.2.2  Dirichlet-Neumann  Boundary  Conditions 


The  Green’s  function  satisfying  Equation  (15)  with  Dirichlet  boundary  con¬ 
ditions  at  y  =  0  and  Neumann  boundary  conditions  at  y  =  ft  is  denoted  by 
Gdn(x\t').  Gdn  and  its  y  derivative  (in  the  only  special  cases  for  which  they 
are  not  identically  zero)  are 


Gdn(X)h,  ZJ  X  ,h,z) 


9Gdn 

dy 


(s,  0,2;  s',  ft,*') 


0,z;x',0,z') 

9Gnd  /i  l  n  i\ 

— 0^-(x,h,z;x  ,0,2 ). 


(45) 

(46) 


4.2.3  Potential  Green’s  Functions  with  PMC  wall  at  z  =  0 

The  various  potential  Green’s  functions  for  a  PMC/PEC  parallel  plate  wave¬ 
guide  with  a  PMC  wall  at  z  =  0  are  given  by  the  following  expressions. 


<?*(r;  r')  =  r1  (Gw(r;  r')  +  Gnj(r;  r"))  (47) 

C?/4,(r;r')  =  n0(Gnd(r,r')  +  G„<i(r;r"))  (48) 

Gvt^rjr')  =  noiGdniTir^  +  Gdnfar"))  (49) 

G^rjr')  =  fio(Gnd(r;r')-Gnd(r\r"))  (50) 

G'Pt(r;  r')  =  fio1  (Gdn(r;  r')  -  Gdn{r\  r"))  (51) 

Gg( r;r')  =  -e(Gdnfrt')  -  Gdn(r;r"))  (52) 

GfcfoiO  =  -€(Gnrf(r;r')-Gnd(r;r"))  (53) 

Gjf(r;rO  =  (Gdn(*\ r')  +  (7jn(r;  r")) .  (54) 
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As  before,  superscript  “int”  used  above  distinguishes  these  interior  magnetic 
source  Green’s  functions  from  the  exterior  region  Green’s  functions  described 
*  in  Section  4.3.  The  subscript  describes  which  potential  quantity  the  Green’s 

function  represents. 

»  4.3  Exterior  Region  Green’s  Functions 

The  remaining  set  of  Green’s  functions  needed  in  this  analysis  gives  the  fields 
radiated  into  the  half  space  y  >  h  by  magnetic  currents  and  charges  impressed 
on  a  PEC  ground  plane  at  y  =  h.  The  partial  differential  equations  satisfied 
by  these  potentials  are  identical  to  (11)  and  (12)  if  e  and  k  are  replaced  by 
their  free-space  values  Co  and  ko,  respectively.  The  effect  of  the  PEC  is  easily 
accounted  for  by  using  image  theory.  The  resulting  equations  for  the  exterior 
region  Green’s  functions  are 

p-jk 0R 

G?‘(  r;z',A,2')=—  (55) 

G%( r;  h,z')  =  G%{ r;  re',  h,  z')  =  -  ■  (56) 

where,  as  usual,  R  is  the  distance  between  source  and  observation  points  which 
in  this  case  is  _ 

R  =  \J{?-  s')2  +  {y  “  *)2 z')2.  (57) 

5  Moment  Method  Formulation 

We  now  make  a  slight  change  of  notation  from  the  previous  sections.  Since 
the  only  currents  to  be  treated  are  surface  currents,  we  will  dispense  with  the 
subscript  “s”  used  to  distinguish  surface  density  functions  from  volume  density 
functions.  Thus,  for  example,  the  symbol  “J”  will  be  used  hereafter  in  place 
of  “Js”  and  will  be  understood  to  represent  an  electric  surface  current  with 
units  of  A/m.  V”  will  be  used  henceforth  to  represent  surface  electric  charge 
density  with  units  C/m2.  Similar  remarks  hold  for  the  magnetic  sources. 

■>  Also,  recall  that  Mext  =  -Mint  and  cr™1  =  Therefore  we  formulate 

all  our  equations  in  terms  of  only  M‘nt  and  a}JJl,  and  the  superscript  “int”  will 
be  dropped  from  these  quantities  in  the  remainder  of  this  note. 

5.1  Integral  Equations 

The  equations  determining  the  electric  and  magnetic  surface  currents  are 

-(xx  +  zz)  •  Es(r)  =  (xx  +  zz)  •  E'(r),  r  6  support  of  J  (58) 
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(xx  +  zz)  •  ^Hmt(r)  -  Hcxt(r)J  =  -(xx  +  zz)  •  H’(r), 


r  e  support  of  M. 

(59) 


In  (58)  E'  is  the  incident  electric  field  radiated  by  the  impressed  electric 
current  source  at  z  =  0+,  and  E4  is  the  scattered  electric  field  between  the 
plates  given  by 

E4  =  -ju A  -  V$  +  e_1VxFint.  (GO) 

A,  $,  and  Fint  in  (60)  are  found  by  convolving  the  appropriate  Green’s  func¬ 
tions  and  sources. 


A(r)  = 

JJ  {xJI(r,)G/iI(r;  r')  +  zJg(r')GAt(r\  r#)}  dS' 

v'=o 

(61) 

*00  = 

JJo(r')G<b(r-,  r')  dS' 

(62) 

Fmt(r)  = 

y  — u 

J [  {xMx(r')<jjp‘(r;  r')  +  zM^Gf^t]  r')}  dS'. 

y'=h 

(63) 

The  function  H'  in  Equation  (59)  is  the  incident  magnetic  field  radiated  by 
the  impressed  electric  source  at  z  =  0+  with  the  notch  aroa  filled  in  by  PEC. 
The  remaining  (scattered)  magnetic  fields  in  Equation  (59)  are  expressed  in 

terms  of  potential  functions. 

jjint  _  jjext  _  juF  _  V$  +  /(J’VXA. 

(64) 

The  new  functions  used  on  the  RHS  of  (64)  are 

«(r)  = 

®inl(r)  -  rxt(r)  =  JJ  am(r')Gy(v,r')  dS' 

(65) 

F(r)  = 

y  =ti 

Fint(r)  -  Fext(r) 

~ 

Jj 1  {xMx(r')(?FI(r;r/)  +  zM*(r')(?F,(r;r')}  dS' 

yr=h 

(66) 

where 

G«(: r;0  =  G^(r;iO  +  G*»“‘V;r') 

(67) 

Gf,( t;t')  =  Ggf: r;r')  +  G?>;r') 

(68) 

Gl-,(r:r')  =  Gg(r;r')  +  GKl(rir'). 

(69) 
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5.2  Basis  Functions 


According  to  the  method  of  moments,  the  unknown  electric  and  magnetic 
currents  are  expanded  in  a  set  of  known  basis  functions  with  unknown  coef¬ 
ficients,  which  are  then  solved  for  by  enforcing  equations  (58)  and  (59)  in  an 
approximate  manner.  We  use  the  triangle  subdomain  basis  functions  of  Rao 
et  a/.[12],  which  are  particularly  suitable  for  approximating  arbitrarily  shaped 
regions.  We  proceed  by  partitioning  the  support  of  J  and  the  support  of  M 
into  triangles.  A  basis  function  is  then  defined  over  each  pair  of  triangles 
which  share  a  common  edge.  Figure  2  shows  two  such  triangles,  T+  and  T~, 
which  share  edge  n  of  the  triangulated  surface.  Points  in  T*  may  be  desig¬ 
nated  either  by  the  position  vector  r,  which  locates  them  with  respect  to  the 
global  origin,  or  by  which  is  defined  with  respect  to  the  free  vertex  of  T+. 
p~  is  defined  similarly  for  points  in  T~,  except  that  it  is  directed  towanl  the 
free  vertex  of  T~.  Which  of  the  two  triangles  is  designated  “plus”  and  which 
is  “minus”  depends  on  the  choice  of  positive  current  reference  direction  for 
the  nth  edge.  Positive  reference  current  flows  across  the  edge  from  T*  to  T~. 
The  basis  function  associated  with  the  nth  edge  is  then  defined  as 


fn(r) 


fn  4. 

2 AtPn 

In  _ 

2  AnPn 


ifr  €T+ 
if  r  G  T~ 


0  otherwise 


(70) 


where  /„  is  the  length  of  the  nth  edge,  and  A*  is  the  area  of  triangle  T*.  The 
properties  of  these  basis  functions  are  discussed  in  [12].  The  unknown  electric 
and  magnetic  currents  are  expanded  in  series  of  these  basis  functions. 


NM 

M(r)  =  EAUfV)  (71) 

n=l 

Nj 

J(r)  =  £j„tf(r).  (72) 

n=  1 

The  f'Mictions  |f,f  j  and  |f are  as  defined  in  (70),  but  the  former  have 
their  support  in  the  plane  y  =  0  while  the  latter  are  nonzero  only  in  the 
plane  y  =  h.  The  number  of  electric  current  basis  functions.  Nj.  is  equal  to 
the  number  of  interior  edges  of  the  triangulated  surface  in  the  y  =  0  plane 
(plus  edges  at  y  =  z  =  0).  Similarly,  the  number  of  magnetic  current  basis 
functions,  Nm,  is  equal  to  the  number  of  interior  edges  in  the  triangulation 
of  the  notch  area  in  the  y  -  h  plane. 
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Figure  2:  Triangular  basis  function  geometry  showing  two  triangles  with  a 
common  edge  (from  [12]). 

5.3  Reduction  to  Matrix  Equation 

We  now  define  an  inner  product  for  use  in  testing  equations  (58)  and  (59)  to 
be 

(f,g)  =  JJf-gdS.  (73) 

We  choose  our  testing  functions  to  be  the  same  as  our  basis  functions  (Galer- 
kin’s  method).  Then  taking  the  inner  product  of  (58)  with  each  basis  function 
in  turn  yields 

m  =  1,2,...  ,1V/.  (74) 

Following  [12]  we  can  approximate  the  second  term  on  the  L1IS  of  (74)  by 

(V*,£)  «  -lm  [#(#)  -  $(C)]  (75) 

where  i\~  is  the  centroid  of  triangle  .  The  remaining  terms  are  also  ap¬ 
proximated  as  in  [12]. 

(A,#)  «  y  [A(r^+)  •  p7c+  +  A(C)  •  Pc~x]  (76) 
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1 .  OF'"1  , 


‘m 

T 


.  c/Fint  c+. 

yxTr<^) 


npint 

Pm  +yx-^— (O'Pr,: 


% 


(77) 


pc-~ 

rm 


is  the  vector  p*  evaluated  at  the  centroid  of  the  corresponding  triangle. 
Returning  now  to  Equation  (59),  we  again  compute  the  moment  with  each 
basis  function  in  turn. 

3«(F,f")-(v*,f")  +  ^‘(yx|^,f")  =  (-H'.f"), 

m=  1,2,...,!^.  (78) 

We  approximate  these  inner  products  in  exactly  the  same  way  as  before. 


(v®,f«}«-/„!*(rS-)-*(C)) 
(F.f'1)  »  If  |F(r;f )  •  pC  +  F(C )  •  Pt) 


a  d  A  fcf 

yx  vf- 


J  yX|^(rm  )  •  Pm  +  yX^(C)  •  Pm 


(79) 

(80) 
(81) 


Note  that  7?r  in  (75),  (76),  and  (77)  refers  to  the  common  edge  of  the  triangles 
comprising  the  support  of  the  mth  electric  current  basis  function,  while  the 
same  index  was  used  in  (79),  (80),  and  (81)  to  designate  an  edge  in  the 
support  of  the  ?ntl1  magnetic  current  basis  function.  This  abuse  of  notation 
should  cause  no  confusion  since  m  will  always  be  used  consistently  in  any 
given  expression,  and  the  meaning  should  be  clear  from  the  context. 

After  substituting  the  expansions  of  the  currents  into  (74)  and  (78)  we 
obtain  the  desired  matrix  equation 


AX  =  B. 


(82) 


A  is  the  so-called  interaction  matrix,  a  square  matrix  of  order  Nj  +  Nm •  X 
is  the  column  matrix  of  unknown  current  coefficients  and  B  is  the  excitation 
matrix.  It  is  illuminating  to  write  these  matrices  in  partitioned  form,  as 
follows. 


(83) 


where  £  =  [£mj  and  J  =  {Jm)  are  column  matrices  of  iength  Nj,  ’’H  =  [Hm j 
and  M.  -  [Mm)  are  column  matrices  of  length  Nm,  and  W  =  [>Vmn],  U  = 
[Z/mn],  y  =  [Tmn]>  and  Z  =  [Zmn]  are  matrices  of  dimensions  Nj  x  Nm,  Nm  x 


■  y 

u  ’ 

’  M  ‘ 

'  H  ' 

w 

z 

J 

£ 
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Nj,  iVjvf  X  NM,  and  Nj  x  iVj,  respectively.  The  elements  of  the  interaction 
matrix  are  given  below. 

Wmn  =  --  fyx9Pt.  ■  # + yx8F;m  ■  (84) 


where 

8F*„  =  //f"(p')'G8F«fy)«' 

=  'f  -jj  JJpi-Gs,.i^-<r')dS'+jz  //p:  ■  a, Hr?; 


and  the  dyadic  Green’s  function  G#/?  is  defined  by 

GaHr;r')  =  xx^(r,0  +  z^(r;r'). 


2mn  —  lm  &  (  A*m  •  ~~  +  Amn  '  2  ^>mn  ^mn 


where 


A*„  =  //  t^r')-6A(^,r'),IS' 

=  |  -j  JJpi-GA(ri-y)dS’+j,  jJp;-GArt-y)dS’ 


=  i  JJ  x>:-t^v')Gtt,v£y)ds' 


-  V  3F  JJw&W-k  //<**«')«' 
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G/i(r;  r')  =  xxGax{v,  O  +  zzGa,{  r;  r') 
and  we  have  used  the  fact  that 

4  if  r  G  T+ 

At 

Vs-f„(r)=  |  -A:  ifrGT" 

/ln 

0  otherwise. 


y,nn  =  lm  (Fmn  ‘  +  F»  '  J  + 


where 


F*„  =  //  f"(r') •  GF(r? ,r')<iS' 

=  I  3f  r(r?;r')«' 

[  "  ri  T ■ 

$±n  =  l  JJ  V's-  f*/(r')G«i-(r^,r,)<f5' 

=  JjG^^dS'-jrJjG^i^dS' 

^  ,  rn  — 


where 


G/Kr;  r')  =  xxG>,(r;  r')  +  zzGV^r;  r'). 

W,n»  =  ~  yX^A+„  •  4"  +  yxdA, mn  ■  2 

L  “ 

=  JJ  f,f(r')  •  Ga^fr^f ,  r')  dS' 
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in 

2 


//rf-e»j(r?;r')rfS'+-i-  JJp: -GeArf-, r')dS' 

T,t  U  T,7 


(97) 


G$a(*\  r')  =  xx^(r;  r')  +  zz^^-(r;  r7).  (98) 

5.4  Efficient  Evaluation  of  Interaction  .Matrix  Entries 

As  discussed  in  reference  [12],  efficient  evaluation  of  interaction  matrix  entries 
requires  consideration  of  face-pair  contributions  rather  than  edge-pair  contri¬ 
butions.  The  face-pair  integrals  are  of  two  basic  types:  vector  and  scalar.  We 
describe  here  how  to  efficiently  calculate  the  entries  in  Z.  Similar  techniques 
are  used  for  the  other  partitions.  Consider  first  a.  typical  vector  integral  needed 
to  evaluate  Zmn.  The  quantity 

Af'  =  ~  JJ Pi  ‘  G'/i(rcp;r7)  dS'  (99) 

yq 


represents  the  magnetic  vector  potential  at  the  centroid  of  triangle  Tp  due 
to  the  ith  basis  function  defined  on  triangle  Tq.  In  (99)  we  have  temporarily 
assumed  a  local  numbering  scheme  on  the  source  triangle,  with  vertices  at  iq, 
i’2,  and  r3,  and  where 


Pi  =  ±(r7  -  r,-),  i  =  1,2,3.  (100) 


The  plus  sign  above  is  selected  if  the  associated  basis  function  defines  positive 
current  to  flow  out  of  the  tth  edge  of  the  triangle,  which  lies  opposite  to  the 
ith  vertex. 

One  problem  in  evaluating  (99)  is  that  the  integrand  is  singular  when 
triangles  p  and  q  are  the  same.  The  singularity,  which  occurs  in  G^,  is  the 
same  as  that  of  the  free-space  Green’s  function  (multiplied  by  2  to  account 
for  the  image  in  the  PMC  plane  in  the  bilateral  case) 


G  4fr  v'\  — 

’  7  Anil 


f  1  A  »  V 

Ovi; 


where  I  is  the  dyad  xx  +  zz  and  where  s  =  1  in  the  unilateral  case  and  s  =  2 
in  the  bilateral  case.  There  are  several  possible  ways  of  handling  this  singu¬ 
larity.  The  first  is  to  explicitly  extract  and  integrate  the  contribution  of  the 
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singularity  in  closed  form.  The  second  is  to  use  an  adaptive  cubature  routine 
which  is  capable  of  accurately  treating  integrable  si ngulari ties[  1 3] .  Although 
the  latter  method  is  simpler  to  implement,  it  is  far  more  computationally 
expensive,  requiring  hundreds  of  integrand  evaluations  to  achieve  only  about 
0.1%  relative  accuracy.  Therefore,  we  continue  now  with  a  description  of  the 
first  method.  It  is  convenient  for  the  later  numerical  procedures  to  define  the 
regular  part  of  the  dyadic  Green’s  function  as 


r>r  _  r<  sl‘°  x 

Sir 

We  can  then  write 

/  \ 

(102) 

JSjdS' 

y  Jq 

► .  (103) 

The  first  integral  in  (103)  is  well  behaved  in  all  cases  and  can  be  integrated 
numerically,  using  a  low  order  Gaussian  scheme  designed  especially  for  inte¬ 
gration  over  triangular  domains  (see  (14),  Table  8.3).  The  second  integral  is 
singular,  but  can  be  integrated  analytically  in  closed  form(15).  These  evalua¬ 
tions  are  examined  in  detail  in  the  following  paragraphs. 

The  dyadic  in  Equation  (103)  contains  only  two  nonzero  elements,  G\ 
and  GAt,  which  are  defined  by 

G"i  =  xxGrAx  +  ZzGAr 

(104) 

The  desired  vector  integral  can  thus  be  written  as  the  sum  of  two  scalar 
integrals. 

jj Pi  ■  G,i(rcp;  r')  (IS*  =  xlix  +  zliz  (105) 

'VO 

I  ■ 

where  the  scalar  integrals  are  given  by 

Ux  =  JJ(x- Pi)  GrAx(vcp\ r')  dS' 

Tn 

(106) 

x  9 

liz  =  Jj(z-Pi)GrAi(tcp-y)  dS'. 

Tq 

(107) 

The  integrals  over  the  triangular  regions  are  most  efficiently  accomplished 
by  using  normalized  area  coordinates(12,14]  £  and  i].  The  source  point  can  be 
written  in  terms  of  area,  coordinates  as 


r'  =  £r2  +  ??r3  +  (1  -  £  -  rj)ri  (108) 
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where  it  is  recalled  that  the  {r,}  are  the  vertices  of  the  source  triangle  T7. 
Similarly, 

Pi  =  ±(r'  -  r.)  =  ±  (fr2  +  7?r3  +  (1  -  ^  -  j/)n  -  r,) .  (109) 

The  area  element  in  normalized  area  coordinates  is  d,S'  =  2A7  d£  dr/. 

Let  t  stand  for  either  x  or  z.  Then  the  integrals  In  can  be  expressed  in 
terms  of  integrals  which  are  independent  of  i  in  the  following  way. 

Iti  =  ±t  •  [r2It(  +  r 3/(r)  +  rju;  -  r,/«]  (110) 

where 

Iti  =  2Aq  f  f  £GA  (rcp;£r2  +  J?r3  +  (1  -  £  -  ?/)ri)  di)  (111) 

Jo  Jo 

hrt  =  2  Aq  [  [  vGrAt(vcp-,^r2  +  'i]T3  +  (l-^-ij)r1)  d£  dr/  (112) 

JO  Jo 

Iti  =  2Aq  [  i  (l-^-7?)G^,(rcp;^r2  +  7/r3  +  (l-^- J?)ri)  dr) 

(113) 

It  =  2 Aq  [  [  Gta ,(rcp;fr2  +  ??r3  +  (l-f-7?)ri)  d£dr}.  (114) 

Jo  Jo 

However,  the  above  integrals  are  not  all  independent.  We  also  have  that 

It  =  Iti  +  +  hi  (H5) 

so  that  in  general  there  are  six  integrals  to  be  numerically  evaluated  for  each 
observation  point  /  source  triangle  pair.  These  integrals  contribute  to  as  many 
as  nine  entries  of  2,  depending  on  the  number  of  basis  functions  defined  in 
each  triangle.  They  are  computed  using  the  formula  from  [14], 

2  /  /  /(£>*/)  d£  dr)  «  (116) 

Jo  Jo  tel 

The  sample  points  and  weights  are  given  in  Table  2.  As  mentioned  previously, 
the  second  integral  in  (103)  is  evaluated  using  the  formulas  of  [15]. 

In  addition  to  the  vector  integrals  that  were  just  discussed,  we  also  must 
calculate  scalar  integrals  of  the  form 

v?=±£kJJG*(r°,rtdS'-  <U7> 

jq 
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k 

& 

Vk 

Wk 

0.53333333 

0.33333333 

0.22500000 

0.05971587 

0.47014206 

0.13239415 

0.47014206 

0.05971587 

0.13239415 

4 

0.47014206 

0.47014206 

0.13239415 

5 

0.79742699 

0.10128651 

0.12593918 

6 

0.10128651 

0.79742699 

0.12593918 

0.10128651 

0.10128651 

0.12593918 

* 


Table  2:  Sample  Points  and  Weights  for  Seven  Point  Triangular  Cubature 


<I,f?  is  the  electric  scalar  potential  at  the  centroid  of  triangle  Tp  due  to  the 
surface  electric  charge  density  associated  with  the  /ll‘  basis  function  defined 
on  triangle  Tq.  As  before,  the  sign  is  taken  to  be  positive  if  the  basis  function 
in  question  defines  positive  reference  current  to  flow  out  of  T'1  through  the 
edge.  Very  little  additional  work  is  needed  to  compute  this  integral,  since 
G',tx  =  /ioc6V  Using  this  fact,  we  can  then  write  (117)  as 


«J >?«  =  ±-  Hi 

'  i^nocA'i 


T< 


s/(p  r  r  d&_ 

4tT  JJ  R 


f  \ 


(118) 


Computation  of  Ix  has  already  been  described,  and  computation  of  the  last 
integral  above  is  accomplished  using  the  formulas  of  [15].  We  note  again  that 
the  quantity  d>f7//;  is  common  to  as  many  as  nine  entries  of  Z. 


5.5  The  Effective  Excitation  Matrix 

The  excitation  matrix  entries  are 

£m  =  (E'.fi)  (119) 

Hm  =  (-HSf*/).  (120) 

First,  since  the  impressed  electric  current  radiates  only  evanescent  waves  into 
the  parallel  plate  region,  our  assumption  that  the  slot  is  far  from  the  origin 
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implies  that  'Hm  =  0  for  all  m.  Second,  the  impressed  currents  are  transverse 
to  2,  the  direction  of  the  induced  currents  near  the  origin.  Thus,  the  incident 
magnetic  vector  potential  makes  no  contribution  to  (119).  The  only  remaining 
contributor  to  (119)  is  the  incident  electric  scalar  potential,  whose  source  is 
impressed  charge  density.  One  can  show  that  there  is  no  impressed  surface 
charge  at  2  =  0,  leaving  only  the  impressed  line  charge  at  2  =  y  =  0  which 
is  due  to  the  discontinuity  in  Jy  there.  Because  the  impressed  and  induced 
currents  must  be  continuous  at  2  =  0,  it  is  also  easy  to  show  that  there  is 
an  induced  line  charge  at  2  =  y  =  0  which  is  exactly  the  negative  of  the 
impressed  one.  The  sum  of  the  scalar  potentials  due  to  these  two  sources  is 
exactly  zero.  Therefore,  we  may  choose  not  to  compute  either  impressed  or 
induced  potentials  due  to  these  sources  when  filling  our  various  matrices,  as 
long  as  we  force  the  induced  currents  near  z  =  0  to  assume  their  prescribed 
values.  We  are  then  left  with  a  homogeneous  matrix  equation  where  some  of 
the  unknowns  (those  corresponding  to  currents  crossing  the  2  =  0  plane)  are 
actually  known. 

We  now  partition  the  “unknowns”  in  a  different  way  than  was  done  in 
the  previous  section.  There,  the  partition  was  defined  in  a  way  convenient 
for  filling  the  matrices.  Here,  we  partition  the  unknowns  according  to  the 
physical  location  of  the  corresponding  basis  function.  We  number  the  electric 
current  basis  functions  so  that  those  corresponding  to  triangle  edges  at  2  =  0 
are  numbered  last.  Then  X\  is  the  column  matrix  of  all  magnetic  current 
coefficients  and  all  electric  current  coefficients  for  edges  not  at  2  =  0.  X 2  is 
the  matrix  of  coefficients  for  electric  current  basis  functions  defined  for  edges 
at  2  =  0.  If  we  partition  the  remaining  matrices  in  a  conformable  fashion  then 
the  matrix  Equation  (82)  becomes 


The  matrix  X2  is  known.  The  set  of  remaining  unknowns,  X\,  can  be  found 
by  solving 

AUX 1  =  -13'  (122) 

where  13'  =  A12X2. 

5.6  Scaling  and  Matrix  Condition  Number 

The  matrix  equation  (122)  exhibits  very  large  condition  numbers.  This  prob¬ 
lem  can  be  traced  to  the  fact  that  the  equations  and  unknowns  are  each 
expressed  in  two  different  sets  of  units.  To  insure  that  all  unknowns  have 
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units  of  A/m  we  define  the  auxiliary  unknown 


M'  =  —M. 

(123) 

Vo 

Our  set  of 

matrix  equations  then  becomes 

y'M'  +  U'J  =  W 

(124) 

WM'  +  ZJ  =  £ 

(125) 

where 

y  =  vly<  w'  =  7/ow,  w  =  % u,  w  =  % n. 

(126) 

In  this  way  all  matrix  entries  have  units  of  ft-ni.  These  equations  exhibit 
reasonable  condition  numbers  and  are  amenable  to  accurate  solution.  After 
solving  them,  M.  is  computed  from  M!  using  (123). 


6  Stripline/Slotline  Transition  Design 

A  crucial  element  in  the  design  of  a  notch  radiator  is  the  stripline/slotline  tran¬ 
sition.  The  transition  performance  can  be  the  limiting  factor  in  determining 
impedance  bandwidth.  Figure  3  shows  a  schematic  diagram  of  a  typical  test 
fixture  used  for  measuring  transition  performance.  It  consists  of  two  tran¬ 
sitions  arranged  back-to-back  to  facilitate  easy  measurement  using  stripline 
connectors.  Swept  frequency  measurements  of  return  loss  and  insertion  loss 
for  this  symmetrical  two-port  provide  information  on  the  useful  bandwidth  of 
the  transitions. 


We  use  the  usual  technique  of  twice  analyzing  half  the  structure,  placing 
first  PMC  and  then  PEC  walls  at  the  symmetry  plane.  If  P^./  and  F \?  are  the 


23 

499 


reflection  coefficients  calculated  in  the  presence  of  the  magnetic  and  electric 
wall,  respectively,  then  the  two-port  scattering  parameters  are  given  by 


5ii  =  522  =  i(rA/  +  rB) 

(127) 

S12  =  =  |(rA/  -  rE)- 

(128) 

The  only  changes  needed  to  implement  this  capability  are  the  addition  of 
two  more  image  points  to  the  Green’s  functions  and  defining  unknowns  for 
edges  along  the  symmetry  plane  in  the  notch  triangulation.  These  unknowns 
are  set  to  zero  in  the  PEC  case. 

7  Results 

In  this  section  we  compare  program  predictions  with  theoretical  and  measured 
results.  All  the  cases  reported  here  assume  er  =  2.2  and  h  =  0.015  inch.  In 
all  cases  when  a  notch  is  present,  it  is  bilateral. 

7.1  Uniform  Width  Strip 

The  simplest  test  case  is  a  uniform  width  stripline  without  any  slot  in  its 
groundplane.  The  strip  dimensions  are  0.022  by  1.5  inches.  The  program  was 
run  using  a  coarse  and  a  fine  triangulation  with  Nj  =  30  and  60,  respectively. 
The  program  predictions  were  compared  to  simple  transmission  line  theory 
assuming  that  the  line  was  terminated  in  a  perfect  open  circuit,  neglecting 
any  end  effects.  The  magnitude  and  phase  of  the  input  reflection  coefficient 
were  computed  at  0.2  GHz  intervals  between  1  and  20  GHz.  The  computed 
amplitude  was  equal  to  1  (to  4  decimal  places)  for  both  triangulations  at  all 
frequencies.  The  difference  between  computed  and  theoretical  phases  was  at 
most  56.0  and  8.1  degrees  for  the  coarse  and  fine  triangulations,  respectively. 

For  the  fine  triangulation  computing  the  entries  in  the  fill  matrix  required 
14.5  seconds  per  frequency  and  solving  the  matrix  equation  using  double  pre¬ 
cision  Gaussian  elimination  required  0.59  seconds  per  frequency.  Since  matrix 
fill  time  dominates  the  cost  of  the  procedure,  we  implemented  the  interpola¬ 
tion  technique  of  Newman[16].  Under  this  scheme  the  matrix  is  filled  at  some 
fixed  number  of  logarithmically  spaced  frequencies  (called  knots)  covering  the 
band  of  interest.  The  matrix  is  reconstructed  for  intermediate  frequencies  by 
quadratically  interpolating  the  matrix  values  at  the  closest  three  knots.  Using 
this  technique  only  10  knots  were  required  to  cover  the  20:1  frequency  band. 
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0  10  20 
Frequency  (GHz) 


Figure  4:  Reflection  coefficient  phase  error  for  fine  (o)  and  coarse  (*)  triangu¬ 
lations  of  the  stepped-width  stripline.  Knots  are  marked  with  bullets. 


Figure  5:  Reflection  coefficient  phase  error  for  fine  triangulation  using  10  (o) 
and  19  (*)  knots.  Knots  used  in  the  19  point  scheme  are  marked  with  bullets. 
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7.2  Stepped-width  Strip 

The  next  case  consists  of  a  stripline  witli  a  step  in  width.  As  in  the  previous 
case,  there  is  no  notch,  and  the  main  line  width  is  0.022  inch.  At  a  point  one 
inch  from  the  feed  port,  the  line  increases  in  width  to  0.044  inch.  The  length  of 
the  wide  section  is  O.f)  inch.  Again,  two  triangulations  were  examined,  using  45 
and  111  unknowns.  Again  the  reflection  coefficient  magnitude  was  correctly 
computed  to  be  unity  over  the  band  for  both  cases.  The  phase  prediction 
was  compared  to  a  Touchstone  circuit  model  and  the  error  with  respect  to 
this  model  is  plotted  in  Figure  4.  Ten  knots  were  used  in  the  matrix  fill 
interpolation  procedure  for  both  cases.  It  is  seen  that  the  error  generally 
increases  with  frequency  and  the  fine  triangulation  is  superior  to  the  coarse 
one.  Jumps  in  the  plotted  values  occur  between  knots  where  a  change  in 
frequencies  used  within  the  interpolation  scheme  occurs.  By  increasing  the 
number  of  knots  the  magnitude  of  these  jumps  can  be  reduced.  This  can  be 
seen  in  Figure  5  where  10  and  19  knot  interpolation  schemes  are  compared 
for  the  stepped-width  strip. 


-0. 


0 

X  (IN) 


0.  25 


Figure  6:  Detail  of  the  triangulation  used  to  model  a  stepped-width  strip 
feeding  a  rectangular  notch  in  it  groundplanes.  Uniform  sized  triangles  are 
used  on  the  portions  of  the  structure  not  shown. 


7.3  Stripline-fed  Notch 

Here  a  stepped-width  stripline  feeds  a  rectangular  notch  as  shown  in  Figure 
6.  The  notch  measures  2  by  0.06  inches.  The  0.022  by  1.25  inch  uniform  strip 
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Figure  7:  Comparison  of  measured  (solid)  and  computed  (dotted)  return  loss 
for  the  stripline-fed  rectangular  notch  of  Figure  6.  Knots  are  marked  with 
bullets. 

is  terminated  with  a  0.15  inch  square  stub.  Notice  that  the  width  step  occurs 
directly  under  one  edge  of  the  notch,  adding  to  the  difficulty  in  modeling  this 
transition. 

A  test  fixture  was  built  and  measured  to  provide  a  comparison  with  the 
numerical  results  of  the  program.  Measured  and  computed  return  loss  are 
displayed  in  Figure  7.  For  this  case  we  had  Nj  =  105  and  Nm  =  79.  Fill  time 
and  matrix  equation  solution  time  were  96  and  10.4  seconds  per  frequency, 
respectively.  The  good  agreement  between  measured  and  computed  results 
verifies  the  accuracy  of  the  theory. 

7.4  Two-port  Transition  Test  Fixture 

A  two-port  test  fixture  of  the  type  shown  in  Figure  3  was  constructed  and 
measured.  Good  agreement  was  obtained  for  insertion  loss  measured  and 
calculated  between  the  two  ports.  These  results  will  be  presented  at  the 
symposium. 

8  Conclusion 

An  efficient  technique  has  been  developed  to  analyze  arbitrarily  shaped,  strip- 
line-fed  notches.  Measurements  verify  the  accuracy  of  the  method.  A  simple 
extension  of  the  analysis  gives  a  useful  tool  for  designing  stripline/slotline 
transitions. 
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A  TECHNIQUE  TO  MEASURE  LARGE  ANTENNA  ARRAYS 
THAT  CONTAIN  RANDOM  AMPLITUDE  ERRORS. 

Peter  R.  Franchi 
Harvey  Tobin 

Hanscom  Air  Force  Base,  Massachusetts  01731 


ABSTRACT: 

The  measurement  of  two  dimensional  planar  arrays  that  are 
electrically  long  is  currently  a  difficult  task.  The  accurate 

measurement  of  the  antenna  pattern  is  limited  by  large  antenna 

range  requirements  for  far  field  ranges,  the--  number  of  sample 

points  that  must  be  measured  for  near  field  characterization,  and 

by  the  great  size  and  very  high  tolerance  of  the  reflector  and 

feed  in  a  compact  range.  At  the  1988  Allerton  Symposium  we 

presented  a  paper  that  described  how  to  correct  for  the  error 

level  of  the  far  field  patterns  using  an  ideal  amplitude 

distribution.  This  paper  is  a  follow  on  to  the  1988  presentation. 

It  describes  how  to  measure  and  correct  for  error  levels  of  a 

large  aperture  antenna  that  contains  random  errors  in  amplitude 

and  phase  by  applying  a  quadratic  and  higher  order  phase 

correction  to  the  phased  array. 

The  errors  that  result  from  this  approach  are  analyzed  as  a 
function  of  range,  antenna  electrical  length,  sidelobe  level  and 
degree  of  correction.  Both  computer  simulations  and  analysis  are 
used  to  examine  the  limitations  of  this  approach. 

INTRODUCTION: 


In  the  future  the  United  States  Air  Force  will  be  mounting 
large  lightweight  phased  arrays  on  the  fuselages  of  long  aircraft. 
The  size  of  the  antenna  arrays  considered  are  approximately  60A  to 
300A  in  length.  Todays  conventional  antenna  ranges  are  not 
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sufficient  to  measure  accurately  the  performance  of  low  sidelobe 
antennas  in  the  far  field.  This  is  based  on  the  accepted  2D  /A 
criterion.  For  example,  the  far  field  distance  for  a  10' x  100' 
aperture  at  C-band  would  be  23  miles  away.  If  low  sidelobes  are 
required,  then  the  distance  would  be  even  greater  to  prevent 
distortion  within  the  far  field  patterns.  To  use  compact  ranges 
would  require  a  massive  high  tolerance  reflector,  and  the 
necessary  near  field  ranges  to  house  the  reflectors  would  be 
extremely  expensive. 

This  paper  discusses  a  method  to  help  overcome  the  limitations 
of  conventional  antenna  ranges.  For  a  typical  two  dimensional 
planar  array  shown  in  Fig.l,  to  measure  the  far  field  performance 
in  the  smaller  dimension  the  antenna  can  be  considered,  to  be  in 
the  far  field.  For  the  larger  dimension  the  antenna  can  be 
considered  to  be  in  the  region  of  the  Fresnel  field  resulting  in 
apparent  far  field  pattern  distortion.  An  accurate  far  field 
measurement  of  the  antenna,  however  can  be  reclaimed  from  this 
Fresnel  region  measurement  by  putting  in  a  phase  correction  over 
the  aperture. 

At  the  1988  Allerton  Symposium  we  presented  a  paper  that 
considered  a  undisturbed  50db  Taylor  amplitude  distribution,  and 
explained  how  to  approximate  the  far  field  pattern  employing  repeatec 
phase  corrections  of  quadratic  and  cubic  terms. 

This  report  is  an  extension  of  the  88  Allerton  report 
considering  how  random  amplitude  errors  added  to  the  50db  Taylor 
distribution  affect  the  validity  of  the  quadratic  phase 
correction.  These  errors  have  been  inserted  by  perturbing  the 
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amplitude  distribution  in  a  statistically  random  way  with  a  zero, 
mean  and  fixed  variance.  Random  phase  errors  were  not 

employed  since  the  random  amplitude  errors  cause  typically  the  same 
pattern  distortions.  These  amplitude  errors  would  be  typical  of 

*  beamformer  network  errors  or  measurement  errors  but  do  not 
represent  errors  caused  by  coarsely  quantized  phase  shifters. 

f  Basically  the  phase  shifter  setting  errors  here  are  considered 

much  smaller  than  the  beamformer  network  errors, 

Finally,  we  must  define  sidelobe  levels  that  are  caused  by 
different  sources  and  which  can  be  separated  for  low  sidelobe 
arrays.  The  theoretical  or  ideal  sidelobe  were  treated  in  the 
previous  report.  Error  sidelobes  result  from  amplitude (and/or 
phase)  error  in  the  aperture  illumination  function  due  to  network 
errors.  In  this  paper  ideal  sidelobes  are  considered  to  be  set 
much  lower  than  the  error  sidelobe  levels,  and  do  not  affect  the 
sidelobe  considerations.  The  error  sidelobes  are  the  principle 
concern  of  our  measurement  technique.  Measurement  of  the  pattern 
in  the  Fresnel  zone  results  in  sidelobes  that  are  incrementally 
higher  than  the  error  sidelobes.  This  incremental  increase  is 
called  add-on  error.  In  order  to  obtain  an  accurate  representation 
of  the  far  field  pattern  with  a  Fresnel  zone  measurement  we  must 
keep  this  add  on  error  very  small.  The  increase  in  sidelobe  level 

* 

can  be  held  within  an  acceptable  bound  by  applying  the  phase 

^  correction  as  many  times  as  needed. 

DISCUSSION: 

In  the  1988  Allerton  report  we  studied  how  to  measure  large 
antenna  arrays  by  applying  a  phase  correction  to  the  antenna. 
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The  results  showed  that  by  applying  a  second  order  correction  to 
the  phase  term  the  add-on  error  can  be  reduced  significantly 
the  sidelobe  levels.  We  also  showed  that  we  can  do  even  better  by 


applying  a  third  order  term. The  equations  are: 


2  2 
k  cos  8  D 


For  the  quadratic  term  (1) 


--- — For  the  cubic  term  (2. 
16R0 

where, 

k=2n/A 

A=free  space  wavelength 

D=aperture  width 

RQ=distance  from  tranmsit 
recieve  (test)  antenna. 

0Q=angle  of  scan  for  which  correction 
is  made. 

The  results  are  plotted  as  add  on  error  versus  R/A  for  several 

D  values,  and  the  number  of  phase  settings.  This  is  shown  in  Figure 

For  this  application  in  which  the  sidelobe  level  is  primarily 

determined  by  aperture  errors,  we  consider  a  low  sidelobe  antenna 

whose  geometry  is  given  by  Figure  3.  To  find  and  expression  of  R  w 

consider  an  arbitrary  triangle  and  from  the  distance  formula, 

R={  (b-x)2+(h-0)2}1/2  (3) 

squaring  this  term  we  obtain 

2  2  2 

R  =(b-x)  +h  (4) 

2  2  2  2 
R  =b  -2bx+x^+hz  (5) 

we  note  that; 

2  2  2 

RQ=b  +hz  (6) 


4 

510 


> 


2  2 

therefore  substituting  for  b  +h  , 

2  2  2 
R  =Rq+x  -2bx 


but, 


cosa=b/Rf 


b=RQCOSa 


(7) 


(8) 


A 


therefore  replacing  this  gives; 


R2=Rq+x2-2R0xcos<x 

then, 

R=[R2+x2-2RQxcosa] 1/2  (9) 

The  expression  for  the  field  in  the  Fresnel  region  of  a 

aperture  antenna  is  given  by, 

E(u)=C0X  f  (x)  ejkRdx  (10) 

R 

considering  R  much  grater  than  D  (10)  can  be  expressed  as, 

E(u)=CnX  f  (x)  ejkRdx  (11) 

R0 

where, 


CQ  is  arbitrary  constant 
1/Rq  is  the  amplitude  term 

e-)kR  is  the  phase  term 

f(x)  is  the  aperture  illumination  functio: 
Expanding  out  R  (in  Eq.9)  in  the  phase  term  we  obtain; 

R=Ro-xsine+x2/2Rocos20+. . . . .H.O.T  (12) 

The  term  we  are  most  concerned  with  is  the  quadratic  term, 
therefore  we  want  to  express  (10)  as, 

E(u)  =  T  f(x)e^k(ux^x2)dx  (13) 

where, 
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u=kdsine 

/3=kd2  cos2e 
SR' 

For  amplitude  and  phase  errors  in  f  (x)  (13)  can  be  written 


where, 


E  (u)  Cj/f  (x)  e' (“<x)  +35<x)  1  e3@x2e3uxdx  (14) 


a (x) -random  amplitude  error  in  f (x) 
5  (x) -random  phase  error  in  f (x) 


for  small  error, 

E  (u)  =/  f  (x)  [1- j/3x2]  e^UXdx  +  S  f  (x)  (cc(x)-j3(x)+ja(x)0x2 

-S(x) )ejuxdx  (15) 

the  first  term  is  the  ideal  theoretical  pattern. 

We  then  assume  that  the  deterministic  sidelobes  are  low  enough 

that  they  do  not  contribute  to  the  overall  pattern,  then  (15)  can 

be  expressed  as, 

E(u)=C1J'f  (x)  [-a  (x)  -  jS  (x)  +  ja  (x)  /3x2-3  (x)  /3x2]  e^UXd  (16) 
writing  this  expression  as  a  sum; 

g  (u)  =Zfn  [-an-jSn+jan0nAx2-5n+jan/3nAx2-5nPnAx2]e^unAx  (17) 


converting  this  into  power; 

P  (u)=  |E  (u)  E*  (u)  |  C12  I 

=C  2Zf  [-a  (1- j/3nAx2)  -  j6n  (1- j/3n2Ax2)  e^unAx 


*Zf v ( 1  +  j/3m2Ax2)  j/3m2Ax2 ) )  e" :umAx  (19) 


multiplying  this  out  we  obtain; 
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multiplying  this  out  we  obtain; 


=CfctVmV‘m(1+'!  1  1  ix  S  2x  23{n 


fmfn  (_a)  - D^2Ax2] e juAx  (n“m) 


(20) 


=C2Z[f  £2^l+e2n4Ax4)+f2£2  ^02n4Ax4)]  (21) 

P (u)  =C2Zf 2 («2n+«2n) +f2n  («2n+62) P2n4Ax4  (22) 


.  2  2. 

if  an=5nis  independent  of  n 


because  the  amplitude  and  phase  errors  do  not  contribute  any  type 


of  different  errors,  then  we  can  write  the  expression  as 

(23) 


22  22442  2 

P  (u)  =  (2a^+2a^Ax  In  f^) 


and  this  is  equivalent  to, 

P (u) =2a2Jf 2 (x) dx+2a2p2fx4 f 2 (x)  dx  (24) 


therefore  we  can  express  the  add  on  error  as, 

AE=  c2=  P  (u) 
cQ  PqTu) 


(25) 


2  2 
2cc*Zjr 
n  n 


1+  Ax402Zn4f2 


(26) 


Zf * 


n 


for  a  continuous  aperture, 


2  4  2 

1+  rs*  f  (x)dx 
Jf2  (x)  dx 


(27) 


The  add  on  error  is  now  defined  by  equation  (28)  relating  it 
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to  the  aperture  illumination  function  f  (x)  and  0  which  includes 
the  scan  angle,  aperture  size,  wavelength  and  Fresnel  range.  This 
equation  along  with  equations  (10)  and  (11)  allow  the 
experimentalist  to  set  acceptable  error  limits. 


SIMULATION: 

The  theoretical  results  of  the  last  section  represent  the 
average  behavior  of  the  add  on  error  when  an  infinite  set  of 
antenna  patterns  are  taken, corresponding  to  Gaussian  distributed 
amplitude  and  phase  errors.  In  an  actual  measurement  only  one 
pattern  or  a  very  small  number  of  sample  patterns  will  be 
taken  (for  one  scan  condition  or  frequency,  etc.).  This 

represents  the  basic  problem  of  characterizing  sidelobe  behavior 
when  phase/amplitude  errors  are  dominant.  In  other  words,  the 
actual  sidelobes  over  any  particular  angular  region  will  not 
represent  the  average  level  or  the  RMS  level.  The  simulation 
results  are  an  average  of  eleven  different  patterns  taken  for  the 
same  statistical  errors.  While  these  patterns  are  statically  the 
same  they  are  different  in  specific  sidelobe  behavior. 

Another  unavoidable  issue  is  what  is  the  best  method  of 
characterizing  the  error  sidelobe  levels.  In  all  cases  the  pattern 
with  and  without  quadratic  phase  correction  is  compared  with  the 
ideal  far  field  pattern.  Three  methods  seem  obvious.  One  is  a 
comparison  of  the  average  sidelobe  level  (in  power)  over  a  finite 
angular  region.  The  second  is  the  RMS  average  of  a  point  by  point 
difference (both  of  these  over  the  same  angular  region)  and  the 
third  is  a  comparison  of  the  changes  on  the  peak  sidelobes.  Since 
this  last  approach  seems  to  be  the  most  commonly  used 
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experimentalists,  it  is  used  here.  Peak  sidelobes  are  defined  by 
the  experimentalist  rule  of  thumb  of  5db  or  more  above  the  mean 
level  of  the  far  field  pattern.  Using  this  criterion  the  simulation 
gives  an  average  of  the  number  of  sidelobes  that  are  above  the  mean 
value  of  -47.5db  within  a  sector  of  width  eight  degrees  or  about 
.14  in  sine  space. 

Figure  4  is  the  far  field  pattern  of  the  array  with  the  first 
set  of  errors  over  4°  to  12°.  Figure  5  is  the  pattern  at  Ro=150' 
without  correction.  Figure  6  is  the  pattern  with  quadratic 
phase  correction  for  0q=O°.  It  is  quite  clear  that  with  the 
quadratic  correction  the  near  field  pattern  approximates  the  far 
field  pattern  well. 

Looking  at  cumulative  probability  curves  of  sidelobes  for  the 
same  cases  (percentage  of  sidelobes  above  a  given  sidelobe 
level)  there  is  a  similar  picture.  Figure  7  represents 
the  curve  of  the  far  field  distribution  for  the  pattern  in  Figure 
4.  Likewise  Figures  8  and  9  represent  the  curves  of  the  field 
distributions  of  Figures  5  and  6. 

Taking  all  eleven  data  sets  and  averaging  the  absolute  value 
of  the  add-on  error  gives  the  plots  in  Figure  10.  The  average  of  the 
eleven  simulations  follows  the  theoretical  behavior  for  the 
three  cases  but  it  is  clear  more  averaging  would  be  needed  for  a 
better  comparison. 

CONCLUSION: 

The  pattern  distortion  error  that  results  from  the  use 
of  phase  correction  in  measuring  very  large  arrays  in  their  near 
field  or  Fresnel  field  have  been  derived  theoretically.  These 
expressions  have  been  compared  to  ideal  pattern  disto.  from 
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simulations  in  the  previous  Allerton  paper  and  there  is  a  very 
close  match  and  to  error  antenna  patterns  in  this  paper  with  a 
less  accurate  but  quite  good  match.  It  should  be  possible  to  limit 
the  add  on  measurement  error  by  the  applications  of  these 
equations  and  applying  both  quadratic  and  cubic  phase  correction 
as  many  times  as  the  conditions,  Rq,D,X,  the  illumination  taper  and 
the  expected  phase  and  amplitude  error  call  for. 
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CUBIC  PHASE  ERROR  N=4  I  *  CUBIC  PHASE  ERROR  N= 6 


FIG.2  AMOUNT  OF  PHASE  SHIFTER  CHANGE 
N-PHASE  SHIFTER  SETTINGS 


FIG.5  PATTERN  W.OUT  CORRECTION  4-12(DEG) 


G.6 


G.7  F.FIELD  AVERAGE  SLL  4-12(DEG) 


AN  INCIDENT  FIELD  SENSOR  FOR  EMP  MEASUREMENTS 


Everett  G.  Farr  and  Joseph  S.  Hofstra 
The  BDM  Corporation 


Abstract 


When  making  field  measurements  at  EMP  simulators,  one  needs  broadband  anten¬ 
nas  that  are  flat  to  above  100  MHz.  In  additional,  it  would  be  helpful  if  these  antennas 
were  directional,  in  order  to  reject  reflections  from  directly  behind  the  sensor.  This  becomes 
important  in  extrapolation  procedures  that  correct  for  differences  between  the  field  gener¬ 
ated  by  the  simulator  and  an  actual  threat  field.  Currently  available  antennas,  however, 
are  short  dipoles  and  small  loops,  which  are  omnidirectional  in  one  plane.  The  purpose  of 
this  paper  is  to  report  on  the  development  Gf  an  antenna  that  solves  this  problem  with  a 
1  +  cos#  antenna  pattern  in  both  the  E-  and  H-planes.  We  call  this  antenna  the  Balanced 
Transmission-line  Wave  (BTW)  sensor. 

The  design  of  the  sensor  is  based  on  the  idea  that  if  one  could  sum  the  antenna 
patterns  of  a  short  dipole  and  a  small  loop,  one  could  achieve  a  directional  sensor.  Although 
a  technique  for  combining  these  patterns  in  a  single  sensor  has  been  known  for  some  time, 
to  our  knowledge,  a  working  model  has  never  been  built  and  tested. 

In  this  presentation,  a  simple  theory  of  operation  is  discussed.  Furthermore,  ex¬ 
perimental  results  are  presented  for  a  prototype  sensor.  A  Front /Back  ratio  somewhere 
between  20  and  29  dB  was  achieved.  These  results  confirm  the  principal  of  operation  of  the 
sensor,  and  encourage  further  engineering  development. 
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1.  INTRODUCTION 


This  paper  describes  the  development  of  a  novel  antenna  that  is  both  electrically 
small  and  directional.  It  was  developed  to  meet  the  need  for  a  directional  antenna  when 
making  EMP  field  measurements  whose  response  would  be  flat  to  100  MHz  [1].  The  antennas 
that  are  typically  used  in  EMP  field  measurements  are  small  dipoles  and  loops,  but  these  are 
not  directional.  The  antenna  we  have  developed  is  called  the  Balanced  Transmission-line 
Wave  (BTW)  sensor. 

The  design  of  the  sensor  is  based  on  an  idea  developed  for  a  unidirectional  EMP 
transmitter  to  be  used  in  a  simulator  [2].  Although  the  original  device  was  proposed  some 
time  ago,  it  has  not  to  our  knowledge  been  built  or  tested. 

In  this  note  a  simplified  theory  of  operation  is  discussed.  Furthermore,  experimental 
results  are  presented  for  two  prototype  sensors  in  a  number  of  different  types  of  tests.  These 
results  confirm  the  principle  of  operation  of  the  sensor  and  encourage  further  engineering 
■development  of  the  BTW. 


.2.  APPROACH 

A  simple  view  of  a  BTW  is  shown  in  Figure  1(a).  It  consists  of  a  100-Q  transmission 
line  that  tapers  to  a  100-12  load  on  either  end.  This  load  can  be  either  a  100-12  resistor  or  a 
100-12  twinax  cable.  For  structural  reasons,’  an  extra  ground  plane  is  normally  inserted  at 
the  symmetry  plane.  In  doing  so,  the  100-11  resistors  become  two  50-12  resistors,  and  the 
100-12  twinax  terminations  become  two  50-12  coax  transmission  lines. 

The  structure  of  Figure  1(a)  has  been  analyzed  in  detail  by  Jiunn  Yu  et  al  in  Sensor 
and  Simulation  Note  (SSN)  243.  The  treatment  in  SSN  243,  however,  is  more  complicated 
than  what  is  required  for  a  simple  understanding  of  the  BTW.  Therefore  we  begin  with  a 
simplified  low-frequency  explanation  of  how  and  why  the  device  operates. 

The  BTW  can  be  described,  most  simply  in  terms  of  a  parallel  plate  transmission 
line.  This  is  later  generalized  to  an  arbitrary  transmission-line  configuration.  We  begin, 
with  the  parallel  plate  model. 


2.1  Parallel  Plate  Transmission  Line  Model 

Consider  the  diagram  of  Figure  1(b).  It  consists  simply  of  two  closely  spaced  parallel 
plates  that  form  a  transmission  line.  Each  end  of  the  line  is  assumed  to  be  terminated  in 
the  characteristic  impedance  of  the  line.  We  can  describe  the  characteristics  of  the  device 
in  terms  of  a  simple  equivalent  circuit,  shown  in  Figure  1(c).  This  is  a  current  and  voltage 
source  feeding  into  two  matched  loads.  The  voltage  and  current  sources  are  derived  from 
inductive  and  capacitive  coupling,  respectively. 
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First,  we  derive  the  voltage  source  due  to  inductive  coupling.  If  a  plane  wave  is 
incident  upon  the  sensor  as  shown  in  Figure  1(b),  then  the  induced  voltage  through  the 
loop  in  the  y-z  plane  is 

m=-A-  f  m 

where  /I  =  Ihx.  Since  B  =  fi0Hxx,  we  can  express  the  above  as 

V(t)  =  -fi0lh (2) 
where  fi0  =  4;r  x  10-7H/m.  This  behavior  is  directly  analogous  to  a  B-dot  sensor. 


Next,  we  derive  the  current  source  due  to  capacitive  coupling.  For  the  configuration 
of  Figure  1(b),  the  current  induced  between  the  plates  is 


l(t)  =  C 


dV 

dt 


(3) 


If  we  assume  a  uniform  electric  field  between  the  plates,  then  the  voltage  between  the  plates 
is 

V(t)  =  -hy-E(t)  (4) 

where  E(t)  is  the  electric  field  due  to  the  incident  field.  If  we  also  assume  that  fringing  is 
negligible  at  the  edge  of  the  plates,  then 


C  = 


€plw 

h 


(5) 


where  e0  =  8.85  x  10-12F/m.  Combining  thv  above  three  equations,  we  find  the  current 
source  as 

dF 

At)  =  -e0  Iw  (6) 

It  turns  out  that  this  current  source  can  be  simply  expressed  in  terms  of  the  voltage  source, 
V(<).  This  is  accomplished  by  assuming  Ev/Hx  =  r/  =  where  77  =  377  0.  is  the 

impedance  of  free  space.  By  doing  this,  we  have 


l(t)  =  -e07]lw 


dHr 

dt 


/i0  lh 


dHx 

dt 


(7) 


Since  the  characteristic  impedance  of  a  parallel  plate  transmission  line  is  Z0  =  r\h/w  we 
have,  by  combining  Equations  (2)  and  (7) 


J(t)  = 


m 


(8) 


This  is  the  desired  result. 
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It  is  dear  from  Figure  1(c)  that  if  Equation  (8)  is  true,  then  the  voltage  at  Port  1 
of  the  sensor  is 

Vi  (t)  =  V(t)  (9) 


whereas  the  voltage  at  Port  2  is 


*4(0  =  0 


(10) 


This  is  the  behavior  required  for  a  directional  sensor.  If  the  wave  is  incident  from  the 
-f -direction,  it  is  simple  to  show  with  symmetry  that  Port  2  is  excited  and  Port  1  has  no 
response. 


2*2  Arbitrary  Transmission  Line  Model 


It  is  not  necessary  to  have  a  parallel  plate  transmission  line  to  observe  the  effect 
described  in  the  previous  section.  We  now  demonstrate  that  the  phenomenon  remains  valid 
with  an  arbitrary  transmission  line. 


The  critical  result  that  needs  to  be  proven  is  Equation  (8).  We  need  a  more  general 
expression  for  X(£).  In  general, 


1(0  =  Cil 


mt) 

dt 


(ii) 


where  Ci  is  tho  capacitance  per  unit  length  of  the  transmission  line,  l  is  the  length  of  the 
line,  and  V(t)  is  a:  defined  in  Equation  (4).  Using  Equation  (4),  we  get 


I(t)  =  -hCily 


dm  hirdEy  hirdHx 


(12) 


Furthermore,  V(t)  is  unchanged  from  Equation  (2),  since  the  area  of  the  loop  does  not 
change.  Thus,  combining  Equations  (2)  and  (12),  we  find 


V(£)  _  "o  _  \/MofO  _  2 
m  "  Vci  ci  ~  0 


(13) 


The  last  equation  follows  from  a  definition  of  characteristic  impedance  tha.  holds  for  trans¬ 
mission  lines  with  air  dielectric.  This  confirms  that  Equation  (8)  is  valid  for  a  transmission 
line  of  arbitrary  cross  section  provided  that  its  phase  velocity  equal  to  the  speed  of  light. 
This  allows  us  flexibility  in  the  final  design  of  the  sensor. 


2.3  Sensor  Bandwidth 

The  frequency  bandwidth  is  of  concern  when  attempting  to  use  the  BTW  at  higher 
frequencies.  It  is,  however,  a  somewhat  difficult  number  to  define.  For  our  purposes,  a 
reasonable  definition  would  be  the  frequency  ban^  over  which  a  minimum  Front/Back  ratio 
is  maintained.  This  turns  out  to  be  somewhat  difficult  to  measure,  since  our  measurements 
of  Front/Back  ratio  are  taken  in  pulse  mode  and  not  in  CW  mode. 
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With  this  in  mind,  our  best  definition  of  bandwidth  is  probably  given  by  Carl  Baum 
in  Sensor  and  Simulation  Note  8  [3].  This  was  developed  originally  for  a  Multi-Gap  Loop 
(MGL)  sensor,  and  since  there  is  an  analogy  between  the  performance  of  a  BTW  and  MGL, 
this  definition  may  be  reasonable.  According  to  SSN  8,  the  maximum  frequency  of  operation 

iS  1 

fmax  =  ^  (14) 

where  %  is  the  round  trip  transit  time  of  a  signal  on  the  sensor.  If  there  is  no  dielectric  in 
the  transmission  line,  then  %  =  2//c,  and 


/max  -  2jfi 


(15) 


As  an  example,  the  BTW  we  designed  could  be  built  to  have  a  length  about  0.35  m  if  some 
extra  lengths  of  transmission  line  were  removed.  For  this  length,  the  bandwidth  would  be 
about  135  MHz.  Of  course,  the  bandwidth  could  always  be  increased  by  decreasing  the  size 
of  the  BTW,  which  in  turn  would  decrease  the  effective  area  of  the  BTW. 


2.4  Design  Considerations 

A  diagram  of  the  BTW  is  shown  in  Figure  2.  The  effective  area  was  designed  to  be 
close  to  that  of  an  MGL-2,  0.01  m2.  This  is  the  same  as  the  area  of  the  loop  one  sees  when 
looking  at  the  side  view  of  the  sensor.  The  ground  plane  acts  as  a  structural  support,  and 
does  not  affect  the  fields  due  to  symmetry.  The  characteristic  impedance  of  each  side  of  the 
parallel  plate  transmission  line  is  50-fi,  and  it  was  calculated  using  reference  [4]. 

The  output  of  the  BTW  is  two  SMA  connectors  at  Port  1,  which  are  each  connected 
by  cable  to  one  of  the  inputs  of  a  differential  mode  fiber  optic  transmitter.  In  fact,  it  is  pos¬ 
sible  to  measure  fields  from  both  directions  concurrently.  This  can  be  done  by  replacing  the 
two  50-11  loads  at  Port  2  with  cables  that  are  connected  to  a  second  fiber  optic  transmitter. 

An  incident  field  excites  the  end  of  the  transmission  line  structure  port  closest  to 
the  direction  of  incidence.  The  cables  from  Port  1  are  routed  to  the  back  end  of  the  sensor 
in  order  to  keep  instrumentation  behind  the  incident  field. 


3.  LABORATORY  MEASUREMENTS 

Two  experiments  were  performed  in  the  laboratory  to  check  the  performance  of  the 
BTW.  First,  a  Time  Domain  Reflectometry  (TDR)  measurement  was  performed  to  check 
the  impedance  matching  between  different  sections  of  the  sensor.  We  connected  the  TDR 
input  to  one  side  of  Port  1,  leaving  the  remaining  connector  of  Port  1  and  both  connectors 
of  Port  2  matched  to  50-D  loads. 

Results  for  the  TDR  are  shown  in  Figure  3.  One  can  relate  bumps  in  the  TDR 
measurements  to  physical  features  of  the  structure.  The  maximum  deviation  from  the  zero 
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line  is  about  50  mp,  or  about  5  percent  of  the  maximum  possible  reflection  coefficient,  which 
is  1. 


Next,  TEM  cell  measurements  were  performed.  In  derivative  sensors,  one  expects  to 
see  the  measured  response  increase  proportionally  with  frequency,  or  20  dB/decade.  TEM 
cell  measurements  were  made  in  order  to  confirm  this. 

A  diagram  of  the  experimental  setup  is  shown  in  Figure  4.  The  BTW  was  placed 
in  a  3-ft  TEM  cell  and  the  output  was  fed  through  a  balun  into  a  network  analyzer.  The 
output  of  the  sensor  was  referenced  to  the  total  throughput  of  the  TEM  cell.  The  BTW 
was  placed  in  a  manner  that  would  maximize  the  coupling  into  it. 

Results  are  shown  in  Figure  5.  The  response  increases  at  20  dB/decade,  exactly  as  it 
should.  There  are  a  few  small  bumps  in  the  data  in  the  decade  between  10  and  100  MHz,  but 
these  are  limitations  of  the  TEM  cell  technique  that  normally  occur  at  higher  frequencies. 
This  is  simply  demonstrated  by  replacing  the  cables  with  ones  of  a  slightly  different  length 
and  watching  the  bumps  move  to  a  different  frequency. 


4.  EXPERIMENTS  AT  HPD 

A  number  of  tests  were  performed  at  HPD  to  verify  the  performance  of  the  BTW. 
HPD  is  an  EMP  test  facilty  and  it  is  shown  in  Figure  *3.  It  consists  of  a  pulser  30  m  above 
the  ground  that  is  connected  to  an  elliptical  loop  antenna  that  connects  to  the  ground  at 
either  end.  This  was  the  facility  that  spurred  the  development  of  the  BTW  sensor;  we 
wanted  to  measure  the  incident  field  by  itself,  while  rejecting  the  reflected  field. 

Two  tests  were  performed  at  HPD.  First,  measurements  were  made  comparing  a 
BTW  to  an  MGL  when  used  under  the  pulser  in  a  normal  mode  of  operation.  Since  the 
MGL  is  a  loop  sensor,  one  would  expect  to  see  both  incident  and  reflected  fields  with  the 
MGL,  and  only  the  incident  field  with  the  BTW.  Next,  the  BTW  was  pointed  down  in 
order  to  observe  how  well  it  rejected  the  incident  field. 

In  the  first  experiment,  the  response  of  the  BTW  was  compared  to  that  of  an  MGL- 
2  on  the  same  shot,  while  both  were  positioned  5  m  above  the  ground,  directly  under  the 
HPD  pulser.  This  is  the  normal  mode  of  operation  for  the  BTW.  We  expected  to  see  the 
BTW  reject  the  reflected  field  while  the  MGL  would  measure  it  The  BTW  was  directly 
under  the  pulser,  while  the  MGL  was  2  m  outside  the  plane  of  the  loop  of  the  HPD  antenna. 
Since  the  HPD  pulser  is  30  m  above  the  ground,  this  is  effectively  directly  under  the  pulser. 

For  this  experiment,  the  output  of  the  BTW  and  MGL-2  were  each  fed  directly 
into  a  differential  mode  Nanofast  fiber  optic  link  with  1  jis  integrators  turned  on  inside 
the  transmitters.  Since  the  output  of  the  MGL  is  a  twinax  connector,  a  twinax-to-coax  tee 
(TCT)  was  used  to  convert  the  twinax  to  two  coaxial  lines.  The  receivers  fed  into  LeCroy 
6880  A/D  converters.  Noise  filters  were  turned  on  in  the  LeCroy  digitizers  at  150  MHz. 
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These  data  were  processed  using  standard  signal  processing  techniques.  The  zero 
line  was  found  by  averaging  the  first  two  hundred  points.  The  effect  of  the  fiber  optic  links 
was  deconvolved  from  the  output  using  standard  Fourier  transform  techniques.  Although 
we  did  not  know  exactly  what  the  effective  area  of  the  BTW  would  be,  we  had  to  assume 
an  effective  area  in  order  to  get  the  data  through  the  standard  processing  algorithms.  For 
this  reason,  we  assumed  during  our  data  processing  that  the  effective  area  of  the  BTW  was 
0.01  m2,  as  it  was  designed.  Later,  during  testing  at  ALECS,  we  conducted  an  experiment 
to  measure  this  quantity  (see  Section  5.1). 

The  waveforms  measured  with  the  MGL-2  and  the  BTW  are  shown  in  Figure  6. 
In  comparing  the  two  results,  it  is  clear  that  the  second  large  bump  in  the  response  of  the 
MGL  is  the  reflected  field.  It  occurs  at  about  33  ns  after  the  first  peak,  which  is  the  two-way 
transit  time  between  the  sensor  and  the  ground.  The  response  of  the  BTW,  however,  is 
missing  this  second  bump  at  33  ns. 

Next,  we  placed  the  BTW  pointing  down,  5  m  above  the  ground,  in  order  to  observe 
the  rejection  of  the  incident  field.  The  same  instrumentation  and  signal  processing  was  used 
as  before,  except  that  a  150  MHz  filter  was  introduced  in  data  processing  rather  than  in 
data  acquisition. 

A  plot  of  the  measured  field  is  shown  in  Figure  7.  Ideally,  we  expect  to  see  the  first 
33  ns  of  the  signal  to  be  noise,  with  a  large  signal  following  that  being  the  reflected  field. 
Our  measurment  shows  essentially  this;  the  first  33  ns  of  the  signal  is  almost  in  the  noise. 


5.  EXPERIMENTS  AT  ALECS 

The  final  set  of  tests  took  place  at  ALECS.  This  facility  was  used  because  it  is  a 
transmission-line  type  simulator.  It  is  shown  in  Figure  9.  It  provides  a  TEM  field  with 
almost  no  reflections.  First,  we  attempted  to  calibrate  the  effective  area  of  the  BTW 
by  comparing  its  response  to  that  of  an  MGL.  This  was  done  with  both  integrated  and 
derivative  data.  Next,  we  measured  the  Front/Back  ratio  of  the  BTW,  again  using  both 
integrated  and  derivative  data.  Finally,  we  made  measurements  of  the  antenna  pattern  of 
the  BTW  in  the  E-  and  H-planes. 


5.1  Effective  Area  Calibration 

The  first  test  was  to  compare  the  effective  area  of  the  BTW  tapered  design  to  that 
of  nn  MGL-2.  To  do  this,  we  compare  the  ratio  of  peak  magnitudes  of  the  MGL-2  and  the 
BTW  on  the  same  shot.  This  technique  requires  the  assumption  that  the  effective  area  of 
the  MGL-2  is  0.01m2,  and  its  accuracy  is  limited  by  this  assumption. 

An  MGL-2  and  the  BTW  were  placed  2  m  apart  in  the  center  of  the  testing  volume. 
These  sensors  were  each  fed  into  a  differential  mode  fiber  optic  link  with  the  1  fxs  integrators 
turned  on.  The  connection  to  the  link  for  the  MGL  required  a  TCT  to  split  the  twinax 


output  into  into  two  coaxial  lines.  The  output  of  the  fiber  optic  receivers  was  then  fed  into 
an  active  signal  divider,  which  splits  the  signal  into  three  parts  and  sends  it  on  to  three 
Tektronix  7912  digitizers. 

Signal  processing  involved  the  normal  process  of  finding  the  zero  line  by  averaging 
the  first  50  points  and  then  deconvolving  the  effect  of  the  fiber  optic  transmitter.  As  before, 
we  assumed  an  effective  area  of  0.01  m2  for  the  BTW  to  get  the  data  through  processing. 
We  can  calculate  the  correction  needed  in  the  effective  area  if  we  start  with  an  assumed 
value  for  it. 

The  peak  magnitude  as  measured  with  \ie  BTW  was  compared  to  that  measured 
with  the  MGL-2  for  four  shots.  The  average  normalized  peak  magnitude  of  the  field  mea¬ 
sured  with  the  BTW  was  1.1.  This  means  that  the  effective  area  of  the  BTW  is  10  percent 
larger  than  that  of  the  MGL-2,  or  0.011  m2.  A  summary  of  this  data  appears  in  Table  1. 

Table  1.  Effective  Area  and  Front/Ba^  Ratio  at  ALECS 

Data  T7Pe  active  Area  .*ra2>  rraat/Iv:  &  'dB^ 

Integra  tad  Q.QUO  IQ 

DarL7ati7«  0,  CHS  IS 

This  process  was  then  repeated  with  the  integrators  turned  off.  The  signal  processing 
involved  only  time  tying,  finding  the  baseline,  and  then  linear  scaling  by  a  factor  of  1  /(/i0Ae). 
Again,  a  normalized  peak  was  calculated  by  comparing  the  peak  of  the  BTW  data  to  that 
of  the  MGL-2  data.  The  average  normalized  peak  magnitude  of  the  field  measured  with  the 
BTW  was  1.15.  This  means  that  the  effective  area  of  the  BTW  is  15  percent  larger  than 
that  of  the  MGL-2,  or  0.0115  m2. 

There  is  some  question  as  to  which  technique  is  more  accurate,  since  we  get  a 
5  percent  difference  in  effective  areas  with  the  two  techniques.  This  is  especially  true  since 
the  data  in  general  do  not  have  a  wide  scatter  about  their  means.  Normally,  one  would 
prefer  to  use  the  data  from  the  measurment  technique  with  the  fewest  signal  conditioners 
to  be  deconvolved  out.  This  suggests  that  the  derivative  data  are  more  accurate,  since  an 
integrator  is  not  used.  This  does  not  explain,  however,  why  there  should  be  a  discrepancy 
between  the  integrated  and  derivative  data.  Perhaps  a  reasonable  compromise  would  be  to 
take  the  average  of  the  two  techniques,  which  is  0.012.  This  is  the  effective  area  we  used  in 
the  Front/Back  ratio  measurements. 

5.2  Front/Back  Ratio  Measurements 

The  next  experiment  involved  a  measurement  of  the  Front /Back  ratio.  This  was 
done  with  both  derivative  and  integrated  measurements,  using  identical  test  setups  and 
signal  processing  as  in  the  HPD  experiment.  The  only  difference  was  that  the  BTW  was 
pointed  away  from  the  pulser  instead  of  toward  it.  By  using  the  peak  magnitudes  and 
correcting  for  the  difference  in  effective  areas  between  the  MGL  and  BTW,  one  can  calculate 
a  Front/Back  ratio. 

The  peak  magnitude  of  a  waveform  with  the  BTW  pointing  backward  was  compared 
to  the  peak  magnitude  for  an  MGL  over  several  shots.  After  taking  into  account  the  different 
effective  areas  between  the  MGL  and  BTW,  there  was  an  average  20  dB  Front/Back  ratio 
with  integrated  data. 


3 


534 


The  same  experiment  was  repeated  with  derivative  data.  This  time  we  found  a 
Front/Back  ratio  of  29  dB.  As  in  the  previous  experiment,  we  obtained  different  results 
for  the  integrated  and  derivative  data.  This  time,  however,  the  difference  in  results  was 
more  significant,  9  dB.  One  way  to  resolve  the  issue  might  be  to  assume  that  the  simplest 
experiment,  i.e.,  the  one  with  derivative  data,  is  the  most  accurate.  It  is  a  little  difficult  to 
say  this  definitively,  however,  and  it  will  have  to  be  left  for  future  investigation. 


5.3  Anterna  Pattern  Cuts 

In  the  final  experiment,  we  measured  the  antenna  pattern  of  the  BTW.  For  far  field 
measurements,  the  BTW  should  have  an  antenna  pattern  that  is  1  +  cos  9,  according  to 
SSN  243  [2].  This  behavior  should  hold  both  in  the  E-  and  H-planes. 

The  data  sets  taken  were  all  integrated,  with  the  instrumentation  and  signal  pro¬ 
cessing  the  same  as  that  for  all  previous  integrated  data  taken  at  ALECS.  Plots  of  the  cuts 
are  given  in  Figures  8  and  9.  It  is  clear  that  the  sensor  performs  about  as  we  expect,  with 
the  exception  near  180  degrees.  We  should  remember,  however,  that  the  Front/Back  ratio 
with  derivative  data  was  9  dB  better  than  that  for  integrated  data,  so  if  we  had  made  these 
plots  with  derivative  data  we  would  expect  better  agreement  at  180  degrees.  Finally,  we 
note  that  each  of  the  points  on  these  graphs  was  generated  with  a  single  shot,  not  an  aver¬ 
age  over  several  shots.  Therefore,  one  might  expect  a  few  points  would  be  slightly  different 
from  theory. 


0.  CONCLUSIONS 

A  broadband  directional  antenna  has  been  developed  that  should  simplify  and  im¬ 
prove  the  measurent  of  the  incident  field  at  hybrid  EMP  simulators.  This  antenna,  the 
Balanced  Transmission-line  Wave  sensor,  has  the  property  of  being  able  to  measure  the  in¬ 
cident  field 'while  rejecting  the  reflected  field.  Preliminary  measurements  of  the  Front/Back 
ratio  suggest  an  isolation  of  29  dB  with  derivative -data  and  20  dB  with  integrated  data. 

There  remain  two  areas  where  further  investigation  would  be  useful.  The  first  is 
whether  or  not  the  Front/Back  ratio  could  be  improved  with  further  experimentation.  There 
are  a  number  of  design  parameters  that  could  be  varied  while  maintaining  the  same  overall 
response  of  the  sensor.  A  systematic  study  of  these  variables  may  reveal  designs  with  a  larger 
Front/Back  ratio.  The  second  area  where  further  investigation  would  be  useful  concerns 
the  discrepancy  between  integrated  and  derivative  measurements  of'Front/Back  ratio.  At 
some  point,  it  will  be  necessary  to  generate  a  commonly  accepted  method  of  measuring  this 
parameter. 

An  interesting  aspect  of  this  work  is  that  it  may  have  applications  outside  the  area 
of  EMP.  It  would  seem  there  could  be  a  number  of  uses  for  a  broadband  antenna  that  is 
both  electrically  small  and  directional. 
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Simplified  BTW  Diagram  (Cone 
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Instrumentation  Setup  for  TEM  Cell  Measurements 
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Figure  7.  Fields  Measured  Under  the  1IPI)  Pulser  With  an  MGI.-2  and  the  BTW 


Figure  H.  Rejection  of  the  Incident  Field  at  HPI)  With  the  BTW 
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Figure  10.  11-Plane  Cut  of  the  BTW  Measured  at  ALECS 
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Figure  11.  E-Plane  Cut  of  the  BTW  Measured  at  A 
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ABSTRACT 

Quantifying  Simultaneous  V  and  H  Coverage  from  an  Underside  Aircraft 

Location 

By:  George  Monser  and  Leonard  T.  Surette 

This  paper  presents  the  results  of  a  measurement  program  conducted  to  evaluate 
antenna  coverage  from  a  constrained-space  airborne  platform. 

Two  levels  of  coverage  are  reported:  Level  l  details  coverage  without  the  aircraft, 
but  over  a  small  groundplane,  while  Level  2  includes  aircraft  effects.  Coverage 
and  gain  for  each  polarization  was  evaluated  by  quantizing  the  data  and  computer 
processing  to  obtain  a  measure  of  goodness  vs  frequency.  In  particular,  horizontal 
polarization  showed  considerable  changes,  as  expected,  compared  to  Level  1 
performance.  Vertical  polarization  showed  smaller,  but  significant  changes. 
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Coverage,  Under-Aircraft,  Simultaneous 
By:  George  J.  Monser  and  Leonard  T.  Surette 


SUMMARY: 

This  paper  presents  the  results  of  a  measurement  program  conducted  to  evaluate 
antenna  coverage  from  a  constrained-space  airborne  platform. 

Two  levels  of  coverage  are  reported:  Level  1  details  coverage  without  the 
aircraft,  but  over  a  small  ground-plane,  while  Level  2  includes  aircraft  effects. 
Coverage  and  gain  for  each  polarization  was  evaluated  by  quantizing  the  data  and 
computer  processing  to  obtain  a  measure  of  goodness  vs  frequency.  In  particular, 
horizontal  polarization  showed  considerable  changes,  as  expected,  compared  to 
Level  1  performance.  Vertical  polarization  showed  smaller,  but  significant 
changes. 

1.0  INTRODUCTION 


This  paper  presents  the  results  of  a  measurement  program  conducted  to  evaluate 
coverage  from  a  under-aircraft  platform.  In  particular,  the  antenna  installation 
was  small,  in  terms  of  wavelengths,  the  bandwidth  was  greater  than  3/1, 


smiuuaneous,  < 


azimuth  by  300  in  elevation. 
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2.0  INSTALLATION 


The  complete  antenna  system  was  constrained  to  fit  inside  the  volume 
dimensioned  «n  wavelengths  shown  in  Figure  1 .  In  addition,  the  bottom  surface, 
which  was  metal,  would  be  positioned  under-side  the  aircraft. 

3.0  DESIGN  CONSIDERATIONS 

The  first  consideration  is  how  to  attain  simultaneous,  field-of-view,  dual¬ 
polarization  from  the  volume  shown  in  Figure  1.  For  vertical  polarization  (i.o., 
perpendicular  to  the  aircraft  surface  underside)  the  best  candidate  appears  to  be 
some  fonn  of  conical  monopole.  This  type  of  antenna  can  be  designed  to 
efficiently  cover  3/1  bandwidth  and  the  radiation  field-of-view  approximates  the 
requirement.  In  addition,  the  radiating  structure  is  1/2-height  for  the  ground- 
plane  image  forms  the  other  half.  This  is  an  important  consideration  approaching 
Flow  where  the  allowable  height  (see  Figure  1)  is  only  0.13  wavelengths. 

For  horizontal  polarization  (i.e.,  parallel  to  the  aircraft  surface)  some  form  of 
horizontal  turnstile  fan-dipoles  can  be  used.  Here  it  is  known  that  the  effect  of  the 
aircraft  surface  will  result  in  tilting  the  radiation  downward  with  the  maximum 
radiation  occurring  approximately  at  a  depression  angle  below  the  aircraft  of 
A/4H  where  X  is  the  wavelength  and  H  is  the  antenna  phase-center  spacing  with 
respect  to  the  aircraft  surface.  From  these  considerations  and  the  space 
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Figure  1,  the  arrangement  shown  in  Figure  2  was  devised  with  the  fan  dipoles 
placed  flat  in  the  upper  surface  of  the  installation  volume. 

Installing  both  the  vertical  and  horizontal  antennas  in  a  common  volume  leads  to 
antenna  shadowing  and  cross-talk  effects  both  of  which  were  considered 
unacceptable. 

One  solution  is  to  integrate  the  vertical  and  horizontal  radiating  structure  into  a 
common  single-input  structure  as  shown  in  Figure  3.  Here,  all  elements  of  the 
vertical  monopole  are  retained  and  an  asymmetrical  top-hat  is  added  to  the 
monopole  yielding  some  degree  of  horizontal  polarization.  Thus,  omnidirectional 
coverage  is  maintained  for  vertical  polarization  and  a  quasi-omni  horizontal 
polarization  results. 

4.0  PRELIMINARY  TESTS 


Before  installing  on  the  aircraft  model  an  antenna  was  configured  to  fit  the 
volume  constraints  of  Figure  1.  Figure  4  shows  a  photograph  of  one  of  many 
antenna  assemblies  which  were  tested.  Here,  the  vertical  monopole  was 
positioned  in  the  center  of  the  metal  base  plate.  To  approximate  a  conical 
monopole  two  triangular  plates  were  used,  one  with  a  flare  angle  of  90°  and  the 
second  with  a  flare  angle  of  65°.  Next  unbalanced  top-loading  was  added  to 
achieve  horizontal  polarization. 
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Tests  with  a  balanced  and  unbalanced  top-hat  showed  that  in  the  un-balanced 
configuration  that  significant  horizontal  polarization  could  be  provided  without 
affecting  the  vertical  omni  coverage. 

5.0  MODEL  TESTS 

Figure  5  shows  the  1/10  scale-model  aircraft  with  the  antenna  system  installed  on 
the  underside  of  the  aircraft. 

Radiation  patterns  and  gain  were  obtained  for  a  number  of  test  frequencies  for 
four  conic  cuts  and  two  polarizations.  Specifically,  the  elevation  cuts  were  0,  -10, 
-20,  and  -30°  in  elevation  over  360°  in  azimuth. 

6.0  TYPICAL  INSTALLED  PATTEPNS 

Figure  6  shows  one  of  many  patterns  recorded  during  the  tests.  In  Figure  6(A)  the 
two  principal  polarizations  at  FLOW  are  plotted  for  0o  in  elevation.  While  the 
horizontal  polarization  is  not  particularly  good,  it  improves  as  shown  in  Figure 
6(B)  as  the  depression  angle  approaches  -30°.  In  Figure  6(B)  two  elevation  cuts 
are  displayed,  one  taken  Fore  and  Aft  and  the  other  Broadside. 
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7.0  INSTALLED  PERFORMANCE 


Data  from  plots  similar  to  Figure  6  were  used  to  compile  gain  vs  frequency 
profiles  for  each  polarization  as  shown  in  Figure  7.  Each  point  in  Figure  7 
represents  an  average  for  360°  in  azimuth  sampled  at  10°  intervals  over  four 
conic  cuts.of  0,  -10,  -20,  and  -30°. 

Figure  8  shows  how  the  average  gain  over  360°  in  azimuth  varies  with  elevation 
for  each  polarization  for  three  frequencies. 

8.0  DISCUSSION 

Data  in  the  preceding  section  illustrates  what  can  be  achieved  over  a  wide 
bandwidth  for  two  polarizations  from  a  constrained  space. 

At  first  consideration  the  gain  values  given  in  Figures  7  and  8  may  appear  low, 
particularly  in  Figure  7.  The  following  assessment  is  offered: 

1.  Power  division  between  the  two  polarizations;  ideally  3  dB. 

2.  VSWR  loss  1  to  3  dB,  based  upon  measured  data. 

3.  Radiated  energy  outside  the  defined  field-of-view  of  0  to  -30®  in  elevation. 
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From  the  above  tally  only  an  improvement  of  gain  by  improving  VSWR  appears 
relevant.  At  most,  about  2  dB  is  recoverable  so  that  the  design  for  vertical 
polarization  can  be  made  to  approach  isotropic  gain.  For  horizontal  polarization 
the  average  gain  can  be  improved  to  an  average  value  of  about  -3  to  -4  dB.  As  an 
alternative  the  design  can  be  re-configured  to  enhance  horizontal  gain  (i.e., 
unequal  power  split)  which  should  lead  to  about  -  2  dB  for  each  polarization. 

9.0  CONCLUDING  REMARKS 

The  paper  has  shown,  by  measurements,  what  level  of  performance  can  be 
expected  over  a  broad  bandwidth  from  a  small-volume-installation  positioned 
under  an  aircraft . 

When  the  installation  space  is  constrained  to  approximately  0.1  x  0.1  x  0.5 
wavelengths  balanced  average  gains  on  the  order  of  0  to  -3  dB  should  be 
realizable  over  360°  in  azimuth  by  0  to  -30°  in  elevation. 
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Installation  Space  in  Wavelengths  at  FlOW 
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Figure  2.  Antenna  Concepts 


Figure  3.  Conceptual  Details  of  Integrated  Antenna 


Figure  4.  Test  Model 
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Figures.  1/10  Scale  Aircraft  (Antenna  Installed) 


Azimuth  Plot 
<g>  0°  Elevation 
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Figure  6(A).  A/imuth  Pattern,  0°  Elevation 
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Figure  7.  Average  Installed  Cain  versus  Installation  Height 


Average  Installed  Gain 
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e  8.  Average  Installed  Cain  versus  Elevation  Depression  Angle 


APPENDIX  A 


LEVEL  1  -  UNINSTALLED  PERFORMANCE 


This  appendix  presents  a  brief  summary  of  tests  for  the  dual-polarized  antenna 
before  installing  on  the  aircraft. 

Figure  A1  shows  return  loss  for  the  dual-polarized  antenna.  Here,  the  return  loss 
is  typically  below  -7  dB  except  at  the  lower  frequencies  where  the  loss 
approaches  -4  dB  at  FlOW-  Thus  over  most  of  the  band  a  transmission 
efficiency  of  better  than  80%  results. 

Figure  A2  shows  typical  the  azimuth  coverage  for  tvi  o  polarizations.  Here,  the 
test  frequency  corresponds  to  upper  mid-band  and  the  gains  for  each  polarization 
are  referenced  to  a  gain-standard  level  as  indicated. 

Figure  A3  shows  azimuth  coverage  for  two  polarizations  for  a  -30°  comic 
taken  parallel  to  the  azimuthal  plane.  Somewhat  higher  response,  as  expected,  is 
evident  for  horizontal  polarization  than  shown  in  Figure  A2. 

Figure  A4  shows  a  compendium  of  gain  plots  taken  to  establish  expected  gain  for 
this  type  of  antenna.  All  of  these  models  v/ere  symmetrical  in  construction  and, 
as  such,  represent  the  achievable  gains  for  vertical  polarization.  Using  the  data 
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from  the  VSWR  plot  and  estimating  the  radiation  efficiency  for  the  antennas  the 
following  gain  tally  results: 

1.  For  one  antenna  tested,  height  equal  to  1/12  wavelength  the 
measured  gain  was  -2.9  dB. 

2.  The  theoreticai  gain  was  about  +3dB. 

3.  The  difference  is  about  6  dB. 

4.  The  VSWR  transmission  loss  was  approximately  3  dB. 

5.  The  radiation  efficiency  was  computed  as  50%  or  -3  dB. 

6.  Thus  the  sum  of  the  VSWR  and  radiation  efficiency  yield  6  dB 
which  is  in  good  agreement  with  the  difference  noted  in  step  3. 


i 
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Figure  A2.  Azimuth  Patterns,  Two  Polarizations  (0°  Elevation) 
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Patterns,  Two  Polarizations  (-30°  Conic) 
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POLARIZATION  DETERMINATION  IN  A  DF  SYSTEM 
AT  MM  WAVELENGTHS  WITHOUT 
A  ROTARY  JOINT  * 

Harry  B.  Sef+on  Jr, 


ABSTRACT 

With  wider  use  of  the  mm  bands  18-26  and  26-40  GHz,  it  is  only  natural  that 
the  surveillance  community  would  wish  to  detect,  locate,  and  characterize  such 
emitters.  The  characteristic  which  this  paper  addresses  is  the  determination 
of  the  emitted  polarization  as  part  of  the  Direction  Finding  (DF)  process.  In 
response  to  this  challenge,  a  technique  of  polarization  identification  in  the 
mm  wave  bands  was  devised  by  the  writer,  and  implemented  by  TECOM  as  part  of 
production  of  1-40  GHz  spinning  direction  finding  systems.  The  technique  is 
based  on  the  use  of  orthogonal  grid  sets  which  rotate  synchronously  with  a  45° 
inclined  electromagnetic  mirror  in  such  a  manner  as  to  act  as  polarization 
filters  for  pairs  of  receiving  horns  positioned  behind  the  grids.  System  test 
results  of  actual  hardware,  taken  inside  an  A-sandwich  radome,  show 
identification  of  the  correct  emitted  polarization  with  a  margin  of  12-20  dB 
across  the  subject  mm  wave  bands. 


1  .  INTRODUCTION 

TECOM  Industries  has  been  a  leading  supplier  of  wideband  direction  finding 
systems,  up  to  8  GHz,  consisting  of  a  shaped  ELPAR  reflector  and  feed,  a 
rotator  and  a  controller.  (The  ELPAR  reflector  has  an  ELliptical  shape  in 
elevation  and  a  PARabolic  shape  in  azimuth,  producing  a  vertically  oriented 
fan  beam.)  The  elliptical  contour  results  in  a  defocused  elevation  beamwidth 
whose  asymtotic  value  at  the  high  end  of  the  band  can  be  adjusted  to  any 


*  The  subject  of  U.S.  Patent  No.  4,831,384 
dated  May  16,  1989 
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convenient  values ,  typically  15°  or  35°,  to  receive  signals  from  both  airborne 
and  ground  emitters.  With  a  conventional  parabolic  contour  in  azimuth,  the 
energy  is  focused  for  increased  gain  and  a  narrow  beamwidth  so  as  to  permit 
precise  determination  of  line-of-bearing  to  the  emitter.  Using  multi-octave 
feeds*  such  as  log  periodic  dipole  arrays  or  cavity  backed  spirals,  having 
typical  bandwidths  of  1.0  to  18.0  GHz,  the  energy  from  the  spinning 
reflector/feed  antenna  subsystem  is  routed  through  a  dc  to  18.0  GHz  coaxial 
rotary  joint  to  the  base  of  the  pedestal. 

To  add  the  mm  wave  bands  of  18-26  and  26-40  GHz  to  such  a  spinning  DF 
system  we  are  faced  with  the  non-availability  of  critical  components. 

Although  single  channel  coaxial  rotary  joints  operating  from  dc  to  40  GHz 
exist,  there  are  no  dual  channel  joints  and  no  broadband  feeds  covering 
1.0-40  GHz.  Any  proposed  configuration  using  separate  antenna  bands,  such  as 
1-18,  18-26,  and  26-40  GHz,  in  order  to  circumvent  the  lack  of  a  full  range 
antenna,  requires  that  the  outputs  either  be  multiplexed  in  frequency  or 
sequentially  sampled  in  time  due  to  the  restriction  of  a  single  channel  rotary 
joint. 

For  frequency  separation,  two  triplexers  are  required  —  one  before  and  one 
after  the  rotary  joint.  This  results  in  disadvantages  which  are  a  9-11  dB 
loss  at  crossover,  along  with  1-2  dB  in-band  insertion  loss  and  high  cascaded 
VSWR  ripple  effects. 

If,  on  the  other  hand,  the  outputs  from  the  three  bands  are  time  sampled 
through  a  single  channel  rotary  joint,  many  emitters  may  be  lost  due  to  such  a 
time  sharing  bottle  neck.  If  dual  polarized  antennas  are  employed  in  the  mm 
wave  bands  the  effects  of  time  sharing  further  reduces  the  viewing  time  for 
1-18.0  GHz  frequencies  due  to  such  sharing.  The  switching  configuration  also 
requires  that  slip  rings  be  added  to  furnish  power  to  the  coaxial  or  waveguide 
switches . 
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What  is  clearly  needed  is  a  configuration  that  embodies  a  means  of 
providing  the  DF  function,  including  polarization  determination,  in  the  mm 
region  that  does  not  require  the  use  of  rotary  joints  or  multiplexers  to 
separate  the  18-26  and  26-40  GHz  bands.  Use  of  an  all  waveguide  system,  where 
feasible,  would  greatly  reduce  losses  and  improve  the  amplitude  of  signals 
delivered  to  the  receiver. 

2.  SYSTEM  DESIGN 

In  1986  the  writer  designed  and  had  built  such  a  system  for  resolving 
emitter  polarization  in  the  mm  wave  bands  of  18-26  and  26-40  GHz  by  piggy-back 
of  the  mm  wave  system  on  an  ELPAR  reflector/feed  system  operating  from 
1-18  GHz. 

Figure  1  shows  the  principal  components  forming  the  polarization-sensitive 
microwave  radiation  receiver.  A  horizontal,  circular  filter  disk  is  mounted 
for  rotation  in  the  sense  of  arrow  "R",  about  a  vertical  axis  of  rotation 
A-A* .  The  disk  comprises  a  thin  membrane  material  that  is  substantially 
transparent  to  radiation  in  the  microwave  frequency  band.  (It  should  be  noted 
that  the  terms  "vertical"  and  "horizontal"  in  this  context  need  not  refer  to 
gravitationally  defined  directions  —  merely  to  orthogonally  oriented  axes.) 

The  filter  disk  is  subdivided  into  3  concentric  regions.  The  outer  2  are 
defined  by  a  plurality  of  parallel  conductors  attached  to  the  disk,  and  the 
central  region  is  clear.  The  sets  of  parallel  conductors  in  the  outer  region 
are  orthogonal  to  the  array  of  similar  parallel  conductors  in  the  inner 
region.  Microwave  receiver  horns  are  statically  positioned  above  and  close  to 
the  surface  of  the  3  concentric  regions  of  the  disc. 

Also  mounted  for  rotation  on  the  axis  A-A*  is  a  flat  mirror  so  constructed 
that  it  acts  as  an  electro-optical  reflector.  In  Figure  1,  the  mirror  is  a 
planar  ellipse,  mounted  at  substantially  45°  with  respect  to  a  horizontal 
plane  orthogonal  to  the  axis  of  rotation  A-A* . 


FIGURE  1.  CONCENTRIC  DISK  DESIGN  WITH  GRIDDED  AREAS,  RECEIVING  HORNS,  AND  45‘ 
MIRROR 
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The  mirror  is  synchronously  driven  with  the  filter  disk,  so  that  a 
horizontally  aligned  minor  axis  B-Bf  of  the  mirror  remains  parallel,  at  all 
angular  positions  of  rotation,  with  the  array  of  parallel  conductors  defining 
the  inner  filter  region,  and  an  array  of  conductors  in  the  outer  filter 
region,  aligned  with  vertical  elements  in  the  mirror,  including  major  axis 
C-C*  therein. 


The  mirror  is  shown  in  Figure  1  as  lying  below  the  filter  disk  and  being 
inclined  at  substantially  45°  so  as  to  reflect  a  horizontally  incident  beam  of 
radiation  on  its  surface  vertically  upward  through  the  filter  disk. 

As  radiation  originating  on  the  horizon  is  reflected  through  the  filter 
-..sk,  vertically  polarized  components  therein  will  tend  to  pass  through  the 
conductor  array  of  the  inner  region  with  little  attenuation,  and  horizontally 
polarized  components  of  the  incident  microwave  beam  will  readily  pass  through 
the  outer  region.  The  converse,  in  both  regions,  is  not  so  and  radiation  with 
electrical  field  vectors  parallel  to  the  array  of  conductors  will  be  greatly 
reflected/attenuated.  Consequently,  the  microwave  receiving  sensors 
statically  positioned  above  the  filter  regions  will  be  exposed  to  two 
orthogonally  polarized  components  of  the  incoming  beam. 

Measuring  the  signal  strengths  collected  by  these  receivers  and  comparing 
them  will  clearly  indicate  the  relative  polarization  of  the  received  signal. 
(The  ratio  of  the  signal  strengths  to  which  the  receiver  horns  are  exposed  is 
independent  of  any  angular  differential  between  them.)  They  are  exposed  to 
radiation  which,  although  derived  from  a  signal  with  a  particular  polarization 
plane,  varies  in  direction  of  polarization  as  a  result  of  the  rotation  of  the 
mirror  and  of  the  filter  disk.  Consequently  it  is  important  that  these 


sensors  be  insensitive  to  such  variation, 


and  be  selected  from  the  class  of 


microwave  receivers  (e.g,,  antennas)  having  circular  polarization  and  low 


axial  ratios. 
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The  central  region  of  the  disk,  with  no  obstacle  to  the  transmission  of  any 
component  in  the  signal  reflected  from  the  mirror,  and  the  two  outer  regions, 
defined  by  arrays  of  orthogonal  grid  lines,  along  with  the  third  microwave 
receiver  sensor  complements  the  filter,  permitting  the  measurement  of  the 
total  signal  emanating  from  the  source  being  analyzed. 

Figure  2  shows  the  piggy* back  dual  range  MM  wave  system,  atop  a 
conventional  ELPAR  antenna  reflector,  and  a  planar  mirror.  The  antenna  and 
the  mirror  are  joined  by  supports  and  both  rotate  about  vertical  axis  A-A*  . 

The  mirror  is  secured  to  the  rear  of  the  ELPAR  reflector  and  is  exposed  to  any 
incoming  microwave  beam  during  a  complete  rotation  of  the  sensor  assembly 
around  the  axis  A  A* .  Rotational  movement  for  the  ELPAR  antenna  and  mirror  is 
derived  from  an  azimuth  rotator,  which  impels  the  turntable  onto  which  the 
antenna,  mirror,  and  antenna  feed  are  mounted. 

Incident  electromagnetic  radiation  striking  the  focusing  surface  of  the 
ELPAR  reflector  is  collected  by  the  feed,  located  at  the  focus  of  the 
parabolic  curve  defining  the  horizontal  shape  of  the  surface.  The  collected 
electromagnetic  signal  is  transferred  to  an  output  port  via  a  coaxial 
transmission  line,  with  a  rotating  coaxial  joint  installed  into  the  bearing 
supporting  the  turntable. 

Due  to  the  synchronization  in  the  rotational  motions  of  the  mirror  and  of 
the  filter  disk,  the  relative  strengths  of  the  three  signals  derived  from  the 
three  receiver  horns  do  not  vary  due  to  polarization  as  the  rotator  impels  the 
mirror  and  shaped  reflector  into  its  circulatory  motion.  The  strength  of  the 
signal  will  be  greatest  when  the  mirror  axis  B-B*  is  normal  to  the  point  of 
origin  of  the  incoming  microwave  beam. 

The  microwave  energy  captured  by  the  receivers  is  suitably  directed  toward 
analytical  instrumentation  *  -  capable  of  determining  the  relative  strengths  of 
the  three  signals. 
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FIGURE  2.  ELEVATION  VIEW  OF  POLARIZATION  DISK,  2  BAND  FEED  HORNS,  45°  MIRROR 
AND  ELPAR  REFLECTOR  AND  FEED 
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Figure  3  shows  in  detail  the  filter  disk,  the  adjacent,  concentric  regions 
formed  by  parallel  conductors  printed  onto  the  surface  of  disk  in  mutually 
orthogonal  arrays  and  the  6  receiving  horns  behind  the  three  zones  of  the 
disk.  It  is  enclosed  in  a  stationary  housing  constructed  from  a  dielectric 
material  that  forms  the  top  of  the  cylindrical  radome. 

One  pair  of  microwave  receiver  horns  is  mounted  onto  the  inner  surface  of 
the  cover  and  project  downward  through  the  clear  center  region  of  the  disk  and 
face,  the  reflecting  surface  of  the  mirror.  The  two  microwave  horns  are 
circularly  polarized  and  are  equally  sensitive  to  microwave  radiation  in  all 
polarized  states,  except  for  opposite- sense  circular  polarization.  Suitably, 
the  smaller  horn  is  sensitive  to  the  26-40  GHz  band,  while  the  second  horn  is 
sensitive  in  the  18-26  GHz  band. 

The  outer  region  of  the  filter  disk  is  surmounted  by  2  horns  (one  each  for 
18  to  26  and  26  to  40  GHz) ,  while  the  intermediate  region  is  surmounted  by  an 
identical  pair  of  horns .  These  horn  pairs  are  best  seen  in  Figure  3 .  These 
horns  are  also  supported  by  the  cover  and  are  positioned  with  the  lips  of 
their  sensing  horns  proximate  to,  but  not  touching,  the  moving  surface  of  the 
filter  disk.  The  energy  derived  from  an  incident  microwave  beam  reflected  by 
the  mirror  into  the  horns  is  conveyed  through  a  network  of  wave  guides  and 
cavity  switches  to  any  desired  measuring  equipment,  such  as  a  mm  wave  receiver 
subsystem,  for  further  processing. 

3.  PRODUCTION  HARDWARE  / 

A  series  of  over  13  DF  sets  were  fabricated  and  tested  using  the  polarizing 
discrimination  scheme  described,  for  the  mm  bands  of  18-26  and  26-40  GHz. 

Figure  4  is  a  photo  taken  looking  down  on  the  ELPAR  and  planar  spiral  feed 
to  the  back  of  which  is  mounted  the  mm  wave  mirror.  Atop  the  assembly  and 
perpendicular  to  the  axis  of  rotation  is  the  3-zoned  disk,  the  outer  two  zones 
of  which  consist  of  orthogonal  sets  of  copper  grids  etched  on  a  thin 
dielectric  membrane.  (Compare  the  photo  with  previous  Figures  2  and  3.) 
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NOTE:  WAVEGUIDE  SWITCHES  ALL  SET  FOR 
VERTICAL  POLARIZATION  BOTH  BANDS. 


FIGURE  3.  PLAN  VIEW  OF  FILTER  DISK  AND  2  BANDS  OF  RECEIVING  HORNS  WITH 
WAVEGUIDE  SWITCHING  AMONG  3  ZONES 
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Figure  5  is  a  photo  of  the  3  pairs  of  mm  band  horns  wherein  one  pair  (i.e., 
18  to  26  and  26-40  GHz)  are  located  in  each  of  the  3  zones.  Also  note  the 
waveguide  switches  and  rigid  waveguide  transmission  lines.  In  the  production 
configuration,  the  base  plate  in  »luding  all  the  waveguides  and  switches,  less 
the  apertures  of  the  3  pairs  of  horns,  are  covered  with  microwave  absorber  to 
suppress  reflections  from  the  mirror  to  mounting  base.  Such  reflections  can 
be  modulated  by  the  spinning  mirror  and  ELPAR  structures  resulting  in  DF  and 
processing  degradations. 

Figure  6  is  an  exterior  view  of  the  radome  enclosed  DF  system  operating 
from  1-40  GHz,  with  polarization  sensing  from  18-26  and  26-40  GHz.  Note  the 
positioning  of  the  biconical  omnidirectional  intercept  antennas  mounted  atop 
the  radome.  These  TECOM  standard  products  cover  the  same  1-18,  18-26  and 
26-40  GHz  bands  as  the  DF  system. 

4.  ELPAR  PERFORMANCE,  1-18  GHz 
The  basic  antenna  subsystem  consists  of  a  shaped  ELPAR  reflector  with  a 
planar  spiral  feed  covering  the  1-18  GHz  frequency  range.  Table  1  summarizes 
the  performance  measured  for  this  4.1  octave  antenna. 

Notice  that  this  ELPAR  reflector  has  an  elevation  beamwidth  which  converges 
to  about  15°  at  the  top  of  the  band.  The  fact  that  the  minimum  value  is  13.5° 
is  accounted  for  by  the  original  15°  design  being  based  on  a  slightly  reduced 
reflector  height  and  log  periodic  dipole  array  feed. 

5.  POLAR i ZAT I  ON  SENSING  PERFORMANCE,  18-26  AND  26-40  GHz 
Polarization  sensing  was  measured  extensively  on  the  first  three  systems. 
Standard  testing  of  VSWR,  patterns,  and  gain  were  supplemented  with  swept 
ratios  of  co-polarization  vs.  cross-polarization  outputs. 


FIGURE  5.  VIEW  OF  3  PAIRS  (18-26  AND  26-40  GHz)  OF  HORNS,  EACH  PAIR  LOCATED 
WITHIN  ONE  OF  3  DISK  ZONES 
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FIGURE  6.  EXTERIOR  VIEW  OF  RADOME  ENCLOSED  ELPAR  AND  MM  DF/POLARIZATION 
SENSING.  NOTE  BICONICAL  INTERCEPT  ANTENNAS  ATOP  THE  RADOME. 
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TABLE  1 


SUMMARY  OF  15°  ELPAR  PERFORMANCE 
1  TO  18  GHz 


Frequency ,  GHz : 

1.0 

2.0 

4.0 

8.0 

12.0 

15.0 

18^_0 

Gain,  dBic:  (1)  (2) 

10.0 

13.0 

16.0 

19.0 

21.0 

22.0 

23.5 

AZ  Beamwidth,  Deg.: 

23.5 

12.0 

6.1 

3.3 

2.1 

1.6 

1.4 

EL  Beamwidth,  Deg.: 

24.0 

15.7 

15.2 

15.0 

17.0 

13.0 

13.5 

Squint,  Deg.: 

.4 

1.3 

.5 

.1 

.1 

.15 

.15 

Sidelobes: 

12 

15 

16 

16 

14 

13 

14 

NOTES : 


(1)  Gain  measured  at  feed. 

(2)  Polarization  RHCP 

(3)  VSWR  2.0:1  maximum 

(4)  Reflector  size  27"  wide  x  30"  high 

(5)  Weight  approximately  20  pounds 


Figure  7  is  a  polar  plot  taken  at  21  GHz,  near  the  center  of  the  lower  nun 


band.  On  this  single  figure  are  plotted  three  separate  patterns,  with  the 
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The  three  patterns  are: 

•  Patterns  generated  for  DF  through  the  clear  zone  of  the  polarization  disk. 

•  Patterns  generated  from  incident  vertical  polarization  and  the  horizontal 
gridded  annulus. 

•  Patterns  generated  from  incident  horizontal  polarization  and  the  vertical 
gridded  annulus. 

To  avoid  confusion  between  the  three  separate  patterns  the  measurement  of 
the  DF  pattern  using  the  ungridded  zone  is  shown  oriented  at  180°  from  the 
other  two  patterns.  Note  that  the  gain  of  this  channel  is  greater  than  the 
gains  for  the  two  polarization  sensing  channels.  Since  the  outputs  of  the  two 
polarization  sensing  channels  are  shown  aligned,  as  in  normal  operation,  the 
ratio  of  desired  to  cross-polarization  is  easily  seen. 

Figure  8,  similar  to  Figure  7 ,  is  a  polar  plot  taken  at  32  GHz,  the  center 
of  the  upper  mm  band.  This  plot  also  includes  three  separate  patterns  for  no 
grid,  vertical  polarization,  and  horizontal  polarization.  As  previously,  the 
measurements  were  through  the  radome. 

Figures  9  and  10  are  swept  plots  from  18-26  GHz  of  the  received  levels  for 
the  two  polarization  outputs  (i.e.,  vertical  and  horizontal)  when  the  "mirror" 
is  rotated  to  peak  the  beam  in  the  direction  of  the  emitting  source.  Figure  9 
is  for  a  transmitted  polarization  of  vertical  while  Figure  10  is  for  a 
transmitted  polarization  of  horizontal.  These  figures  show  the  copolar  signal 
to  be  12-20  dB  greater  than  the  cross-polar  signal. 


FIGURE  9.  SWEEP  FROM  18.0  -  26.5  GHz  USING  VERTICAL  POLARIZATION  VS.  OUTPUTS 
FOR  VERTICAL  AND  HORIZONTAL 


FIGURE  1' 


Figures  11  and  12  are  similar  swept  plots  from  26-40  GHz  of  the  received 
levels  for  the  signals  filtered  by  the  polarization  disk  to  distinguish 
between  vertical  and  horizontally-polarized  emitters.  Figure  11  is  for  a 
verticaliy-polarized  source  while  Figure  12  is  for  a  horizontally-polarized 
source.  Clearly,  the  design  objective  of  discriminating  between  principally 
polarized  signals  sufficiently  to  identify  emitter  polarization  has  been  more 
than  adequately  met. 
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FIGURE  11.  SWEEP  FROM  26.5  -  40.0  GHz  USING  VERTICAL  POLARIZATION  VS.  OUTPUTS 
FOR  VERTICAL  AND  HORIZONTAL 


FIGURE  12.  SWEEP  FROM  26.5  -  40.0  GHz  USING  HORIZONTAL  POLARIZATION  VS 
OUTPUTS  FOR  VERTICAL  AND  HORIZONTAL 
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*  Abstract 

In  this  paper,  a  rigorous  analysis  is  performed  to  analyze  the  effect  of  the 
feed  on  the  input  impedance  of  a  circular  and  annular-ring-loaded  circular 
microstrip  antenna. 

The  traditional  approach  to  calculating  the  input  impedance  of  a  mi¬ 
crostrip  antenna  is  via  the  reaction  formula  (also  known  as  the  induced  EMP 
formula).  However,  this  formula  does  not  include  the  effect  of  feed  geometry 
on  the  input  impedance.  For  example,  for  a  probe-fed  microstrip  antenna, 
the  effect  of  the  gap  capacitance  at  the  base  of  the  probe  is  not  included  in 
the  reaction  formula. 

A  new  variational  expression  for  the  input  impedance  of  an  antenna  driven 
by  a  voltage  across  a  gap  is  derived.  This  formula  incorporates  the  effect  of 
the  gap  dimension  (which  affects  the  gap  capacitance)  on  the  input  impedance 
of  the  antenna.  Such  a  formula  is  used  to  derive  the  input  impedance  of  a 
probe-fed  microstrip  antenna. 

In  addition,  the  solution  of  the  microstrip  antenna  driven  by  a  coax¬ 
ial  probe  is  obtained  using  the  mode-matching  method,  by  matching  the 
modes  inside  the  coaxial  line  with  the  modes  of  the  microstrip  antenna.  The 
mode  matching  solution  is  done  with  and  without  a  diaphragm  at  the  base 
of  the  probe.  If  a  single  mode  is  used  in  the  mode  matching  approach,  the 
impedance  calculated  with  and  without  the  diaphragm  is  the  same  as  the 
input  impedance  obtained  by  the  varational  formula.  However,  the  mode¬ 
matching  solution  is  more  rigorous  than  the  variational  formula,  because,  in 
theory,  the  input  impedance  of  the  probe-fed  antenna  can  be  calculated  to 
any  desire  accuracy  depending  on  computational  resources. 

It  is  found  that  the  effect  of  the  gap-capacitance  on  the  input  impedance 
is  large  when  the  thickness  of  the  substrate  is  large.  This  is  because  the 

*  characteristic  impedance  of  the  parallel  plate  “transmission  line”  is  large 
when  the  substrate  thickness  is  large.  Since  the  gap  capacitance  shunts  this 
characteristic  impedance,  its  effect  is  more  pronounced  for  thick  dielectric 
substrates.  Our  theory  is  found  to  compare  very  well  with  experimental 
measurement. 
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I.  Introduction 

In  recent  years,  much  progress  has  been  made  on  the  input  impedance  of 
the  microstrip  antennas  [1-4].  To  have  a  more  general  solution,  an  integral 
equation,  solvable  by  Galerkin’s  method,  has  to  be  employed  to  solve  the 
problem  [3,  4],  A  simple  way  is  to  consider  a  uniform  current  on  the  probe  to 
be  the  exciting  source,  with  which  the  integral  equation  is  then  formulated  for 
the  unknown  currents  on  the  disk.  The  currents  can  then  be  solved  for  using 
standard  techniques  like  Galerkin’s  method  [5].  Having  found  the  current 
distribution  on  the  disk,  the  reaction  formula*  can  be  used  to  calculate  the 
input  impedance  of  the  antenna.  This  formula  is  [6] 


Z:n  = 


(Er,J) 

P 


(1) 


where  (Ep,  J)  =  fEr-JdV  is  the  reaction  between  E t  and  J.  Ep  is  the  total 

v 

electric  field  generated  by  the  probe  current  and  the  disk  current,  and  J  is 
the  probe  current  plus  the  disk  current.  When  the  disk  current  Jp  is  solved 
for  via  the  Galerkin’s  method,  it  can  be  shown  easily  that  (E p,  J#)  =  0  since 
Galerkin’s  method  sets  this  reaction  to  zero.  Therefore,  (1)  can  alternatively 
be  expressed  as 


Z;n  =  — 


(Er,Jp) 

P 


(E0,Jp)  +  (Ep,JP) 

P 


(2) 


where  we  have  written  Ep  =  E#  +  E p  and  E^  is  the  field  produced  by  the 
disk  current  while  Ep  is  the  field  produced  by  the  probe  current.  In  this 
decomposition,  the  physical  meaning  of  (2)  is  clear:  the  first  term  in  (2)  is 
the  input  impedance  of  due  to  the  disk  current,  whereas  the  second  term  in 
(2)  is  due  to  the  self  inductance  of  the  probe. 

One  thing  is  peculiar  about  the  expression  (2):  it  is  independent  of  the  size 
of  the  gap  at  the  base  of  the  probe.  This  gap  contributes  to  a  capacitance 

*  This  formula  has  also  been  called  the  induced  E.M.F.  formula,  while  the  power  formula 
Zin  =  -(Er,  J')/|J|2  has  also  been  called  the  induced  E.M.F.  formula.  We  call  this  the 
reaction  formula  to  avoid  confusion. 
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that  shunts  the  input  impedance  of  the  microstrip  antenna.  Hence,  it  is 
exceedingly  important  when  the  gap  at  the  base  of  the  probe  is  small,  or 
when  the  input  impedance  of  the  antenna  is  large.  A  small  gap  gives  rise  to  a 
large  capacitance  which  shunts  the  input  impedance  of  the  antenna;  and  the 
impedance  is  more  sensitive  to  this  shunt  when  it  is  large.  When  the  dielectric 
substrate  is  thick,  the  effect  of  this  gap  capacitance  is  more  important  as  we 
shall  show. 

In  this  paper,  we  shall  illustrate  three  ways  of  incorporating  this  gap 
capacitance  into  the  input  impedance  calculation  of  a  probe-fed  antenna. 
The  three  methods,  with  increasing  level  of  complexity  and  sophistication,  are 
the  variational  formula,  the  mode-matching  method,  and  the  mode-matching 
method  with  a  diaphragm. 

II.  Variational  Formula 

When  an  antenna  is  fed  by  a  coaxial  line  as  shown  in  Figure  la,  it  could 
be  replaced  by  an  equivalent  problem  of  Figure  lb  exterior  to  the  coaxial  line. 
In  Figure  lb,  M„  =  Ea  X  n,  i.e.,  the  tangential  component  of  the  electric  field 
at  the  aperture  of  the  coaxial  line  in  Figure  la.  Examples  of  such  antennas 
are  probe-fed  microstrip  antennas  and  coaxial  fed  monopoles. 

In  Figure  lb,  Ma  is  actually  unknown,  but  when  the  aperture  of  the 
coaxial  line  is  small,  its  approximate  distribution  can  be  ascertained.  With 
this  in  mind,  we  shall  assume  that  Ma  in  Figure  lb  is  known  which  is  also  the 
excitation  source  for  the  field.  Consequently,  an  electric  current  J4  is  induced 
on  the  antenna  structure  due  to  the  presence  of  M0.  A  new  variational 
expression  for  the  input  admittance  of  the  antenna  can  be  derived  such  that* 

-fm  —  y2  *  W 

In  the  above,  Hy  and  Ey  are  the  total  magnetic  field  and  electric  field  due  to 
Ma  and  J4.  V  is  the  voltage  across  the  gap  at  the  aperture  of  the  coaxial  line. 

t  It  is  correctly  pointed  out  to  the  first  author  by  R.  E.  Collin  [6]  that  such  a  formula 

does  not  exactly  reproduce  the  input  admittance. 
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Furthermore,  H t  =  Hs  +  Ha,  Ey  =  E3  +  E0,  where  the  subscript  s  denotes 
fields  due  to  J3,  and  the  subscript  a  denotes  fields  due  to  Ma.  Taking  the 
first  variation  of  (3)  by  keeping  Ma  and  V  fixed,  we  have 

VHYin  =  -(M0,mr)  +  (SEt,  J4e)  +  (Er.,«,),  (4) 

where  the  extra  subscript  e  denotes  exact  values.  Since  ii  x  Ey,.  =  0  on  the 
antenna  structure,  the  last  term  in  (4)  vanishes.  Furthermore,  since  only  Ja 
is  allowed  to  vary  but  not  Ma,  (JEy  —  £E.,.  Consequently, 

(5) 

implying  the  vanishing  of  the  first  variation  of  8Y{n  in  (4)  and  the  stationarity 
of  (3 ;  about  the  exact  solution.  It  can  be  shown  that  Equation  (3)  is  also 
stationary  if  Ma  and  V  are  allowed  to  vary. 

When  the  exact  current  J3C  is  known  on  the  antenna,  then  (Jje,  Eje)  =  0 
since  n  x  Ey,.  =  0  on  the  antenna  structure  and  (3)  becomes 

(Ma,HTe)  _  -(Mq,Hq)  -  (Mq,Hje) 

^,nc  yi  y2 

The  first  term  on  the  right-hand  side  of  (6)  is  identifiable  as  the  contribution 
to  the  input  admittance  due  to  the  gap  capacitance,  while  the  second  term 
is  the  contribution  due  to  the  antenna  structure. 

The  numerator  of  (3)  is  the  reaction  between  the  sources  (Ma,  J3)  with 

the  total  field  (Hr,  E^)  generated  by  the  sources.  Hence,  it  falls  into  the 
class  of  variational  formulas  derivable  from  reaction  [7].  When  an  antenna 
is  driven  by  a  voltage  gap  (slice-generator)  whose  equivalent  problem  is  an 
antenna  driven  by  a  magnetic  current,  (3)  is  the  variational  formula  to  be 
used  but  not  (1)  as  has  been  erroneously  suggested  by  some  authors  [8]. 

Before  using  (3),  the  value  of  J3  is  required.  One  way  of  obtaining  J4  is  to 
formulate  an  integral  equation  from  which  J3  could  be  solved  by  Galerkin’s 
or  moment  method.  The  integral  equation  is  obtained  via  the  boundary 
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condition  that  n  x  (E„  4-  E4)  =  0  on  S'  where  S  is  the  surface  of  the  antenna 
structure.  E4  is  generated  by  J4  via  a  dyadic  Green’s  function.  Consequently, 
the  integral  equation  is 

— n  x  E0(r)  =  iu  J  dr'h  x  Ge(r,r')  •  //(r')J4(r'),  r  e  S  (7) 


where  Ge(r,  r')  is  the  electric-type  dyadic  Green’s  function  for  a  general  in¬ 
homogeneous  medium  [8].  To  solve  the  above  via  the  Galerkin’s  method,  we 
let 

N 

J,  =  £>nJ,n  (8) 

n=l 

in  (7)  and  testing  the  resultant  equation  with  n  x  J4m,  m  =  1, •  •  •  ,N.  As  a 
result,  we  have 


N 

(Jann^o)  =  W  ^  m  =  1,  •  •  •  ,  AT.  (9) 

n=l 

where  (Jsm,Ge/z,  J4m>  =  J  dr /dr'J4m(r)  •  Ge(r,  r')/x(r')  •  Js„(r').  Equation 
(9)  is  a  matrix  equation  from  which  an  can  be  solved  for  yielding  J4  from  (8). 

Multiplying  (9)  by  am  and  summing  over  m  from  1  to  N ,  we  can  easily 
show  that  J4  from  (8)  satisfies 

=  fu>(J4,Ge/i,  J4)  =  (J4,E4).  (10) 

Therefore,  J.  obtained  from  Galerkin’s  method  will  make  the  second  term  in 
the  numerator  of  (3)  vanish,  and 

(Ms;Hr)  (Ma,Ha)  +  (Ma,H4) 

Jtn  y2  y2  ' 

We  have  applied  the  above  technique  to  solve  for  the  input  impedance 
of  a  probe-fed,  annular-ring-loaded  (ARL),  circular  microstrip  disk  antenna 
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as  shown  in  Figure  2  with  the  diaphragm  absent.  Due  to  the  presence  of  a 
dielectric  substrate,  a  stratified  medium  Green’s  function  has  to  be  used  to 
calculate  the  field  due  to  the  current  induced  on  the  probe  and  the  disk.  The 
field  distribution  E0  on  the  coaxial  aperture  is  assumed  to  be  that  of  the 
TEM  mode  of  the  coaxial  cable.  The  resultant  input  impedance  calculation 
using  (3)  is  compared  with  previous  calculation  using  (1)  in  Figure  3.  Better 
agreement  with  experiment  is  noted  when  Equation  (3)  is  used.  We  have  not 
used  any  fudge  factor  in  the  theory.  The  integral  equation  is  solved  similar 
to  previous  work  using  vector  Hankel  transforms  [10, 11]  and  the  singularity- 
subtraction  method  [3,  12]  to  accelerate  convergence. 

III.  Mode  Matching  Method* 

In  using  the  variational  formula  (3),  one  has  to  assume  knowledge  of  Ma. 
The  formula  itself  does  not  lend  much  insight  as  to  how  Ma  is  to  be  chosen. 
In  order  to  systematically  solve  for  the  value  of  M„,  the  mode-matching 
method  is  needed.  In  th.  method,  the  reflected  field  inside  the  coaxial  cable 
is  expanded  in  terms  of  the  modes  of  the  cable.  Then,  continuity  of  tangential 
fields  at  the  aperture  is  imposed  to  solve  for  the  reflection  coefficients  of  the 
reflected  modes.  For  long  wavelengths,  only  the  TEM  reflected  mode  will 
be  propagatory  while  the  higher-order  modes  are  evanescent.  The  reflection 
coefficient  of  the  TEM  mode  allows  one  to  unambiguously  derive  the  input 
impedance  of  the  antenna.  We  shall  consider  the  case  where  the  diaphragm 
is  absent. 

Inside  the  coaxial  region,  one  assumes  that 


L-l 


E(r)  =  E0(r3)e{k*  +  £ T,E,(r3)e-{k 
/=0 


(12) 


where  Eo(r3)  is  the  field  of  the  TEM  mode  and  Em(r,),  m  >  0,  are  the  fields 


t  The  mode-matching  method  has  been  used  previous  in  conjunction  with  tie  cavity¬ 
mode  model  [13]. 


of  the  higher-order  modes.  The  corresponding  magnetic  field  is 


H(r)  =  H0(r,)e"“  -  £  r,H,(r,)e-*-. 


Hence,  the  total  electric  field  at  the  aperture  in  accordance  with  (12)  is 


E0  =  e0  +  y]rtEt. 


The  corresponding  magnetic  current  is 


M0  =  Ea  x  n  =  M0  +  ^2 


where  M/  =  E/  x  n.  The  electric  field  in  the  outer  region  generated  by  Ma 
can  be  written  as 

E .=  [  *'«-■(  r)V  x  Gm(r,  r')e(r')  •  M. 

V 

=  f  dr'e~v V  x  Gm(r,  r')e  •  M0  +  Y  T,  [  dr'e"1  V  x  Gm(r,  r')e  •  M, 
l  feo  l  (16) 


where  Gm(r,  r')  is  the  magnetic-type  dyadic  Green’s  function  in  an  inhomo¬ 
geneous  medium  due  to  a  magnetic  current  source.  Using  (16)  in  (9),  we 
have 

L- 1 

-  (J4m,e_1V  x  Gme,M0)  -  x  Gme,  M/)T/ 

/=0 

N 

=  Gc///5  J sn)ani  txi  =  1,  *  •  •  ,  N.  (17) 

n= 1 

The  magnetic  field  outside  the  aperture  is  Hr  =  Ha  +  Hs  where 

H-  =  iu)  J dr'Gm(r,  r')e(r')  •  M„(r'),  (18a) 


H a  =  J dr'/T1  V  x  Ge(r,  r')n  •  J,(r'). 


(18b) 
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With  M0  from  (15)  and  J,  from  (18),  the  above  becomes 
Hy  =  iu  J  dr'Gm(r,  r')e(r')  •  M0(r')  +  iioTt  J  dr'Gm( r,  r')e(r')  •  M,(r') 

V  1=1  1/ 


N  r  ~~ 

+  /  dr'fx~1^  x  Gc(ry)fi  •  Ja»(r'). 

n=l  1 


(19) 


Continuity  of  the  tangential  components  of  the  magnetic  fields  requires  the 
equality  of  the  tangential  components  of  (13)  and  (19).  Equating  (13)  and 
(19),  and  testing  the  result  with  Mp(ra),  p  =  0,  •  •  •  ,  L- 1,  from  (15),  we  have 

L-l  _  L-l 

Mo  -  =  MMP,Gme,M0)  +tw]Tr/<Mp,Gme,M,) 

/=o  1=0 

N 

4*  'y  ] an(Mp,  p  V  x  Ge/z,  Jin)>  P  =  0,  •  •  •  , L  —  1,  (20) 

n=l 


where  A/  =  (M/,  H/)  is  a  normalization  constant  and  mode  orthogonality  [14] 
has  been  used  in  arriving  at  (20).  Equations  (17)  and  (20)  form  a  set  of 
N  +  L  linear  algebraic  equations  from  which  T/,  /  -  -  0,  •  ■  •  ,  L  —  1,  and  an, 
n  —  1,“*  ,N  can  be  found.  Once  r0  is  known,  v  input  admittance  of  the 
antenna  is  found  from 

r<n = (21) 

t  +  j.  0 

where  Yq  is  the  characteristic  admittance  of  the  coaxial  line. 

If  only  a  single  TEM  mode  is  used  to  approximate  the  field  inside  the 
coaxial  cable,  one  can  readily  show  that  the  input  admittance  obtained  via 
(21)  is  identical  to  that  using  the  variational  formula  (3)  where  a  TEM  field 
is  assumed  for  Ea  at  the  aperture.  Hence,  the  variational  formula  is  a  special 
case  of  the  mode-matching  method  under  the  single-mode  approximation. 

We  have  implemented  the  above  formulation  for  the  annular-ring-loaded 
microstrip  antenna  using  a  stratified  medium  Green’s  function.  Figure  4 
shows  the  comparison  for  the  ARL  microstrip  disk  antenna  between  the  mode 
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matching  theory  and  experiment  for  different  V s  (L  =  1,  being  the  single¬ 
mode  approximation,  is  equivalent  to  the  variational  formula).  It  is  seen  that 
the  need  for  higher-order  mode  is  not  that  critical  in  this  case. 

IV.  Mode  Matching  with  a  Diaphragm 

Sometimes,  a  coaxial-fed  antenna  is  made  with  a  diaphragm  at  the  aper¬ 
ture  as  shown  in  Figure  2.  It  will  be  interesting  to  study  the  effect  of  this 
diaphragm  on  the  gap  capacitance,  and  hence,  the  input  impedance  of  an 
antenna.  This  method  proceeds  very  much  like  the  mode-matching  method 
except  that  the  aperture  field  now  is  expanded  in  terms  of 


<3-1 

E  „  =  (22) 

9=0 

where  iy9  is  a  basis  function  with  support  only  over  the  aperture.  Hence, 


0-i 

Ma  =  E0xn  =  ^69M,  (23) 

9=0 

where  M,  =  i |,xn,  The  equivalence  of  (17)  is  then 


0-i  N 

x  Gme,M,)6,  =ik>Bjam,Ge/b  J™K,  m  =  l,--  ,N. 

q-0  n=l 


(24) 

Again,  the  fields  inside  the  coaxial  line  are  expanded  as  (12)  and  (13).  We 
equate  (15)  with  (23)  to  establish  the  continuity  of  the  electric  field,  and  test 
the  resulting  equation  with  Up,  p  —  0,  •  •  •  ,  L  —  1  to  yield 


L- 1  O-i 

Vo  +  Yj r,6'»A'  =  Y  M1^)’  P  =  0,  •  •  *  ,  L  -  1.  (25) 

1=0  9=0 
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The  equivalence  of  (19)  is 


Q~l  r  _ 

Hr  =  ita  Ylbl>  dr'Gm(r, r')e(r')  •  M,(r') 

9=0  * 


+ 


(26) 


Equating  (26)  with  the  magnetic  field  inside  the  waveguide  given  by  (13), 
and  testing  the  resultant  equation  with  Mr,  r  =  0,  •  •  •  ,  Q  —  1,  we  have 
L- 1  <5-i 

(Mr,H„)  -£r,(MT)H,)  =;W£6,(Mr,G„e,M,) 

1=0  7=0 

N 

+  J]«„(Mri/j'1VxGe/i)JJ„))  r  =  (),•••  ,Q  -  1.  (27) 

n=l 

Equations  (24),  (25),  and  (27)  constitute  N+L+Q  linear  algebraic  equations 
from  which  the  unknown  an,  &9,  and  T|  can  be  found.  Once  r0  is  known,  the 
input  admittance  is  given  by  (21). 

Strangely  enough,  if  a  single  mode  is  used  in  (12)  and  (13),  and  only 
one  basis  function  is  used  in  (22),  the  input  admittance  found  this  way  can 
again  be  shown  to  be  identical  to  that  derived  using  the  variational  formula. 
However,  the  use  of  a  single  mode  in  the  coax  will  not  meet  the  boundary 
condition  on  the  diaphragm,  and  the  result  is  inaccurate.  In  other  words,  the 
variational  formula  is  inaccurate  for  the  geometry  shown  in  Figure  2  where  a 
diaphragm  exists. 

This  theory  has  been  used  to  calculate  the  input  impedance  of  a  probe- 
fed  microstrip  antenna  with  a  diaphragm.  Figure  5  shows  the  comparison 
of  the  theory  with  a  diaphragm  and  experiment  for  an  ARL  microstrip  disk 
antenna.  Excellent  agreement  is  observed.  Figure  6a  shows  the  effect  of 
varying  the  diaphragm  size  for  a  thick  substrate  case,  while  Figure  6b  shows 
the  effect  of  varying  the  diaphragm  size  for  the  thin  .substrate  case.  Varying 
the  diaphragm  size  alters  the  gap  capacitance.  Hence,  the  effect  of  the  gap 


10 

606 


1 


l 


1 


capacitance  on  the  thick  substrate  case  is  more  pronounced  than  the  effect 
on  the  thin  substrate  case. 

The  integral  equations  are  solved  using  the  method  of  reference  12.  With 
the  notation  of  reference  12,  in  Figure  3,  (M,  P,  L)  Q )  —  (2, 1,4, 1)  in  curve  2, 
and  (M,  P,  L ,  Q)  =  (2, 1, 8, 1)  in  curve  3.  In  Figures  4,  5,  and  6,  (M,  P,  £,  Q) 
=  (1,0, 2,0)  foi  all  curves. 

V.  Conclusions 

Rigorous  calculations  of  a  probe-fed,  annular-ring-loaded,  circular-disk 
microstrip  antenna  is  presented.  The  first  calculation  uses  a  new  variational 
formula  which  improves  on  the  input  impedance  of  an  antenna  driven  by  a 
voltage  gap.  The  variational  formula  is  shown  to  be  good  for  the  case  when 
the  aperture  at  the  coaxial  line  does  not  have  a  diaphragm. 

The  second  calculation  uses  a  mode-matching  method  where  the  field  in 
the  coaxial  line  is  expanded  in  terms  of  the  waveguide  modes.  It  can  be  shown 
that  if  only  a  single  mode  is  used  to  expand  the  field  inside  the  waveguide, 
the  mode  matching  solution  reduces  to  that  of  the  variational  formula.  It  is 
found  that  when  the  wavelength  is  long  compared  to  the  transverse  dimension 
of  the  coaxial  line,  the  use  of  higher  order  modes  is  not  crucial  in  improving 
the  input  impedance  calculation.  1 

The  last  calculation  uses  a  mode-matching  method  with  a  diaphragm. 
Although  the  theory  under  the  single  mode  approximation  agrees  with  the 
variational  formula,  the  single  mode  approximation  is  bad  in  this  case  because 
the  boundary  condition  on  the  diaphragm  cannot  be  met.  Therefore,  one 
concludes  that  the  variational  formula  cannot  be  used  when  the  aperture 
consists  of  a  diaphragm.  This  is  consistent  with  the  fact  the  variational 
formula  (3)  need  not  exactly  represent  the  input  admittance. 


4- V>  0J7A  fn  olf  a**  IVa  CfriTi  f 
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the  input  impedance  and  the  resonant  frquency  of  the  microstrip  antenna  is 
more  sensitive  to  the  influence  of  this  shunting  gap  capacitance  when  the 
dielectric  substrate  is  thick  as  evidenced  by  Figure  6.  This  can  be  under- 
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stood  from  an  approximate  equivalent  circuit  of  the  microstrip  antenna  close 
to  resonance  as  shown  in  Figure  7  (we  ignore  the  probe  inductance  in  this 
discussion.).  When  the  dielectric  substrate  is  thick,  the  effective  C  in  the 
equivalent  circuit  is  smaller  due  to  the  larger  separation  between  the  top 
plate  and  the  bottom  plate.  Since  Cg  is  in  parallel  with  C ,  its  influence  on 
the  resonance  frequency  of  the  circuit  is  more  pronounced  as  shown  in  Figure 
6a. 
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Figure  1.  (a)  A  probe  fed  antenna,  (b)  The  equivalent  problem  exterior  to  the 
coaxial  line  where  Ma  =  E  x  n. 


Figure  2.  An  annular- ring-loaded  microstrip  antenna  fed  by  a  coaxial  probe.  In 
Figures  3,  4,  5,  and  Ga,  a  =  1.2cm,  b  =  2.4cm,  ax  ~  0.815cm,  ef  =  2.65, 
V  =  0.6cm,  R  a  0.0635cm,  ki  =  0.2cm,  d  =  0.3175cm. 
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Figure  8.  «m  thtwi}  using  ^1),  (3),  aim  «;a.|>ci iinctic,  vmm’H*  hi v<*  i 

experiment,  curve  2— theory  using  Equation  (3),  and  <urvc  3-  wiemy 
using  Equation  (1). 


Zin  (ohm) 


Figure  4.  Comparison  of  theory  using  mode-matching  method  and  experiment,  for 
an  ARl  microstrip  disk  antenna  where  the  curve  with  full-circle  symbols 
is  the  experiment,  and  the  curves  with  open-circle,  triangle,  and  square 
symbols  are  mode-matclung  theory  with  1  TEM  incident  mode,  and  1  3 
and  5  TM  reflected  modes,  respectively,  for  the  coaxial  line. 
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Figure  5.  Comparison  of  theory  using  mode  matching  will,  a  diaphragm  and  ex¬ 
periment  for  an  annular-ring-loaded  microstrip  antenna,  where  the  curve 
with  full  symbols  is  experiment  while  the  curve  with  open  symbols  is  the 
theory.  The  outer  radius  of  the  aperture  is  b  =  0.12cm. 
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Figure  fla.  Effect  of  diaphragm  size  on  the  input  impedance  of  an  ARL  microstrip 
disk  antenna  with  a  thick  substrate,  where  curves  1,  2,  3,  4,  and  5  are  the¬ 
oretical  results  with  substrate  thickness  d  =  0.1At  and  various  diaphragm 
sizes  where  the  outer  radii  of  the  aperture  arc  ht  =  0.10,0.12,0.14,0  1G, 
and  0.18cm,  respectively. 
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Figure  Ob.  Effect  of  diaphragm  size  on  the  input  impedance  of  a  circular-disk  mi¬ 
crostrip  antenna  with  thin  substrate,  where  curves  1,  2,  3,  4,  and  5 
are  theoretical  results  with  substrate  thickness  d  =  0.023 A,  and  var¬ 
ious  diophragm  sizes  where  the  outer  radii  of  the  aperture  are  b\  = 
0.10,0.12,0.14,0.16,  and  0.18cm,  repectively.  a  =  4.0cm,  b'  —  2.5cm, 
R  =  0.0G35cm,  &2  =  0.2cm,  d  =  0.3175cm  in  this  case. 
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Figure  7  An  approximate  equivalent  circuit  of  the  microstrip  antenna  to  explain 
the  more  severe  influence  of  the  gap  capacitance  Ct  on  the  resonant  fre¬ 
quency  when  the  substrate  is  thick  or  C  is  small. 
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Command,  Control,  Communications  and  Intelligence 
{ C 3 1 }  activities.  Technical  and  engineering 
support  within  areas  of  competence  Is  provided  to 
ESV  Program  Calces  ( POs )  and  other  ESV  elements 
to  perform  elective  acquisition  o f  C$1  systems. 

The  areas  of  technical  competence  Include 
communications,  command  and  control,  battle 
management,  information  processing,  surveillance 
sensors ,,  Intelligence  data  collection  and  handling, 
solid  state  sciences,  electromagnetics ,  and 
propagation,  and  electronic,  maintainability , 
and  compatibility. 


